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1 .0 INTRODUCTION 


The use of any new and unfamiliar tool is often accompanied by errors from 
misunderstanding or simply from lack of experience. The TETRA program is no 
exception. However, care has been taken to minimize these problems in the use 
of TETRA as a computational tool for engine dynamics. 

Based on the modal synthesis approach, the component element method 
employed in TETRA follows a modular or building block scheme both in the con- 
struction of the mathematical model of an engine to be analyzed and in the 
architecture of the program structure, subroutines and nomenclature. It 
cannot be overemphasized the importance in keeping this idea of modular 
construction in mind while preparing the schematic of the engine model and the 
inputs for TETRA. 

The turbine engine is described by a reduced system of second order 
differential equations and a solution for the transient response is obtained 
through an explicit numerical integration scheme - the central finite differ- 
ence method. Global stiffness and damping matrices are not assembled and only 
the right hand side of the system of equations, i.e*, the forces, is updated 
at each time step. The differential equations are formulated in terms of 
generalized coordinates and model elastic and rigid body elements and elements 
that represent physical connections. In the former case, the elements are 
obtained through the coordinate transformations that are associated with the 
free-free modes and partially constrained modes computed for engine subsystem 
structures (rotors, case, pylon). In the latter case, the elements describe 
the physical connections between the subsystem structures. These connections 
can be nonlinear and are defined with bearing/frame springs and dampers, 
engine support elements, link elements, rotor-case rub springs or stop ele- 
ments, and gyroscopic cross-axis coupling elements. The rotor-case rub 
springs are used to model the additional load path between the rotor and the 
case that exists for the large rotor excursions caused by high rotor unbalance. 
For this element, the effect of the force dead band associated with the 
structural clearance is included in the formulation for the effective restor- 
ing forces. 



One should begin with a schematic of the structure to be analyzed; this 
would probably have had its normal modes and frequencies calculated. Using 
a unified coordinate <^xis, one would identify the various structural subsys- 
tems into which the whole structure could be broken down, to be followed by 
establishing the connecting elements, their locations, subsystems joined and 
their mechanical characteristics (spring rates, damping). The normal modes of 
each structural subsystem may then be calculated in a program such as GE's 
VAST or finite element programs with eigenvalue capability such as SAP4, 
NASTRAN, etc. From these modal calculations, one chooses the modes and number 
of modes to represent each subsystem, and then generate the modal input 
file. 

Next, the engine operating conditions such as speed, amount and location 
of unbalance, time interval or external excitation are chosen. Finally, the 
joints where loads and displacements are to be printed out are defined. 

The assembly of the entire engine is made by the sequence and type of 
inputs and the NAMELIST input provides ease and flexibility in input prepara- 
tion. Following the input instructions closely and in the sequence in which 
these are presented in this manual, will minimize the errors and any confusion 
in using this new program. 
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2.0 DESCRIPTION OF TETRA 


The transient dynamic analysis method used in the TETRA (Turbine Engine 
Transient Response Analysis) program is based on a component element approach. 
The component elements consist of elastic and rigid body elements described by 
generalized coordinates obtained by coordinate transformations, and physical 
connecting elements that model bearing/frame springs and dampers, rotor-case 
rub springs and gyroscopic cross-axis coupling effects. 

The generalized coordinates are based on the free-free modes and par- 
tially constrained modes associated with engine subsystem structures (rotor, 
casing, pylon). The method extends the conventional modal analysis procedure 
to account for physical damping and symmetric stiffness terms, and rotor-case 
rubs including the effects of the force deadband associated with the struc- 
tural clearance. 

An efficient numerical time integration scheme - the central finite dif- 
ference method - is used to obtain the solution through an explicity step-by- 
step calculation. 

For each generalized coordinate, a differential equation is formed which 
relates the current generalized accelerations to the current generalized 
forces as follows: 

®i^i * Qi » i 2, 3, 4, (1) 

qi are the current generalized accelerations 

Qi are the current generalized forces 

m£ are the generalized masses 

Figure 1 shows a schematic of a subsystem and physical component element 

definition for an engine system. The nonlinearities can be treated with the 
physical connecting elements, and the subsystem natural modes are used to 
define regions of the engine system which are expected to remain linear. As 
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Figure 1. Typical Subsystem and Physical Component Representation of an Engine System, 



long as the nonlinearities are not severe enough to induce operating regimes 
where the original modes are totally unrepresentative, then the use of the 
natural modes is effective. 

The major TETRA subroutines where the transformations and the numerical 
integration of the generalized differential equations is performed are iden- 
tified in Figure 2. 

The task of establishing the current . total generalized forces acting 
on each of the generalized coordinates is accomplished in the subroutine 
GEN. The current physical displacements and velocities, needed to establish 
the physical coupling forces used by GEN in the formation of the generalized 
forces, are determined in subroutine CURRT which transforms the modal displace- 
ments qj[^ into real space at each time step. It will be noted that in comput- 
ing the generalized forces from the physical displacements and velocities, 
that a transformation from real space to generalized space is performed by 
subroutine GEN for each time step. The generalized forces also include the 
effects of applied physical forces which are not displacement or velocity 
dependent, such as unbalance forces. The generalized forces and the general- 
ized masses are utilized in subroutine TILOOP to compute the generalized 
accelerations using equation*(6) , and to compute the future values of the 
generalized displacements using the central finite difference method as shown 
in equation*(9) . 


*These equations are given as labeled in Section 3.0 Volume II of the "Blade Loss 
Turbine Engine Dynamic Analysis Program." 
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Loop Subroutine 
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Figure 2. 


TETRA Rou 
Numerical 
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The modal subsystem numbers, types, and the regions they represent are 
identified in Figure 3. The number of directions (physical degrees of freedom) 
and the direction ntxmbers associated with these degrees of freedom are listed 
in Figure 4. 

Eleven triple subscripted arrays are used to store the subsystem mode 
shape data. 

Sf (I, J, K) “ Mode shape modal displacements for the i-th 
modal subsystem. 

i « 1, 2, 3, . . . . , 11 (see Figure 3) 

I identifies the generalized coordinate 
J identifies the location or point on the mode shape 
K identifies the direction (see Figure 4) 

The S£ array data are used to define the coupling ratios, r, d needed 
to perform the transformations between real and generalized space. The 
coupling ratio r, d is the deflection at point r in the direction d for a 
unit value of the i-th generalized coordinate. 

The physical displacement and velocity data for the current time step 
are stored in the arrays X(I, J) and VEL(I, J), respectively. 

I = The point number 

J ^ The direction number (direction numbers 1-6) 

The generalized coordinate values are stored in the array Z(I, J). 

Z(I, J) Present and previous values for the generalized coordinates 
I = Generalized coordinate number 

J “ 1 for current time (0), 2 for one previous time step (-1), 3 for two 
previous time steps (-2). 

The generalized masses, stiffnesses and damping coefficients for the 
governing differential equations are stored in the arrays ZM, ZK, and ZC. 

ZM(I) ^ Generalized masses 

ZK(I) « Generalized stiffnesses 

ZC(I) ^ Generalized damping coefficients 
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Figure 3. Modal Subsystems. 
SUBSYSTEM 

VERTICAL PLANE ROTOR- 1 MODEL 
HORIZONTAL PLANE ROTOR- 1 MODEL 
3D RIGID BODY ROTOR-1 MODEL 
VERTICAL PLANE ROTOR-2 MODEL 
HORIZONTAL PLANE ROTOR-2 MODEL 
3D RIGID BODY ROTOR-2 MODEL 
VERTICAL PLANE CASE MODEL 
HORIZONTAL PLANE CASE MODEL 
3D RIGID BODY CASE MODEL 
TORSIONAL CASE MODEL 
3D FLEXIBLE PYLON (WING) MODEL 


TYPE 

FLEXIBLE-PLANAR 

FLEXIBLE-PLANAR 

RIGID WITH 5D0F 

FLEXIBLE-PLANAR 

FLEXIBLE-PLANAR 

RIGID WITH 5D0F 

FLEXIBLE-PLANAR 

FLEXIBLE-PLANAR 

RIGID WITH 6 DOF 

FLEXIBLE-TWIST 

FLEXIBLE-3D 



SUBSYSTEM MODE L 


SUBSYSTEM DESCR I PTION 

1 ROTOR- I VERT. PLANE 

FLEX. MODEL 

2 ROTOR- i HOR. PLANE 

FLEX. MODEL 

3 ROTOR- 1 RIGID BODY 

MODEL 


A ROTOR-2 VERT. PLANE 

FLEX. MODEL 

5 ROTOR-2 HOR. PLANE 

FLEX. MODEL 

6 ROTOR-2 RIGID BODY 

MODEL 


7 CASE VERT. PLANE 

FLEX. MODEL 

8 CASE HOP,. PLANE 

FLEX. MODEL 

9 CASE RIGID BODY 
MODEL 

10 CASE- TORSIONAL FLEX. MODEL 

!1 PYL0N-3D FLEX. MODEL 


(C 


Figure 4 


Global 



Degrees of Freedom. 


The physical coupling forces used in the formation of the generalized 
forces are computed by multiplying the physical displacements and velocities 
into the connecting element stiffness and damping matrices, respectively and 
then adding the results. The current values for the physical connecting 
element data are stored in the arrays AKE^ and AD££. 

AKE£(I,J,K) = Stiffness array for the i-th physical connecting element 
type. 

ADE£(I,J,K) = Damping array for the i-th physical connecting 
element type. 

i = 1,2, .... Identifies the physical connecting element type. 

For example, the generalized spring-damper element, 
the rub element, etc. 

I Identifies the physical connecting element number. 

J Identifies the row number associated with a force for 
a given point and direction. 

K Identifies the column number associated with a displacement for 
a given point and direction. ^ 

The information needed to identify the locations of the AKE and ADE 
array elements is provided in the following arrays. 

ICOMPJ (I) = Number of points for element I. 

ICOMPE (I,J) ® The global point numbers for element I associated 

with the J-th point . 

ICOMPD (I,J) = The number of directions for element I associated 

with the J-th point . 

ICOMPN (I,J,K) = The direction numbers for element I associated 

with the J-th point and the K-th direction. 


2.1 OVERALL PROGRAM STRUCTURE: DESCRIPTION OF TETRA AND FLOW CHARTS 

The TETRA program consists of a main routine, a function subprogram, 
and twenty-nine subroutines. A brief description of the function subprogram 
and each of the twenty-nine subroutines is given in Figure 5. Figure 6 shows 
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IROUND 

INIT 

SYBSYS 

FLEX 

RBODY 

CONEL 

ELEMl 

ELEM2 

ELEM3 

ELEM4 

ELEM5 

STIFFE 

STIFFt 

INVERT 

MATH 

DEAL 

SINCOS 

FORHIS 

GYROE 

PLOTD 

SCAN 

TILOOP 

ROPROP 

CURRT 

FMODES 

FORCE 

APFOR 

GEN 

MODES 

LISTPP 


Function subprogram which rounds off float inp. point miiuhcrs 

Initializes variables and arrays 

Processes data for the nodal subsystems 

Finds flexible subsystem node shapes 

Computes rigid body mode shapes 

Processes data for the physical connecting elements 

Processes spring-damper (type 1) physical connecting element 
data 

Processes llnk-danper (type 2) physical connecting element data 

Processes rtjb (type 3) physical connecting element data 

Processes engine support-links (type 4) physical connecting 
element data 

Processes data for the uncoupled point spring-damper (type 5) 
physical connecting elements 

Coaq>utes stiffness matrix for engine support element 

Computes stiffness matrix for link elements that are to be 
combined with engine support element 

Matrix Inversion and determinant calculation for engine support' 
links (type 4) physical connecting element 

Matrix multiplication engine support-links Ctype 4) physical 
connecting element 

Processes unbalance load data 

Processes Pcos wt and Psin wt load data 

Processes force-time history load data 

Processes gyroscopic load data 

Processes data for output plot file 

Establishes element/subsystem connections 

Time Integration loop 

(hlculates rotor properties (speed, acceleration, and angular 
displacement) 

Computes current physical displacements, velocities, and modal 
forces 

Provides modal displacements and modal forces 

Computes physical connecting element and gyro element forces 

Computer applied forces 

Conqiutes generalized forces 

Computes modal displacements 

Prints at least a partial listing of the output plot file 


Figure 5. TETRA Subprograms. 




Figure 6. TETRA Program 











the hierarchy of the program, and a flow chart of the overall program is 
given in Figure 7. Flow charts of the main routine, function subprogram, 
and all the subroutines are given in Figures 8 through 38. 

In general, left to right order in Figure 6 is the order in which 
the subroutines are executed. First, subroutine INIT, which initializes 
various program arrays and variables is executed. Then comes the subroutine 
SUBSYS, which processes data for each of the modal subsystems in turn. If a 
particular subsystem is a rigid body subsystem (subsystem 3, 6, or 9), subrou- 
tine RBODY is called by subroutine SUBSYS to calculate the rigid body mode 
shapes. If the subsystem is not a rigid body subsystem (subsystem 1, 2, 4, 5, 

7, 8, 10, or 11), however, subroutine FLEX is called to find the mode shapes. 

Next, subroutine CONEL is called, which processes data for each of the 
physical connecting elements in turn. For each physical connecting element, 
subroutine ELEMl, ELEM2, ELEM3, ELEM4, or ELEM5 will be called from subroutine 
CONEL depending whether the physical connecting element type is 1, 2, 3, A, or 
5, respectively. For type 4 physical connecting elements, subroutines STIFFE, 
STIFFT, INVERT, and MATM are also called from subroutine ELEM4. 

The next subroutines executed are UBAL, SINCOS, FORHIS, GYROE, PLOTD, and 
SCAN (in that order). Subroutine SCAN is the last subroutine called prior to 
when the time integration is performed. The purpose of the subroutines men- 
tioned thus far is, in general, to process data in preparation for the time 
integration. The program progresses very rapidly through subroutine SCAN, 
since the previously mentioned subroutines are executed only once or a small 
number of times. 

Subroutine TILOOP, which consists of the time integration loop, comes 
next. For each time step the subroutines ROPROP, CURRT, FORCE, APFOR, and GEN 
are called by subroutine TILOOP. In addition, subroutine FMODES gets called many 
times by subroutine CURRT and subroutines MODES gets called many times by sub- 
routine GEN. Nearly all of the process time for the program is used by TILOOP 
and the associated subroutines (unless there are only a very few time steps). 

Finally, after the time integration is completed, subroutine LISTPF is 
called, which lists out part or all of the output plot file. 
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Figure 7. TFTRA Flow Chart. 
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Figure 7. 


TETRA Flow Chart (Continued). 
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Figure 7. TETRA Flow Chart (Continued). 
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Figure 7. TETRA Flow Chart (Concluded) 


1 





Read restart input, time integration j^nput , rotor speed and rate 
input, unbalance lead input, P cos uut and P sin a;t load 

input, time-force history input, ^yro load input, and plot file 
input (LIST4 Namelist Section) 


Yes 


Process Rotor Speed and Rate Data 



Read Starting Data 
from Restart File 



Figure 8. Flow Chart of Main Routine. 


18 







Call Subroutine UBAL to Process Unbalance Load Data 


I , 1 


♦ 


Call Subroutine SINCOS to Process P 

cos ijut and P sin Ji^t Load Data 


Call Subroutine FORHIS to Proce 

ss Time-Force History Load Data 



[ 

Call Subroutine GYROE to Process Gyroscopic Load Data 



Call Subroutine PLOTD to Proc 

ess Data for Output Plot File 


Call Subroutine SCAN to Establish Connections Between the 
Physical Connecting and Gyro Elements and the Modal Subsystems 


t 


Call Subroutine TILOOP to 

Perform Time Integration 



Call Subroutine LISTPF to List a 

t Least Part of Output Plot File 



Figure 8. Flow Chart of Main Routine (Concluded). 
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Figure 9. Flow Chart for Function 
IROUND. 









Figure 11. Flow Chart of Subroutine SUBSYS. 









Start 



Figure 12. Flow Chart of Subroutine FLEX. 
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Figure 14 


Flow Chart for Subroutine CONEL 
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Figure 15. Flow Chart of Subroutine ELEMl. 
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Figure 16. Flow Chart of Subroutine ELEM2. 
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Increment Counter for the Number of Type 3 
Physical Connecting Elements NUMEL3 = NUMEL3 + 1 



Find and Print Values 

of Element Variables 



Print Spring Rate, Dead 

Band, and Damping Value 



Define Values for Spring Rate Array SKRATE 
and Damping Value Array CCRATE 


f 

Convert Dead Band fr 
Define Value for D 

om Mils to Inches and 
ead Band Array SRUB 



Figure 17. 


Flow Chart of Subroutine ELEM3. 





Figure 18. Flow Chart of Subroutine ELEM4. 




the^S. 
of LinKSs^ 
(Truss) Load Pachs 
S. ITRUSP - 0? ^ 


Set Values of 
Element Variables 


Define the Global Point Numbers 
for the I, J, and K Ends of the 
Engine Support Element 


Compute the Stiffness 
Matrix AKE4 


Establish the Number of Distinct Points 
on the Other Side of the Link (Truss) 
Elements and Identify their Point 
Numbers 


Was ^ 

Damping Requested 
SsfQELEM 1< 0)? ^ 


Establish if there are J or K Points 
on the Case to Which Link (Truss) 
Elements are not Attached 


Zero out the 
Damping Matrix ADE4 


Compute the Damping 
Matrix ADE4 Based 

on Q-Factor and 
Frequency 


Add the Contrlubtions of the Link 
(Truss) Elements to the Stiffness 
Matrix (Subroutine STIFFT Called 
for Each Link (Truss) Load Path) 


Return to Subrouting 
S. CONEL 


Define Independent-Direct Stiffness 
Matrix AKll 


Define Cross-Stiffness Matrix AK12 


Define Transpose of AK12 (AK21) 


Define Dependent-Direct Stiffness 
Matrix AK22 


Figure 18, Flow Chart of Subroutine ELEM4 (Continued). 
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Figure 18. 


Flow Chart of Subroutine ELEM4 (Concluded). 





Figure 19. Flow Chart for Subroutine ELEM5. 
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Define the X, Y, and Z Coordinates at the 
I and J Ends of the Link (Truss) Element 


Compute the Length of the Link (Truss) Element 


Compute the Direction Cosines for 
the Link (Truss) Element 


Compute the Matrix of Direction Cosine Terms 
ATRUS for the Link (Truss) Element 


Compute the Spring Rate KRATE for the 
Link (Truss) Element 


Define the Stiffness Matrix ATRUS for the Link (Truss) 
Elements (Overwrites Matrix of Direction Cosine Terms) 


Return to Subroutine 
. ELEM4 ^ 


Figure 21. Subroutine STIFFT. 















Figure 24. Flow Chart of Subroutine UBAL 












Figure 26. Flow Chart of Subroutine FORHIS. 





Was Rotor 
Speed Time History 
Requested? ^ 


Print Message that there 
' ► are no Gyro forces 



Is the 

^ Number of 
Gyroscopic Load 
Locations * 0? 



Return to Main Routine 


Figure 27. Flow Chart of Subroutine GYROE, 


41 




Figure 28. Flow Chart of Subroutine PLOTD 





Establish and Print the Connections Between 
the Modal Subsystems and the Physical 
Connecting Elements 



Establish and Print the Connection^ Between the 
Modal Subsystems and the Gyroscopic Elements 


+« 



Figure 29. Flow Chart of Subroutine SCAN. 





Figure 30. Flow Chart for Subroutine TILOOP 





















Figure 32. Flow Chart of Subroutine CURRT 






















Start 



Figure 33. Flow Chart of Subroutine FMODES. 
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Figure 34. Flow Chart of Subroutine FORCE. 
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Calculate and Print the Forcea for the Type 5 Physical 
Connecting Eleneots (Uncoupled Point Sprlng-Danper 
Elements) 



Figure 34. Flow Chart of Subroutine FORCE (Concluded) 
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Zero Out the Pertinent Elements 
of Applied Force Array PFORCE 



Figure 35. Flow Chart of Subroutine APPOR. 
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Figure 36. Flow Chart of Subroutine GEN. 
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Figure 37. Flow Chart of Subroutine MODES 





Figure 38. Flow Chart of Subroutine LISTPF. 
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3.0 MODULAR ELEMENT MODELING 


Modal representations obtained by transforming the degrees*-o£*£reedom 
from real space to modal space are used to define the structural dynamic 
properties of the subsystems. The real space formulations for the subsystems 
are based on beam-like finite element (discrete mass and stiffness) subsystem 
models of the engine rotor (s) and case and a three-dimensional finite element 
model of the pylon, all described with respect to a ground fixed coordinate 
system. Figure 39 shows typical subsystem normal modes and coupling ratios 
for a rotor finite element model. The flexual vibration characteristics of 
equivalent nonrotating shaft(s) are used to define the rotor(s). The whirl 
phenomenon is addressed by using Euler's law for angular motions to establish 
cross-axis coupling forces that are dependent on angular velocities in two 
planes and proportional to the polar inertia; these forces are applied as 
external forces to the modal coordinates that represent the rotor in two 
planes . 

The free-free undamped modes in both the vertical and horizontal planes 
for the rotor (s) and the case are derived from planar finite element models. 

In addition, a single dimension (single degree-of-freedom at each station) 
torsional model is also used to model the case. The rigid body modes for the 
rotor(s) and case in the vertical and horizontal directions can be defined 
either with the "soft spring" rigid body modes or can be defined with separate 
modal subsystems representing computed rigid body modes based on the mass pro- 
perties and the geometry. The torsional direction rigid body modes for the 
case are represented either with the "soft spring" rigid body modes or with a 
rigid body modal subsystem. It will be noted that the case torsional direc- 
tion is important for modeling the engine mounting system load paths. The 
fore-and-aft motions are represented with rigid body modal subsystems for 
both the rotor (s) and the case. The pylon finite element model is used to 
define a set of three-dimensional cantilevered modes. Figure 40 shows the 
global coordinate system, direction numbers, and a typical engine mounting 
arrangement. Table I identifies the physical degree-of-freedom associated 
with each subsystem. Relative to defining the the rigid body mode shapes 
for separate rigid body subsystems, it will be noted that in general an 
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Figure 39. Engine Rotor Subsystem Normal Modes in Two Orthogonal Planes 



Figure 40. Typical Mounting Arrangement and Coordinate System, 









Table I. Physical Global Degree-of -Freedom and Direction Numbers for 
the Subsystems* 



GLOBAL DEGR 

EES-OF-FREEDOM 

SUBSYSTEMS 

Z 

e, 

Y 


X 

®x 

1 

2 

3 

4 

5 

6 


VERTICAL PLANE FLEXIBLE 
ROTOR (S) 

• 

• 





HORIZONTAL PLANE FLEXIBLE 
ROTOR (S) 



• 

• 



RIGID BODY ROTOR (S) 

• 

• 

• 

• 

• 


VERTICAL PLANE FLEXIBLE 
CASE 

• 

• 





HORIZONTAL PLANE FLEXIBLE 
CASE 


i 

• 

• 



RIGID BODY CASE 

• 

• 

• 

• 

# 

• 

TORSIONAL FLEXIBLE CASE 






• 

3D FLEXIBLE PYLON 

• 


• 


• 



Z 0) VERTICAL 
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unconstrained rigid body is free to move in space in six directions. Such 
"rigid body" motion corresponds, in rectilinear coordinates, to translations 
in X, Y, Z and rotations in e^, 6y> and equations describing these 

motions are as follows: 


Applied 

Physical 

Forces 


Inertia Matrix 


Physical Acceler- 
ations of the C.G. 


Fx 


m 

0 

0 

0 

0 

0 



0 

m 

0 

0 

0 

0 

Fz 


0 

0 

m 

0 

0 

0 

Mx 


0 

0 

0 

Ix 

0 

0 



0 

0 

0 

0 


0 

Mz 

L J 


0 

0 

0 

0 

0 

Iz 



These equations are included in the set defined by equation (1). In 
this special case, m£ corresponds to the physical mass properties m, I^, ly, 
^z» ^i corresponds to the physical coordinates X, Y, Z, 6x> 6y> or Oz 
of the C.G., and the k£ terms are equal to zero. Note that for a general 
(nonprincipal) set of axes, the inertia matrix of equation (17) would include 
off-diagonal terms. For a three dimensional body where the planes of symmetry 
pass through the center of mass there is not inertia coupling and the inertia 
matrix corresponds to that of equation (17). The eigenvectors for the six 
rigid body "modes" are shown in Table II. 

In general, the flexible beam like rotor and case finite element models 
will be used to generate modal data for a single plane, say the vertical 
plane. This modal data can also be used to represent the horizontal plane. 
However, it will be noted that the sign of the modal data in either the Y 
or 6z directions must be changed in order to maintain a sign convention that 
is consistent with a right hand global coordinate system. This is shown in 
Figure 41. 
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Table II. Mode Shapes for the Six Rigid Body Modes for a Subsystem. 


PHYSICAL DIRECTION AND 
COORDINATE FOR POINT r 

DIRECTION i 
NUMBER 

COORDINATE 

5 

Xr, 5 

3 

Xr, 3 

1 

Xr, 1 

6 

Xr, 6 

2 

Xr, 2 


Xr. 4 


GENERALIZED COORDINATE q. 

iW 

z 


m 

©z 

1.0 

0 

0 

0 

n 

I0S9 

0 

1.0 

0 

"^■z.r 

0 

H 

0 

0 

1.0 


0 

0 

0 

0 

1.0 

0 

0 

0 

0 

o' 

0 

1.0 

0 

0 

0 

0 

0 

0 

1.0 

01 

02 

03 

04 

05 

06 


COLUMNS CORRESPOND TO RIGID-BODY MODE SHAPES 


POINT r ON RIGID BODY 





























































Figure 41. Modal Planes. 
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Physical Connecting Elements 

The force-displacement and force-velocity relationships for the physical 
connecting elements that connect the modal subsystems are expressed in terms 
of the stiffness and damping matrices shown in equation (18). 

{F} = - [Kel {X} - [Cg] {X} (18) 

Equation (18) defines the physical forces exerted by the physical con- 
necting element on the modal subsystems to which it is connected. 

Generalized Spring-Damper Element 
(Type 1 Physical Connecting Ele^nt) 

This element is associated with two physical points located at arbitrary 
locations in global space. Each of these points is assigned six degrees of 
freedom, three translational displacements (or velocities), and three rota- 
tional displacements (or velocities). Thus, the dimension for both the stiff- 
ness or damping matrix is equal to twelve. 

This element can be used to represent simple uncoupled point springs and 
dampers that are useful in modeling rolling element bearings. Among the other 
load path configurations that this element can be used to model are engine 
frame/sump structures, hydrodynamic bearing definitions with direct and cross 
stiffness and damping, and cross axis stiffness coupling arising from aero- 
dynamic effects. A sketch of this element is shown in Figure 42. 

Typically, the stiffness elements in [K^] are computed via a finite 
element program, or with a closed form solution for an idealized model, or 
are obtained from static or dynamic testing. The damping elements in [C^] 
are computed if analytical expressions are available, or are obtained from 
dynamic testing, or are based on the assumption that the damping is propor- 
tional to the stiffness. In this latter case, [Cg]a[K 0 ] and the proportionality 
can be based on a specified percent of critical damping at a selected fre- 
quency as shown in equation (19). 
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Z(l) 



Figure 42. General Spring-Damper Element (Typical Physical Connecting 
Element ) . 


(\o 



ic^] 


“Of 


tK_] 


(19) 


Op “ specified Q - factor ® 2 c/c 

ID = selected frequency, radians/second. 


Space Link-Damper Element (Type 2 Physical Connecting Element) 

The space link-damper element sketched in Figure 43 is described by 6th 
order stiffness and damping matrices. 


iKe] = f 
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ml 



2 



2 

-nl 

-Im 

-1 

In 

ml 

1 


( 20 ) 


n, m, and 1 are the direction cosines, and the column and row order c6rrespond 
to directions 1, 3, and 5 at points 1 and J, respectively. 

A = cross section area, in^ 

E = Young's modules 

L ® length 

The damping matrix [C^] can be defined either in terms of translational 
(dashpot) damping directed along the axis of the link or by proportional 
damping. 


In the former case, [Cg] ® C [Kg] 
^ Ib/sec 


( 21 ) 


in 


is a specified scalar damping value. 
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z 



X(5) 


Figure 43. 


Space Link-Damper c: 
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In the latter case^ the form of the damping matrix is that shown in 
equation (19). 

Rub Element (Type 3 Physical Connecting Element) 

The large rotor amplitudes associated with blade loss cause heavy rotor- 
case rubs that are usually accompanied by severe local damage to blades and 
case. It is to some degree fortunate that this local damage does occur be- 
cause the loads acting on the bearings and frames are reduced by the action 
of the additional load path between the rotor and case that is provided 
through the mechanism of the heavy rubs. As a consequence of this action, the 
engine may be capable of withstanding without catastrophic structural failure, 
the transient loading induced by the blade loss. 

The rub element allows the mathematical modeling of the nonlinear tip rub 
that includes the dead band displacement interval prior to closure between the 
rotor and case. Upon closure, an equivalent linear spring, representing the 
local case distortion and blade compliance, is used to define the rotor-to- 
case load path. The net result is a bilinear spring with zero slope over the 
dead band and a finite slope over the region of interference. 

Figure 44 shows the relationships between the rub force and the position 
vectors of the rotor and case centers. Neglecting friction, the rub force 
can be defined as a force that acts in the direction defined by the line of 
centers. The vector rub force acting on the rotor can be written as: 

Fr = -T^ (|A| - eo) K - (Vr - Vq)C For |a| > co 

( 21 ) 

Fr “ 0 14| s eo 

Where K = Radial spring rate representing the local case distortion (repre- 
sented by the bulge shown in Figure 44) and the blade compliance. 

eo ® Structural clearance 

A “ (Dj.-Dq) ® vector difference between the position vectors of the 
rotor and case centers. 
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AT TIME = T, VECTORS. Dc AND Dr 
DEFINE POSITIONS OF CASE S ROTOR 
CENTERS. 

VECTOR DIFFERENCE =A= Dr - Dc 
IF lAf >Co THEN RUB 
STRUCTURAL CLEARANCE 

Fr = VECTOR RUB FORCE =- j^(|A| -tV)K (NEGLECTING DAMPING! 
K = SPRING RATE REPRESENTING THE LOCAL CASE DISTORTION AND BLADE 
COMPLIANCE 

Figure 44. Rub Force Model to Represent the Nonlinear Tip Rub Mechanism 
That Includes the Dead Band Interval Prior to Closure. 
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= unit vector in the direction of vector A. 

|d| Vj( and Vq are the vector velocities of the rotor and case centers, 
respectively. 

C is the damping rate. 

In complex notation, & can be written as: 

A ■ Dr - Dc “ (yp + jZr) - (yc + jZc) * (yp -yc) + j(Zr - Zc> 

Where j = unit vector in the z-direction. 


-/ 


2 2 
|A| = amplitude of A = y(yR~yc) (Zr-Zq) 


( 22 ) 


€ 

Fr “ - RfR-Fc) + J(Zr-Zc)][i.O - -^] K - (Vr-Vc)C 

Fr “ ”(yR~yc) [i-o " T?r] ^ ~ [i-o " J?f] 

- c |(yR-yc)j “ jc [(Zr-Zc)J 

Defining A “ 


(23) 


1 . 0 - 


/^yR~yc^^ (Zr-Zc)2 


1 - 


(24) 


as an effective spring coefficient which is dependent on the initial struc- 
tural clearance and the rotor-case relative displacement and case-blade local 
spring rate. A stiffness matrix can be defined as follows: 


r 


r 

Z 
C 

Cc/ 




-A 0 


0 


-1 r 


0 

A 

0 


-A 0 


0 -A 0 


-A 


(25) 


Zi 


2 V ♦ 

Where and F^ are the forces acting on the rotor center in the z and y 
directions , repsectively . 
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Fg and are the forces acting on the case center in the z and y direc- 
tions, respectively. 

Zr, yR, Zq, y^ are the absolute displacements of the rotor and case cen- 
ter. 

The stiffness matrix in equation (25) is also shown in Figure 45 and 
represents the [K^j* matrix in equation (18). The damping forces can be writ- 
ten in matrix form as: 



• • • • 

Zr, yr, Zq, yQ are the absolute velocities of the rotor and case centers. The 
damping matrix in equation (26) represents the [Cg]* matrix in equation (18). 

It will be noted that the rub element defined by equations (25) and (26) 
can also be used to model rotor-to-rotor rubs as well as rotor-to-case rubs. 

In this case, the inside rotor is identified with the subscript r and the out* 
side rotor is identified with subscript C. 


Engine Support Element 


(Type 4 Physical Connecting Element) 


Real aircraft engine mounting systems are quite complex and must be 
modeled with three dimensional models if accurate simulation is to be ob- 
tained. Figure 46 shows an example of an aircraft engine mounting arrangement 
where three mounting planes are utilized. The forward mounting plane takes 
vertical, side, and axial loads. The mid mounting plane takes side and roll 
loads. The aft mounting plane takes vertical loads only. 


^Equations (25) and (26) represent the forces acting on the connecting sub- 
systems and already incorporate the minus signs shown in equation (18). 

Hence, in this case, the [Kg] and [Cg] matrices shown in equations (25) 

and (26) should be multiplied by (-1) to conform to the form of equation (18). 
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Rotor-Case Relative Displacement 
Vector Direction 
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A 0 -A 
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Spring 
Rub Forces 
Acting on 
Rotor and Case 
When /A/ > Gq 



Superscripts 1 and 3 Represent Directions Z and Y 


Figure 45. Rub Element (Type 3 Physical Connecting Element) 



\ 


Vertical 



Figure 46. An Example of an Aircraft Engine Mounting Arrangement. 
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The engine support element is a multipoint » multidirection element that 
provides the capability to model the complex load paths between the engine 
case and the pylon and thus allows for the direct modeling of actual engine 
mount structures. In addition, this element couples the flexible and rigid 
body^centerline modal subsystems, that represent the engine case, to the 
support links that comprise the mounting system. Thus , this element is really 
a combined engine support-case flexibility element that couples the beam-like 
casing model to the 3D pylon subsystem. Figure 47 shows typical engine 
case-link attachment points. The loads at these attachment points induce 
frame/case distortions that increase the effective flexibility that is seen at 
the engine centerline. In this example, the side direction flexibility at the 
forward mounting plane would be increased by the fan frame/case ovalization 
effect and the vertical direction flexibility at the aft mounting plane would 
be increased by the turbine frame/case ovalization effect. 

Figure 48 shows the attachment and centerline points in global space, and 
the load directions for the two mounting planes corresponding to the typical 
mounting arrangement shown in Figure 40. At the forward mounting plane, 
points L, M, and N are the attachment points on the pylon and point I is at 
the engine centerline. The two attachment points on the engine case are 
identified by J and K. In this case, the engine support element is repre- 
sented with the stiffness matrix relating to the points I, L, M, N. The 
displacement terms associated with points J and K are reduced out. Point I 
has six degrees-of-freedom and points L, M, and N each have three degrees- 
of-freedom. It is necessary to eliminate points J and K because the physical 
connecting elements are only used to define the load paths between modal 
subsystems or between modal subsystems and ground. 

The model definition for the portion of the engine support element that 
pertains to the case flexibility (flexibility between the engine centerline 
and the case attach points) is shown in Figure 49. In this representation, 
all of the case flexibility is lumped at the case attach points (points J 
and K) and rigid members are used to transfer the loads to the centerline 
point (point I). The case flexibility is described by three spring rates, 

Ky, K^, and K^. These rates are reciprocals of the case flexibility values 
obtained by either calculation or measurement. In the former case, a finite 
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TURBINE 

MOUNT 

RINC 


IN THIS ILLUSTRATION, THE LOCAL LOADS 
AT THE ATTACHl^ENT POINTS RESULT IN 
ADDITIONAL MOUNTING FLEXIBILITY AT THE 
FORWARD PLANE, IN THE SIDE DIRECTION 
AND AT THE AFT PLANE, IN THE VERTICAL 
DIRECTION. 


Figure 47. Engine Frame/Case Ovallzation Effects Increase the 
Effective Flexibility at the Mounting Planes. 
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^3 



DIRECTION NUMBERS 


FORWARD 


INE 
FRAME/CASE 


THIS SPRING ELEMENT 
COUPLES THE FLEXIBLE 
& RIGID BODY-CENTERLINE 
SUBSYSTEMS THAT REPRESENT 
THE ENGINE CASE TO THE 
CASE ATTACH POINTS. 

THE LOCATIONS OF THE 
CASE AH ACM POINTS 
ARE ARBITRARY AT THE 
MOUNTING STATION. 


(FWD, LOOKING AFT) 



POINT I. 

Figure 49, Case Flexibility Represented by Multipoint/Direction Physical 
Spring Element (Which Defines a Portion of the Engine Support 
Element) « 
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element model of the frame is typically used to define the flexibility values 
For example 9 would be defined by applying an n * 1 shear flow restraining 
the frame model at the case attach point locations* The vertical displace- 
ment of a horizontal diameter would then be used to establish the spring rate 
Kv, Figure 50 shows the twelfth order stiffness matrix that describes the 
case flexibility* The terms» a» b» c, d, e, and f are multipliers used to 
proportion the spring rates, and the association between the spring rates 
Kh, Kfl, and the lumped spring location points J and K are shown below. 

Case Distortion Spring Rates 
Kv ® Vertical Spring Rate 
Kh - Horizontal Spring Rate 
Ka = Axial Spring Rate 


In conjunction with this scheme, the following restraints must be fol- 
lowed . 


Multipliers 


Spring Rates 
(Ib/in. ) 

Point 

J 

K 

Kv 

a 

c 

Kh 

b 

d 

Ka 

e 

f 


a c *= 1.0 

b d * 1.0 

e + f = 1.0 


For symmetrical mounting, a=b=c®d®e*f®.5. 


The engine support element is a variable geometry element that is formed 
by combining the spring element stiffness matrix (Kf J of Figure 50 with the 
stiffness matrices for connecting link elements. The stiffness matrices for 
the individual link elments have the same form as the stiffness matrix shown 
in equation (20). 


[KcJ 

iKel 


[Kf] + Kf] 

K„ - K,„ k“^ K,,1 

L 11 12 22 2lJ 


hi 1 ‘^12 


(27) 

(28) 


[Kf] 

[Kl] 


spring element stiffness matrix representing the case 
stiffness matrix for the links 
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in to the point numbers and the direction numbers, respectively. 


Matrix for the Physical Spring Element Shown 







































































































(KcJ * totAl unreduced stiffness matrix for the assembly 

[Ke] * reduced stiffness matrix for the engine support element 

The matrix [Kc] in equation (27) represents the total unreduced stiff- 
ness matrix which includes the case points J and R. Equation (28) represents 
the matrix reduction that is performed to eliminate points J and K. Figure 48 
shows two examples of engine support element configurations. It will be noted 
that there are no local moment load paths at points J and R or at the link 
attachments at the pylon. 

The damping matrix [Ce] for the engine support element is defined in 
terms of proporational damping and has the form of the damping matrix shown 
in equation (19). 


Uncoupled Point Spring-Damper Element (Type 5 
^ysical Connecting Element) 

This element is for connection of two points with an uncoupled spring 
and damper. Each point has five degrees of freedom, three translational 
displacements (or velocities) and two rotational displacements (or velocities) 
(see Figure 51). These elements are typically used to connect the centerlines 
of beam - like modal subsystems. Because of the lack of load path coupling, 
good modeling practice infers that the points being connected by this element 
should be coincident in space. 


The equations for the forces at the I and J end points of this element 
become rather simple due to the lack of coupling: 
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I 

Jt 
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e\y y/ e\y y/ 
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I and J Points 
are Coincident 


z (Vertical Direction) 



Figure 51. Uncoupled Point Spring - Damper Element (Type 5 Physical 
Connecting Element). 




where 

^ Force that the element exerts on the modal subsystem or ground. 

subscript i is for direction (x, y, z, 0^^, or By) and superscript 
j (I or J) indicates the I-end point or the J-end point. 

K£ = Spring constant. Subcript i is for direction (x, y, z, 0y, or Bg) . 

C£ = Damping coefficient. Subscript i is for direction (x, y, z, 0^ , 
and 0^ ) . 

• • • • ^ 

X J , 6z ® Displacements in the x, y, z, 6y, and directions, 

respectively. Superscript i (I or J) indicates the I-end point 
or the J-end point . 

•1 i i *i *i 

X , y , z , 6y , ©2 ® Velocities in the x, y, z, 6y, and direc- 

tions respectively. Superscript i (I or J) indicates the I-end 
point or the J-end point . 

As for the generalized spring-damper element (type 1 physical connecting 
element), the spring constants (K^) are computed via a finite element program, 
or with a closed form solution for an idealized model, or are obtained from 
static or dynamic testing. The damping coefficients (C^) are computed if 
analytical expressions are available, or are obtained from dynamic testing, or 
are based on the assumption that the damping is proportional to the stiffness 
(see equation 19). 
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Note that the Type 5 element is a special case of the Type 1 element 
(generalized spring-damper element). The Type 1 element could always be 
used instead of the Type 3 element. However, it is recommended that the 
Type 5 element be used rather than Type 1 wherever possible. The equations 
for the forces for the Type 5 element are short and simple, but the Type 1 
element requires matrix multiplication (see equation 18) (the Type 1 element 
stiffness and damping matrices are each 12 x 12) , which makes the calculation 
of forces more lengthy. Thus, the program runs more efficiently (and saves 
on cost) using Type 5 elements instead of Type 1 where possible. The input 
is also simpler for the Type 5 element than for the Type 1 element. 


Gyro Element 


The gyro element models the cross-axis forces due to Coliolis accelera- 
tion and addresses to general whirl motion defined by the response in two 
planes. These cross-axis forces are mathematically treated as "Right Hand 
Side Forces," or externally applied forces, and couple the vertical and 
horizontal plane rotor subsystem models • Figure 52 shows the damping 
matrix used to define the gyro forces acting on the connecting subsystems. 
This matrix represents -1 times the [Ce] matrix in equation (18) and is 
derived from Euler's equations for rotational motions. For a spinning disk 
on a whirling shaft, Euler’s equations of motion in a fixed frame can be 
written as: 



My and Mz are the applied moments. 


(29) 



the inertia matrix that defines the uncoupled rotary 
interia moments acting in the two planes. 

the damping matrix that defines the coupled-velocity 
dependent gyroscopic moments. 


The affects of the uncoupled rotary inertia are implicitly included in 
the rotor subsystem modal data and the velocity dependent moments are treated 
as the applied physical forces shown in Figure 52. 
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Figure 52. 


Gyro Element - Gyro Forces are Computed with a Velocity-Dependent Element 


4.0 INPUT AND OUTPUT FILES FOR TETRA 


The main input file of the TETRA program consists of namelist input. 

There are four namelist names * LISTl, LIST2, LIST3, and LIST4 (must be ar- 
ranged in that order). First must come the LISTl namelist input (for ident- 
ifying input and input for points not located on the modal subsystems). There 
must be one and only one LISTl section. Next come the LIST2 namelist sections 
(for the modal subsystems). One LIST2 section is required for each modal sub- 
system (there must be at least one modal subsystem and at most 11). The LIST3 
namelist sections (if any) (for the physical connecting elements) come next. 

There is one LIST3 section for each physical connecting element. There may be 
zero physical connecting elements. The current upper limits are a maximum of 
five Type 1 physical connecting elements, ten of Type 2, ten of Type 3, four 
of Type 4, and ten of Type 5. Finally, there must be one and only one LIST4 sec- 
tion (which contains the restart input, the time integration Input, rotor 
speed and rate input, unbalance load input, P cos cut and P sin cut load input, 
time-force history load input, gyro load input, and plot file input). 

The namelist input for the main input file must always start in Column 2. 
Within a given namelist section, the input variables can be in any order. The 
file code for the main input file is the standard 03 file. A listing of the 
main input file is provided at the beginning of the printed output of TETRA for 
user convenience. 

An interface program has been written to read the output file from the 
NASTRAN program, and then to automatically generate all or most of the modal 
subsystem (LIST2) input (see Section 6.0). The modal subsystem input thus 
generated by the interface program must be merged with the main input file 
for TETRA in the correct order as explained above. 

If the run is a restart run, one additional input file is required, 
namely, the restart file generated by the original run. This input restart 
file is read in using file code 22, so the user must assign the input restart 
file to File 22 (@ ASG control card). 
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Two output files nay be generated by a TETRA run. A plot file (File 
Code 23) is generated provided the user selects the IPLOT*! option (see 
Type P-1 namelist inupt sheet). A listing of all or part of the output plot 
file contents is provided at the end of each run as an aid for the develop- 
ment of software to plot the TETRA results. (If 21 or less times are written 
to the output plot file the entire file is listed, but if more than 21 times 
are written to the output plot file, only the first 21 are listed due to space 
limitations.) The following pages define the variables which are on the out- 
put plot file. An output restart file (File Code 24) is also generated. If 
the user wishes to save these output files he must assign them (@ ASG control 
card). Both output files are unformatted output. 

The output restart file is always a very short file and would best be 
saved on disk. The output plot file, however, may be a lengthy file and might 
be saved on tape or disk. 


Plot File Format for TETRA Program 

NP,NRE,NEL,NUMT 
IPPLOT^,IDPLOTi ,XPTi ,YPTi ,ZPTj^ 

IPPLOTnp , IDPLOTnp , XPTnp , ZPTnp 
ILEM 3i , ILEM3i , ILEM3i 

ILEM3nre , ILEMSjjre , ILEM3nre 
ILEM i , IPTi , IDIRi 

ILEMnel . IPTneL » IDIRnel 

TIME ( 1 ) , SPEEDK 1 )STEEDD ( 1 ) , THETAI ( 1 ) , THETAD ( 1 ) 
Xi(1),VELi(1),FM0Di(1) 

t 

Xnp ( 1 ) , VELnp ( 1 ) , FMODnp ( I ) 

DMAGi ( I ) , CLEARi ( 1 ) , FMAGi ( 1 ) 


DMAGnre ( 1 ) . CLE ARnre ( 1 ) , FMAGrre ( 1 ) 

FELEMi ( 1 ) , FELEMi ( 1 ) , FELEMi ( 1 ) 

« 


# 
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TIME(NUMTE>,SPEEDI(NIJMT) ,SPEEDD(NUMT) ,THETAI(NUMT) ,THETAD(NIMT) 
Xi(NUMT),VEL, (NUMT) ,FMOD, (NUMT) 


Xiip(NUMT) ,VELnp(NUMT) ,FMODnp(NUMT) 

DMAGi ( NUMT ) , CLEARi ( NUMT ) , FMAGi ( NUMT ) 

4 

m 

• 

DMAGnre ( NUMT ) , CLE ARnre ( NUMT ) , FMAG^Re ( NUMT ) 
FELEMl(NUMT) , FELEMK NUMT ) ,FELEMi( NUMT) 

* 

FELEMnel ( NUMT ) , FELEMnel ( NUMT ) , FELEMnel ( NUMT ) 


Definitions 


NP = Number of (global point, global direction) pairs for which the displace- 
ment, velocity, and modal force is written on the plot file 

NRE = Number of Type 3 physical connecting elements (rub elements) for which 
the relative displacement magnitude, clearance, and force magnitude are 
written on the plot file 

NEL = Number of (element number, global point, global direction) trios for 
which the physical connecting element or gyro element force is written 
on the plot file 

NUMT = Number of timne steps on the plot file 

IPPLOTi = Global point number for the i-th (global point, global direction) 
pair on the plot file 

IDPLOT^ = Global direction number for the i-th (global point, global direc- 
tion) pair on the plot file 

XPT£ ,YPTi ,ZPTi = X,Y and 2 coordinates (global system) respectively for the 
i-th (global point, global direction) pair on the plot file 

ILEM3i = Element number for the i-th Type 3 physical connecting element (rub 

element) for which the displacement magnitude, clearance, and force magni- 
tude are written on the plot file 

ILEMi = Element number for the i-th physical connecting element or gyro ele- 
ment for which the element force is written on the plot file 

IPTi = Global point number for the i-th physical connecting element or gyro 
element for which the element force is written on the plot file 
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IDIR^ * Global direction number for the i-th physical connecting element or 
gyro element for which the element force is written on the plot file 

TIME(I) “ Time (seconds) for the I-TH time step on the* plot file 

SPEEDl(l) » Independent rotor speed (rpm) for the I-TH time step on the plot 
file 

SPEEDD(l) B Dependent rotor speed (rpm) for the 1-TH time step on the plot 
file 

THETAI(I) = Independent rotor angular displacement (revolutions) for the I-TH 
time step on the plot file 

THETAD(I) = Dependent rotor angular displacement (revolutions) for the I-TH 
time step on the plot file 

X^d) ’= Displacement (inches or radians) for the i-th (global point, global 
direction) pair and for the i-th time step on the plot file 

VELi(l) = Velocity (in/sec or rad/sec) for the i-th (global point, global 
direction) pair and for the I-TH time step on the plot file 

FMODi(I) = Modal force (lb or in-lb) for the i-th (global point, global 
direction) pair and for the I-TH time step on the plot file 

DMAG^d) ° Relative displacement mangitude (inches) for the i-th Type 3 phys- 
ical connecting element (rub element) and for the I-TH time step on the 
plot file 

CLEARi(l) = Clearance (inches) for the i-th Type 3 physical connecting element 
(rub element) and for the I-TH time step on the plot file 

FMAGid) = Force mangitude (pounds) for the i-th Type 3 physical connecting 
element (rub element) and for the I-TH time step on the plot file 

FELEMid) = Force (lb or in-lb) for the i-th (element number, global point 

number, global direction number) trio and for the I-th time step on the 
plot file 
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5.0 INPUT SHEETS 


Following is a discussion of the namelist input sheets for the main input 
file of TETRA. Input data should follow the order given for the input sheets 
(starting with Type A input, then Type B, etc.). Always begin the namelist 
input in Column 2 of each line. 

Identification Sheet Type A 

Each TETRA run, no matter %ihat the engine system or the type of analysis, 
must have one of these sheets. The four lines are filled in with the informa- 
tion indicated. The last two lines give the user the opportunity to include 
a descriptive name for the case. All of these lines are reproduced on the 
output . 
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Page of 

NAMELIST 
Type A 

TETRA 

IDENTIFYING INPUT 


V 

$LIST1 

NAME= ' @DATE© ^©CHARGE© 

ADDRES= ' ^©MAIL DROP© ^©EXT. 

IDENTIC ■ 

IDENT2= ' 


Maximum of 60 characters enclosed within apostrophes for each of the 
above variables. 



TYPE B INPUT 


Physical Points Not Located on Modal Subsystems 

This input sheet defines the boundary conditions of the problem. The 
points which the user desires to ground are defined here. If the code is 
specified as fixed (CODE^l), then the point is fixed in 6 DOF. If the code 
is specified as free (CODE®0), then the point is free to move in 6 DOF. The points 
defined on this sheet may not be located on the modal subsystems but can only 
be located at ground or at the junctions of links and engine support elements. 



Page of 
NAMELTTT 
Type B 

TETRA 


PHYSICAL POINTS NOT LOCATED ON MODAL SUBSYSTEMS 


Z (vertical direction) 



X (axial direction, positive forward) 


V 

PP(1.1)= 


POINT NUMBER 
NOT ON MODAL 
SUBSYSTEM * 


CODE 

0=FREE 

1=FIXED 


COORDINATES RELATIVE 
TO GLOBAL SYSTEM 
X Y 2 



* These points include ground points and points between link and engine 
support elements. It will be noted that points on modal subsystems 
cannot be attached directly to ground, but can be attached to ground 
through physical connecting elements. 




Modal Subsystem Input (Type Cl Through C-11) 


Modal subsystem input is accomplished using the LIST2 namelist section. 
One LIST 2 namelist section is required for each modal subsystem. The number 
of modal subsystems present is counted by the program automatically. 

The user may input a title (this is optional) for each subsystem using 
the TITLE= variable. Up to 60 characters may be enclosed within the 
apostrophes. 

The subsystem number (ISUB variable) is required. The subsystem number 
must be an integer number between 1 and 11. The subsystem numbers are not 
arbitrary - each subsystem number represents a different type of modal sub- 
system with differing degrees of freedom (see Figure 53). The user chooses 
which of the modal subsystems to use, and can arrange these modal subsystems 
in any order. The user must have input for at least one modal subsystem. 

Note that Subsystems 1-3 represent "Rotor 1" and Subsystems 4-6 repre- 
sent "Rotor 2" (see Figure 53). Of course, the user's model may have only one 
rotor, in which case either Rotor 1 (Subsystems 1-3) or Rotor 2 (Subsys- 
tems 4-6) could be used. If the engine is a dual spool engine, however, the 
user might want to use both "Rotor 1" and "Rotor 2" in his model. In this 
instance, "Rotor 1" might represent, say, the low pressure system (fan, low 
pressure compressor, shaft, and low pressure turbine), and Rotor 2 might rep- 
resent the high pressure system (high pressure compressor, shaft, and high 
pressure turbine). Or vice-versa - the choice is up to the user. The user 
has a great deal of flexibility and can model with these subsystems virtually 
any jet engine configuration . 

Each subsystem requires input for points on the subsystem (PTS array) 

(see input sheets C-2, C-5, C-8 and C-11). Each subsystem must have at least 
one point. The current upper limits are 10 for the rotor subsystems (Num- 
bers 1-6), 20 for the case (or housing) subsystems (Numbers 7-10), and 10 for 
the pylon subsystem (Number 11). The rotor and case subsystems (Subsystems 
I-IO) represent engine centerline models. Thus, if the user chooses his 
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SUBSYSTEM 

NUMBER 

DESCRIPTION 

NUMBER OF 
DEGREES OF 
FREEDOM 

DIRECTIONS 

1 

Rotor 1 Vertical Plane 

2 

z, 0 

Y 


2 

Rotor 1 Horizontal Plane 

2 



3 

Rotor 1 Rigid Body 

5 

X, Y, e^, Z, 


h 

Rotor 2 Vertical Plane 

2 

z, K 


5 

Rotor 2 Horizontal Plane 

2 

CD 


6 

Rotor 2 Rigid Body 

5 

X, Y, 9^, Z, 


7 

Case Vertical Plane 

2 

>< 

CD 


8 

Case Horizontal Plane 

2 



9 

Case Rigid Body 

6 

X 

CD 

X 

CD 1 

®z 

10 

Case Torsional 

1 

ex 


11 

Pylon 

3 

X, Y, Z 



Figure 53. Modal Subsystem Summary. 
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X**axl8 (axial direction) to coincide with the engine centerline, then the y 
and z coordinates for the points on Subsystems 1-10 will all be zero. This 
is not true for the pylon subystem (Number 11), however. If the subsystem 
referecne point does not coincide with the global origin, the user should in- 
put the coordinates (XREF, YREF and ZREF variables) of the subsystem refer- 
ence point, so the program can make adjustment. 
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PAGE OF 

NAMELIST 
TYPE C-1 


MODAL SUBSYSTEM INPUT FOR VERTICAL AND 
HORIZONTAL PLANE SUBSYSTEMS 

THIS INPUT SHEET APPLIES FOR THE FOLLOWING MODAL SUBSYSTEMS: 

NL'MBER 1 (ROTOR 1 VERTICAL PL/'^E) 

NUMBER 2 (ROTOR 1 HORIZONTAL PLANE) 

NUMBER 4 (ROTOR 2 VERTICAL PLANE) 

NUMBER 5 (ROTOR 2 HORIZONTAL PLANE) 

NUMBER 7 (CASE VERTICAL PLANE) 

NUMBER 8 (CASE HORIZONTAL PLANE) 


Z 



Z (VERTICAL DIRECTION) 



(HORIZONTAL DIRECTION 
POSITIVE FORWARD) 


NOTE - THE PHYSICAL MODEL (MODAL SOURCE) FOR THE SUBSYSTEM IS ALWAYS A 
VERTICAL PLANE MODEL. 



SLIST2 

TITLE= 

ISUB= 


I 


SUBSYSTEM NUMBER (1,2, 4, 5, 7, or 8) 


> 


ENTER THE COORDINATES (INCHES) OF THE SUBSYSTEM REFERENCE POINT; 

XREF= ,X COORDINATE (0 ASSUMED) 

YREF» ,Y COORDINATE (0 ASSUMED) 

ZREF= ,Z COORDINATE (0 ASSUIED) 


Maxltninn of 60 characters enclosed within apostrophes for the title. 
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PAGE OF 

NAMELIST 
TYPE C-2 


MODAL SUBSYSTEM INPUT FOR VERTICAL AND 
HORIZONTAL PLANE SUBSYSTEMS (CONTINUED) 


ENTER THE COORDINATES OF THE TETH.A POINTS. >UXIFIUM OF 10 
POINTS FOR SUBSYSTEMS 1,2, A, AND 5. MAXIMUM OF 20 POINTS FOR 
SUBSYSTEMS 7 AND 8. Y and Z COORDINATES NORMALLY 0. 
n=NUMBER OF POINTS IN THE SUBSYSTEM. 



ENTER THE VALUES BELOW FOR EACH MODE. MAXIMLT1 OF 
30 MODES. 

N=NUMBER OF SUBSYSTEM MODES. 
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NAMELIST 
TYPE C-3 


MODAL SUBSYSTEM INPUT FOR VERTICAL AND 
HORIZONTAL PLANE SUBSYSTEMS (CONTINUED) 


ENTER THE REQUIRED MODE SHAPES BELOW. 
n»NUMBER OF POINTS IN THE SUBSYSTEM. 
N=NL’MBER OF SUBSYSTEM MODES. 


TRANSLATION T 

SLOPE 8 

SHEAR V 

MOMENT M 

INCHES 

RADIANS 

POLTJDS 

IN-LB 




PS 

1 

u 

2 

O 

» Z 

J 

Q 

• 

W XT* 
2 2 

• 



w M 

• 

S 

n 


VH(1,1,2)= 

> 

» 

9 

9 

» 

9 

9 

9 


9 

9 

9 

9 


* 

<1 

9 

9 


PS 

1 


2 

CSI P 

3 

g* 






S £ 

r ^ 


O J 




o 

-3 

n 


VH(1,1,N)= 


cc 

1 

- - .9 

9 

9 

9 

* 1 

2 

_ 9 

9 

9 

9 

3 

9 

9 

9 

9 

Q 

• 

• 

• 





• 

• 





• 

• 

n 

• 

• 

« 

• 

9 

9 

9 




PAGE OF 

NAMELIST 
lYPE C-4 


MODAL SUBSYSTEM INPUT FOR RIGID BODY SUBSYSTEMS 


THIS INPUT SHEET APPLIES FOR MODAL SUBSYSTEM 3 (ROTOR 1 RIGID BODY), MOD.AL 
SUBSYSTEM 6 (ROTOR 2 RIGID BODY), AND MODAL SUBSYSTEM 9 (CASE RIGID BODY). 


Z 



Z (VERTICAL DIRECTION) 



X (AXIAL 
DIRECTION) 


NOTE: MOTION IN THE 3X DIRECTION IS NOT CONSIDERED FOR TOE ROTOR SUBSYSTEMS 

(3 AND 6) . MOTION IN THIS DIRECTION CAN BE CONSIDERED FOR THE CASE SUBSYSTEM 
(9), HOWEVER. 


’W 


$ END 
$LIST2 

TITLE= ’ 

ISUB= , SUBSYSTEM NUMBER (3,6, OR 9) 

ENTER THE COORDINATES (INCWES) OF THE SUBSYSTEM REFERENCE POINT: 


XREF= 

YREF= 

ZREF= 


j X Coordinate (0 assumed) 
^ Y Coordinate (0 assumed) 
J Z Coordinate (0 assumed) 


$ 


Maximum of 60 characters enclosed within apostrophes for the title. 


96 



LOCAL POINT NUMBER 


PAGE OF 

NAMELIST 
TYPE C-5 

MODAL SUBSYSTEM INPUT FOR RIGID BODY SUBSYSTE?.1S (CONTINUED) 


ENTER COORDINATES OF TETRA POINTS. MAXIMUM OF 10 POINTS FOR 
SUBSYSTEMS 3 AND 6 AND MAXIMUM OF 20 POINTS FOR SUBSYSTEM 9. 
Y AND Z COORDINATES NORMALLY 0. 
n= NUMBER OF SUBSYSTEM POINTS. 


TETRA 

COORDINATES (INCHES) WITH RESPECT 

POINT 

TO SUBSYSTEM REFERENCE POINT 

NmSER 

X I Y Z 


n 


PTS(1,1) = 


ENTER THE SUBSYSTEM CENTER OF GRAVITY COORDINATES (INCHES) WITH 
RESPECT TO SUBSYSTEM REFERENCE POINT: 

X0G= , X COORDINATE (0 ASSU:^IED) 

Y0G= , Y COORDINATE (0 ASSUMED) 

ZCG= , Z COORDINATE (0 ASSlT.tED) 
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PAGE OF 

NASIELIST 
TyPE C-6 


MODAL SUBSYSTEM INPUT FOR RIGID BODY SUBSYSTE?.1S (OC»rriNUED) 


ENTER THE SUBSYSTEM PHYSICAL WEIGHT PROPERTIES FOR EACH OF THE SIX 
DIRECTIONS. IF MOTION IN A PARTICULAR DIRECTION IS TO BE NEGLECTED, 
ENTER ZERO FOR THE PHYSICAL WEIGHT PROPERTY IN THAT DIRECTION. 

(SINCE THE 9X DIRECTION IS ALWAYS NEGLECTED FOR SUBSYSTEMS 3 AND 6, 

WMIX SHOULD ALWAYS BE SET TO ZERO OR OMITTED FOR SUBSYSTEMS 3 AND 6), 
LEAST ONE OF THE PHYSICAL WEIGHTS MUST BE NON-ZERO. 




^ 

B 

)~i 

X 

ux= 

Y 

\VY= 

Z 

\VZ = 

•-HOW 
CO M ft, 

0 

U 

0X 

TO1IX = 

> W O 
S a: 

a 

M 

ev 

WMIY = 

CU O4 

Q 

6z 

WMIZ = 


WEIGHT, POUNDS (0 ASSUMED) 

WEIGHT, POUNDS (0 ASSUJIED) 

WEIGHT, POUNDS (0 ASSUMED) 

POLAR MOMENT OF INERTIA, LB-IN^ (0 ASSUMED) 

MOMENT OF INERTIA ABOUT Y AXIS, LB-IN^ (0 ASSUMED) 
MOMENT OF INERTIA ABOUT Z AXIS, LB-IN^ (0 ASSUMED) 
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PAGE OF 

NAMELIST 
TYPE C-7 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 10 
(CASE TORSIONAL SUBSYSTEM) 


Z (VERTICAL DIRECTION) 



V 

$END 

$LIST2 

TITLE= * 

ISUB=10, 

ENTER THE COORDINATES (INCHES) OF THE SUBSYSTDI REFERENCE POINT: 

XREF= ,X COORDINATE (0 ASSUMED) 

YREF= ,Y COORDINATE (0 ASSUTIED) 

ZREF= ,Z COORDINATE (0 ASSUMED) 


Maximum of 60 characters enclosed within apostrophes for the title. 
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PAGE 

NAMELIST 
TYPE C-8 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 10 
(CASE TORSIONAL SUBSYSTEM) (CONTINUED) 


ENTER THE COORDINATES OF THE TETRA POINTS. MAXIMLTI OF 20 POINTS 
X AND Y COORDINATES NORtLALLY 0 
n=NUMBER OF POINTS IN THE SUBSYSTEM 


TETRA 

POINT 

NL»MBER 

COORDINATES (INCHES) WITH RESPECT 
TO SUBSYSTEM REFERZN'CE POINT 

X 

Y 

Z 


V 


PTS(1,1)= 


; 

1 

f 

* 

9 

9 

ca 

M 

2 

9 

9 

9 

9 


3 

9 

9 

9 

9 

H 

2 

5 

S 

n 

9 

m 

> 

9 

9 

9 


ENTER THE VALUES BELOW FOR EACH MODE . MAXIMUM OF 30 
N=NUMBER OF SUBSYSTEM MODES. 


oi 

1 



B 

2 


3 



» 


a 


o 


s 


•ji 


< 


O 


o 

9^ 

N 


FREQUENCY 

(RPM) 

POTENTIAL 

ENERGY 

Q-FACTOR 

MODE TYPE 
0=FLEXIBLE 
1=RIGID BODY 

V 





XM0DES(1,1) 

9 

9 



9 

9 

» 9 



9 

] 

» 9 


f 

9 

9 

9 



9 


MODES . 





PAGE 

NAMELIST 
TYPE C-9 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 10 
(CASE TORSIONAL SUBSYSTEM) (CONTINUED) 


ENTER THE REQUIRED MODE SHAPES BELOW. 
n=NUMBER OF POINTS IN THE SUBSYSTEM. 
N=NUMBER OF SUBSYSTEM MODES. 


Tl«ST 6 

MOMENT 

RADIANS 

IN-LB 


^R(l,l,l) = 



a 

1 

> 

9 


S 

2 

9 

9 


D 

:z 

3 

9 

9 

o 

g 

< 

M 

2 




CJ 

q 






5 

o 

n 

9 

9 


Z- 

« 

1 

1 T0R(1,1.2)= 

9 

9 



2 

9 

9 

1 

:d 

3 

9 

9 






Q 

M 




i 

2 









s 

o 


n 

9 

9 







n 

ft: 

1 

c: 

u . 




2 

i 

p 

3 


H 

• 




9 

Hi 


i 

2 





5 

S 


o 

o 




n 


T0R(1,1,N)= 

9 

■ ■ 1 

> 


> 

> 

♦ 









PAGE OF 

NAMELIST 
T5fPE C-10 


MODAL SUBSYSTEM INPOT FOR SUBSYSTEM 11 
PYLON SUBSYSTEM 



X (AXIAL DIRECTION, POSITIVE FORUARD) 


SEND 

$LIST2 

TITLE= ’ 

ISUB=11, 

ENTER THE COORDINATES (INCHES) OF THE SUBSYSTEM REFERENCE POINT; 
XREF= , X COORDINATE (0 ASSUMED) 

YREF= , Y COORDINATE (0 ASSUMED) 

ZREF= , Z COORDINATE (0 ASSUMED) 


Maximum of 60 characters enclosed within apostrophes for the title. 
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NAMELIST 
TYPE C-11 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 11 
(PYLON SUBSYSTEM) (CONTINUED) 


ENTER THE CXX)RDINATES OF THE TETRA POINTS. MAXIMUM OF 10 POINTS. 
tuzfnJMBER OF POINTS IN THE SUBSYSTEM 


TETRA 

POINT 

NUMBER 

COORDINATES (INCHES) WITH RESPECT 
TO SUBSYSTEM REFERENCE POINT 

X 

Y 

Z 




PTS(1,1) = 




_ _ 9 

9 

9 

9 


1 






2 

9 

9 

9 _ . 

9 

u 

3 

9 

9 

9 

- 1 

H 


« 

• 




s 

g 

• 

• 

• 





• 





a 

• 

• 




o 

# 

" i 

• 

f 

.9 


..9 


ENTER THE VALUES BELOV.' FOR EACH MODE. MAXIMUM OF 30 MODES. 
N=NUMBER OF SUBSYSTEM MODES. 



1 

1 

2 

i 

3 

u 

# 

Q 


O 

0 

s 



0 

s 

0 

s 

N 


FREQUENCY 

(RPM) 

POTENTIAL 

ENERGY 

Q-FACTOR 

MODE TYPE 
0=FLEXIBLE 
I^IGID BODY 

XMODES (1,1)= 

» » 

9 

f 

9 

_ J 

J 


9 

9 

_ . 9 

t 






LOCAL MODE N LOCAL MODE 2 LOCAL MODE 1 

LOCAL POINT NUMBER LOCAL POINT NUMBER LOCAL POINT NUMBER 
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NAMELIST 
TYPE C-12 


1 

2 

3 


1 

2 

3 


n 


1 

2 

3 


MODAL SUBSYSTEM INPUT FOR SUBSYSTEM 11 


PYLON SUBSlfSTEM (CONTINUED) 

ENTER THE REQUIRED MODE SHAPES BELOW. 
n= NUMBER OF POINTS IN THE SUBSYSTEM, 
N= NUMBER OF SUBSYSTEM MODES. 

TRANSLATION 

TRANSLATION 

TRANSLATION 


Ty 


INCHES 

INCHES 

INCHES 

PYL(1,1,1) = 

.3 


9 9 

S 


9 9 

. 3 


9 . 9 


* * ^ 


PYL(1,1,2) = 



_ _ 9 

_ _» 


9 

9 

J . _ 

_ 9 

f 


-• 

PYL(1,1,N) = 

. .3 

9 . _ 

9 

.9 

9 

9 . 

3 

. 9 

... 9 

J _ 


.3 

9 
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Now for a discussion of the modal subsystem input that is particular to 
certain subsystems* First) consider the vertical and horizontal plane subsys-* 
terns (Subsystems 1, 2, 4, 5, 7, and 8). Input sheets C-1, C-2, and C-3 apply 
for these subsystems. The XHODES array (see input sheet C-2) is required to 
supply data for each mode. These subsystems must have data for at least one 
and a maximum of 30 modes. The frequency and potential energy for each mode 
are entered in the XMODES array. These values are used to compute the general"* 
ized weight for the mode. Next, Q-f actor Qf is input for the mode. The Q- 
factor determines the modal damping for the mode (provided the mode type code 
equals 0) and can be defined as follows: 



where C^damping coefficient and CQ^critical damping coefficient. If the 
user wishes to neglect damping) he should input Qf«0) which is a signal to 
the computer to neglect damping. Next, the user inputs a code specifying mode 
type (O^flexible, or l=rigid body). The code 1 is used for "soft spring" 
rigid body modes whch are used to approximate the true rigid body modes. If 
this code equals 1 (rigid body)) then the modal stiffness and modal damping 
are both set to zero (regardless of what was inputted for Q-factor). If) how- 
ever) the mode type code equals 0 (flexible)) the modal stiffness is calcu- 
lated based on the potential energy) and the modal damping is calculated based 
on the Q-f actor, generalized stiffness, and frequency. 


Also the mode shape input (VH array - see page C-3) is required for the 
vertical and horizontal plane subsystems* 

If the user is using the NASTRAN/TETRA interface program to generate in- 
put data for the vertical and horizontal plane subsystems, he should make his 
NASTRAN model a vertical plane model as shorn on input sheet C-1. Since the 
rotors are always rotationally symmetric (and often the case is also), the 


105 



same vertical plane NASTRAN model can usually be used to generate modal sub-* 
system input for both the vertical and horizontal plane subsystems. In this 
instance, the only difference in the TETRA modal subsystem input for the vert- 
ical and horizontal plane subsystems of a given component will be the title 
and the subsystem number. 

Next, consider the rigid body subsystems (Subsystems 3, 6 and 9). Input 
sheets C-4, C-5y and C -6 apply for these subsystems. The center of gravity 
coordinates with respect to the subsystem reference point (XCG, YCG and ZCG 
variables) are required for these subsystems. In addition, the physical 
weight properties should be input (see input sheet C- 6 ), but only for those 
directions for which rigid body motion is to be considered. Generalized co- 
ordinates are assigned for each direction for which the physical weight prop- 
erty is non-zero, but are not assigned if the physical weight property was set 
equal to zero or omitted. Since the Ox direction is always neglected for 
Subsystems 3 and 6 , WMIX should be set to zero or omitted for Subsystems 3 
and 6 . At least one of the physical weights must be non-zero. 

The user is cautioned to avoid doubling the effect of the rigid body 
modes, as would happen if the user included the "soft spring" rigid body modes 
for the vertical and horizontal plane subsystems as obtained from the NASTRAN 
program and then included the same rigid body modes (that is, rigid body modes 
for motion in directions Z, 6y, and 6 z) in the rigid body subsystem. How- 
ever, the user could model rigid body motion for the Z, By, Y, and $2 direc- 
tions using the vertical and horizontal plane subsystems (which cover those 
four degrees of freedom only) and then use the rigid body subsystem just for 
motion in the X and 8 ^ (Subsystem 9 only) directions (the user would zero out 
or omit the physical weight properties in all but the X and By directions for 
the rigid body subsystem). 

From the above, it is evident that the user has the choice of modeling 
rigid body motion for the Z, By, Y, and 62 directions using either "soft 
spring" rigid body modes in the vertical and horizontal plane subsystems or 
the rigid body subsystem. Both ways of modeling rigid body motion for these 
directions give equivalent results. However, it has been found that TETRA 
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runs faster (and thus saves on cost) if the rigid body modes for these direc- 
tions are included in the vertical and horizontal plane subsystems rather 
than the rigid body subsystem. This is due to the fact that the rigid body 
subsystems have more degrees of freedom (five for Subsystems 3 and 6 and six 
for Subsystem 9) than the vertical and horizontal plane subsystems (which 
have two degrees of freedom). Also, using the vertical and horizontal plane 
subsystems for the rigid body modes eliminates the need to input the center 
of gravity coordinates and the physical weights. 

The input for the case torsional subsystem (Number 10) (see input sheets 
C-7, C-8, and C-9) and the pylon subsystem (Number 11) (see input sheets 
C-10, C-11, and C-12) is a similar to the input sheets for the vertical 
and horizontal plane subsystems (Sheets C-1, C-2, and C-3) . However, 
instead of using the VH array to input the mode shapes as is done for the 
vertical and horizontal subsystems, the TOR array (see input sheet C-9) 
is used to input the mode shapes for the torsional subsystem and the PYL 
array (see input sheet C-12) is used to input the mode shapes for the 
pylon subsystem. 
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TYPE D1 AND TYPE D2 INPUTS 


Type 1 Physical Connecting Element 
(Generalized Spring-Damper Element) 

This element is associated with two physical points located at arbitrary 
locations in global space. Each of these points is assigned six degrees of 
freedom, three translational displacements (or velocities), and three rota- 
tional displacemencs (or velocities). 

A full complement of stiffness and damping coefficients can be input to 
allow the modeling of fully coupled load paths. The units of the stiffness 

and damping coefficients are: Ibp/in., Ibp/rad, ln.-lb_/ln., ln.-lb_/rad, 

F F 

Ibp-sec/ in. , Ibp-sec/rad , in.-lbp-sec/in. , and in.-lbp-sec/rad. 

Damping can be specified directly via the coefficient input or can be 
specified in terms of a structural Q-factor and a selected frequency. For the 
latter case» the TETRA program computes the damping matrix by multiplying the 
stiffness matrix by the proportionality term l/mQp, where w is the selected 
frequency in radians/sec, and Qp is the Q-factor. The user inputs the fre- 
quency with the units cycles/sec and TETRA converts this to radians/sec. 
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NAMELIST 
Type D-1 


TYPE 1 PHYSICAL CONNECTING ELEMENT 
(GENERALIZED SPRING-DAMPER ELEMENT) 


Z (vtrtical tfirtctlon) 



Y (horizontal dirtctlon) 


X (axial 
dlractlon, 
positive 
forward 


1 (vertical direction) 


Y (horizontal direction) 


X (axial direction, positive forward) 


V 

$EKD 
SLIST3 
ITYPE»1 , 

ILEM= .element number 

I>end J>end 

point number point number 


STIFFNESS MATRIX 

DEFINITION 

I-END 

J-END 

GLOBA 

L DIRECTION DISPLACEMENT 

GLOBAL DIRECTION DISPLACEMENT 

1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 

Z 


Y 


X 


Z 

*y 

Y 


X 



Global 

Dir. Forces 


1 

F, 


2 


O 

3 


as 

LU 

4 


1 

V-4 

5 



6 

f’^ 

ex 


1 

F, 


2 

F^ 

O 

3 


LU 

1 

4 

F^, 


5 

F*" 


6 

ex 


SPRING(1,1)= 
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Type D-2 


TYPE 1 PHYSICAL CONNECTING ELEMENT 
(GENERALIZED SPRING-DAMPER ELEMENT) (Continued) 


I DAMP* 


1*damp1ng 
0*no damping 


(0 assumed) 


If IDAMP=1 , complete the input for one of the following 
two options: 


Option 1: Tor damping based on Q-factor and selected 

frequency, enter the following: 

q-factor 

QELEM* , 

frequency (hertz) 

QFREQ* . 


Option 2: For damping based on damping matrix definition, 

enter the following: 


DAMPING MATRIX DEFINITION 

I-END 

J-END 

GLOBAL DIRECTION VELOCITY 

GLOBAL DIRECTION VELOCITY 

1 

2 

3 

tm 

5 

6 

1 

2 

3 

4 

5 

6 

• 

Z 

K 

m 

• 

«z 

■ 

m 

H 


Y 

^2 

H 



Global 

Dir. Forces 

DAMP0.1)= 





2 iy 


O 



LU 

t Sz 


1 








- 




O 

3 

§ t • ft ft ft ft ft 

UJ 

1 




5 ^x 



6 F* 

ex 



i 
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TYPE E INPUT 


Type 2 Physical Connecting Element 
(Space Link - Damper Element) 

This element is used to model load paths which have no local moment capa- 
bility. The load paths are associated with two physical points located at 
arbitrary locations in global space. Each of these points is assigned three 
translational degrees of freedom. The user inputs the cross section areas 
(in. 2), and Young's modulus (Ib/in.^), and TETRA uses this information along 
with the coordinates of the two connecting points to calculate the stiffness 
matrix. The damping can be defined either in terms of translational (dashpot) 
damping directed along the axis of the link or by proportional damping. In 
the latter case, the user inputs a structural Q-f actor and a selected fre- 
quency and the TETRA program computes the damping matrix. 


Ill 
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Type E 


TYP; 2 ’>‘!YSTCAL CONNECTING ELEMENT 
(S PACE LINK-DAMPER ELEMENT) 



$LIST3 

ITYPE=2, 

element number 

ILEM= , 

I -end J-end 

point number point number 


2 

Area(in ) 

TAREA= , 

Young's Modulus (psi) 

TY0UNG= , 

1 =damping 
0=no damping 

IDAMP= , (0 assumed) 


(vertical 

direction) 


Space Link (Truss) 


I 

^ Y (horizontal 
direction) 

direction, positive 
forward) 



X (axial 


If IDAMP=1, complete the input for one of the following two options; 

Option 1: For damping based on Q-factor and selected frequency, enter the 

following: 

Q-factor 
QELEM= , 

frequency (hertz) 

QFREQ= . 


Option 2: For damping based on the translational damping coefficient c along 

the axis of the link element, enter the following: 

C (Ib-sec/in) 

TDRATE= 



TYPE F INPUT 


Type 3 Phsycial Connecting Element 
(Rub Element) 

The rub element allows the mathematical modeling of the nonlinear tip 
rub that includes the dead band displacement internal prior to closure between 
the rotor and case (or between rotor and rotor). The user enters the radial 
spring rate K which represents the local case distortion and the blade compli- 
ance, the structural clearance Cq (mils), and damping coefficient C (Ib-sec/in) 
if so desired. This element can be used to model rotor-to-case rubs as well 
as rotor-to-rotor rubs. In the latter case, the I-end must be on the inner 
rotor and the J-end must be on the outer rotor. In the former case, the I-end 
must be on the rotor and the J-end on the case. 
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TYPE 3 PHYSICAL CONNECTING ELEMENT 
(RUB ELEMENT) 



voeaa 1 

band) _i L 

T 


1 

V 

$ END 

$LIST3 

ITYPE=3, 

ILEM= .element number 

I -end J-end 

point number point number 


Note; If rotor-case rub, I-end must be on rotor and J-end must be on case. 

If rotor-rotor rub, I-end must be on inner rotor and J-end must be 
on outer rotor. 

SK= .Case or outer rotor (if rotor-rotor rub) radial spring 

constant K (Ib/in) (becomes active on closure) 

DBAND= _,Radial dead band Eq (niils) between rotor and case (if rotor- 

case rub) or between inner rotor and outer rotor (if rotor- 
rotor rub) 

CC= ,Rub element damping coefficient C (Ib-sec/in) (becomes active 

on closure) 
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TYPE Gl» TYPEG2 and TYPE G3 INPUTS 


Type 4 Physical Connecting Element 
(Engine Support-Links Element) 

The engine support element is a multipoint, multidirection variable geom- 
etry element that provides the capability to model the complex load paths be- 
tween the engine case and the pylon or ground and thus allows for the direct 
modeling of actual engine mount structures. In addition, this element couples 
the flexible and rigid body-centerline modal subsystem that represents the 
engine case to the support links that comprise the mounting system through the 
case flexibility. The case flexibility is described by three flexibility 
rates SKV, SVH and SKA. SKV is the vertical direction case distortion flex- 
ibility rate (in. /lb). 

SVH is the horizontal direction case distortion flexibility rate (in./ 
lb). SKA is the axial direction case distortion flexibility rate (in. /lb). 

These rates are the reciprocals of the case spring rates that represent 
local distortion of the engine case under mount reaction loads. Multipliers, 
defined by the input values AM, BM, CM, roi, EM, and FM, are used to proportion 
the above flexibility rates per the diagram shown on the Type G-2 input sheet. 
If it is desired to cut a load path at a case attach point J or K in a given 
direction, then the applicable multiplier can be set equal to zero. The load 
paths between the engine case and the pylon or ground can be defined with up 
to 8 link elements. In defining these load paths, no more than 4 distinct 
points can be specified on the pylon or ground side. This means that more 

than a single link may be used to connect a case point to a pylon or ground 

point. The coordinates for the points that correspond to the J and K points 
on the case are defined on the Type B input sheet if these points are con- 
nected to links. The damping is defined by a structural Q-factor and a se- 

lected frequency. These data are used by TETRA to compute the damping matrix 
based on stiffness matrix proportionality. 
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Type G-1 

TYPE 4 PHYSICAL CONNECTING ELEMENT (ENGINE SUPPORT-LINKS ELEMENT) 

Some example configurations that can be modeled with the type 4 element: 



Other side attached to 
subsystem or ground 



$END 

$LIST3 

ITYPE=4, 

ILEM= , element number 


I-end J-end 

point number point number 

(engine centerline) (Case point) 


K-end 

point number 
(Case point) 


JT= 


Point I is attached to subsystem. Points J and K can connect to link load paths or 
subsystems or ground. 


SKV= 

, Vertical direction case distortion 
flexibility rate (in/lb) 

SVH= 

, Horizontal direction case distortion 
flexibility rate (in/lb) 

SKA= 

, Axial direction case distortion 
flexibility rate (in/lb) 
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Type G-2 


TYPE 4 PHYSICAL CONNECTING ELEMENTS (ENGINE SUPPORT-LINKS ELEMENT) (Continued) 


Enter multipliers a, b, c, d, e, and f to proportion the spring rates: 


V 

AM= 

BM=' 

CM= 

DM= 

EM= 

FM= 


Rate 

POINT 1 

J 

K 

Kv 

a 

c 

kh 

b 

d 

Ka 

e 

f 


a + c = 1 .0 '> 
b + d = 1.0 < 
e + f = 1.0 y 


This restraint must be followed 
(zero values are permissible). 


Input one line for each link load path (maximum of 8). If no link load paths, 
omit this input. 


Case side point 

Other side point 



number (must 

number (no more 


Young's 

correspond to case 

than 4 distinct 


Modulus 

point J or K) 

point numbers) 

Area (in'^) 

psi 


TLP(1,1)= 
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NAMELIST 
Type G-3 

TYPE 4 PHYSICAL CONNECTING ELEMENT (ENGINE SUPPORT-LINKS ELEMENT) (Continued) 

Enter the following only if you want damping based on Q-factor and selected 
frequency. If no damping desired, omit this input. 

Q-factor 

QELEM= , 

frequency (herz) 

QFREQ= . 
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TYPE HI AND TYPE H2 INPUTS 


Type 5 Physical Connecting Element 
(Uncoupled Point Spring - Damper Element) 

This input allows the connection of two points with a set of uncoupled 
springs and dampers. These spring/dampers are typically used to connect the 
centerlines of beam*like modal subsystems and provide load paths in three 
translational and two rotational directions. Because of the lack of load path 
coupling, good modeling practice infers that the points being connected by 
this element should be coincident in space. Damping can be either specified 
directly via the coefficient input or can be specified in terms of a struc- 
tural Q-factor and a selected frequency. 
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TYPE 5 PHYSICAL CONNECTING ELEMENT 
(UNCOUPLED POINT SPRING-DAMPER ELEMENT) 


Page of 

NAMELIST 
Type H-1 


I and J points are 
coincident . 


^(vertical direction) 



(horizontal 

direction) 


(axial direc- 
tion, positive 
forward) 


Z (vertical direction) 


6 


e 


y (horizontal direction) 


(axial direction, 
positive forward) 


$END 
$LIST3 
ITYPE=5 , 

ILEM= ,element nunber 

I-end J-end 

point number point number 


K (Ib/in) 

X 


Ky(lb/in) 

K (Ib/in) 
/ . in-lb 
*^0y^ rad 

K, 

6z rad 


■) 

•) 


Enter the following spring constants: 

XS = , (0 assumed) 

YS = , (0 assumed) 

ZS = j (0 assumed) 

TYS = , (0 assumed) 

TZS = , (0 assumed) 
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NAMELIST 
Type H-2 


IDAMP= 


1=damping 
0=no damping 


(0 assumed) 


If IDAMP=1 , complete the input for one of the following two 
options: 

Option 1: For damping based on Q-factor and selected 

frequency, enter the following: 


QELEM= 

QFREQ= 


Q-factor 

» 

frequency (hertz) 
» 


r /lb-sec^ 
^"ih ^ 

f. /lb;jec\ 

> ' in ' 

^z ' in ' 
r /in-lb-sec^ 
rad ‘ 
I in-1 b-sec^ 
^ rad ' 


'6z 


Option 2: For damping based on damping coefficient definition, 

enter the following: 


XD = 
YD = 
ZD = 


TYD= 


TZD= 


, (0 assumed) 
, (0 assumed) 
, (0 assumed) 
, (0 assumed) 
, (0 assumed) 
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Restart and Time Integration Input (Type I) 

If the run is a restart run, the user must input the variable ISTART = 1 
(otherwise this variable should be set to zero or omitted). Also, if the run 
is a restart run, an additional input file is required - namely, the restart 
file which was generated by the original run. This restart file is read in on 
file code 22, so the user must assign the restart file to File 22 (@ ASG con- 
trol card) prior to the run. 

The user must choose for this restart time (RTIME input variable) one of 
the times for which output was printed on the original run (these are the only 
times for which restart information was written on the restart file). If the 
user omits the RTIME variable, the program will restart at the final time 
which was printed out on the original run (provided the original run didn't 
terminate prematurely). 

Next, the user must input the time step (DELTA input variable). The time 
step should be made equal to about 1/40 of the smallest period of oscillation. 
The run always begins at time equal zero and then time accumulates. The final 
time (TFINAL input variable) must also be inputted so that the program knows 
when to stop. It is recommended that the user choose TFINAL such that the 
program will do a small number of time steps until the user is sure that his 
or her input is correct, so as to avoid costly no good runs. 

Next, the user must input the print multiple (IPRMUL variable). This 
value governs the number of time steps that get printed out and the number 
that gets written onto the output restart file. If IPRMUL = 100, then one 
out of every 100 time steps computed is printed out and is written onto the 
restart file. Similarly, the plot multiple (IPpiUL) governs the number of 
time steps that get written onto the output plot file. If IPLMUL = 10, 
then one out of every 10 time steps is written onto the plot file. 

Example: If DELTA = 0.0001 second, TFINAL = 0.1 second, IPRMUL = 100, 

and IPLMUL = 5, then computations are made for 1001 times (starting with 
time = 0 and ending with time = 0.1 second). Eleven times are printed and 
written onto the restart file (0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 
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0.08, 0.09, and 0.1), and 201 times are written on the output plot file (start- 
ing at time = 0 and ending at time = 0.1). 

It is recommended that the user pick IPRMUL such that no more than about 
10 times are printed, in order to avoid being buried in printed output. The 
user should pick IPLMUL such that enough times are included to adequately 
define the curve being plotted (several hundred may be needed). Of course, 
the larger the value of IPLMUL the smaller the output plot file will be. 
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NAMELIST 
Type I 


RESTART AND TIME INTEGRATION 


$ END 
$LIST4 

l=Restart Run 
0=New run 

ISTART= , (0 assumed) 

If run is a restart run (ISTART=1), an input restart file (file code 22) 
must be supplied. 

If restart run, enter the restart time. This time must correspond to one of 
the times for which output was printed on the initial run. 

RTIME= , Restart time (Program assumes the final time printed 

for the initial run provided that the initial run 
didn't terminate prematurely). 


DELTA= 

, Time step (seconds) 

TFINAL 

, Final time (seconds) 

IPRMUL= 

, Print multiple 

IPLMUL= 

, Plot multiple 


$END (Include if this is the last card of a run, otherwise omit.) 
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Rotor Speed and Rate Input (Type J) 

Rotor speed and rate input is required if unbalance forces or gyroscopic 
forces are present. If imbalance forces and gyroscopic forces are not present, 
the rotor speed and rate Input is not needed. If rotor speed and rate input 
is not desired, the user can skip the Type J input sheet altogether. 

If rotor speed and rate input is desired, the user must specify which 
rotor is the ''independent” rotor (rotor for which the speed and rate history 

is specified) by setting IROTI to 1 or 2. (Rotor 1 corresponds to subsystems 
1, 2, and 3 and Rotor 2 corresponds to subsystems 4, 5 and 6). Then the user 
must enter the beginning time (BEGTIM) and beginning speed (BEGRPM) for the 
first speed segment and the ending time and rate for all the speed segments 
(TRHIS array). If speed and rate input is present, there must be at least 
one and no more than ten speed segments. 

If another rotor is present, the other rotor is referred to as the "de- 
pendent” rotor because its speed is a function of the "independent” rotor 
speed. If a second rotor is present, the user should input the coefficients 
A, B, C, and D relating the dependent rotor speed Y to the independent 
rotor speed X, where 


Y » Ax3 + Bx2 + C X + D 
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NAMELIST 
Type 0 

ROTOR SPEED AND RATE INPUT . 


This sheet is required if unbalance forces or gyroscopic forces are desired. 

Enter independent rotor number (rotor for which below ending time/rate table 
is input). Permissible values are 0 (rotor speed and rate not considered), 

1 (rotor corresponding to subsystems 1,2, and 3), and 2 (rotor corresponding 
to subsystems 4, 5, and 6). 

V 

IR0TI= , (0 assumed) 

If IR0TI=1 or 2 enter the following: 

BE6TIM= , Beginning time (seconds) for the first speed segment 

(applies to independent rotor) 

BE6RPM= , Beginning speed (rpm) for the first speed segment 

(applies to independent rotor) 

If IR0TI=1 or 2 enter the following table in chronological order (applies to 
independent rotor) (maximum of 10 segments): 

Ending Time Rate 

(seconds) (rpm/sec) 


TRHIS(1,1)= 


Segment 1 


Segment 2 

9 9 

Segment 3 

9 9 

Segment 4 

9 9 

Segment 5 

9 9 

Segment 6 

9 9 

Segment 7 

9 9 

Segment 8 

9 9 

Segment 9 

f 9 

Segment 10 

9 9 


If a second rotor is present, input the following coefficients relating 
the second-(dependent) rotor speed Y to the independent rotor speed X, 
where Y=AX'^ + BX'^ + CX + D: 

A= , (0 assumed) 

B= , (0 assumed) 

C= , (0 assumed) 

0= , (0 assumed) 


$END (Include if this is the last card of a run, otherwise omit.) 


126 





Unbalance Load Input (Type K Input Sheet) 

If unbalance load input is not desired, the Type K input sheet can be 
skipped altogether. To have unbalance load input, rotor speed input (Sheet J) 
is needed - otherwise the unbalance load input is ignored. If unbalance loads 
are desired, the unbalance load input (UNBAL array) must be provided. For 
each unbalance '*birth event" four values must be entered - the time of birth 
(when the unbalance load becomes active), the point number on the unbalanced 
rotor (the point must lie on Rotor 1 or Rotor 2), the magnitude of the unbal- 
ance in gm-in. , and the phase angle 0 which gives the position of the unbal- 
ance relative to the horizontal (Y) axis at the time of birth. There can be 
from zero to a maximum of 20 unbalance birth events. 

The unbalance load input is quiete flexible. The time of birth can be 
the same for different birth events. Also, the same point can be referenced 
more than once if desired. In this way one could model a point (say, the 
engine fan) which starts out with a nominal unbalance, and later a very much 
larger unbalance is introduced (fan blade loss). If more than one unbalance 
birth event is specified for the same points, the unbalance loads for each 
birth event are added together to get the total unbalance loads. 

If the run is a restart run and unbalance loads were present in the orig- 
inal run, then the unbalance load input present in the original run for those 
birth events that become active prior to the restart time should be left in 
for the new run. (Otherwise, the unbalance loads would disappear for the new 
run.) Additional birth events that became active after the start of the re- 
start run may be added to the UNBAL array input for a restart run, however. 
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NAMELIST. “ 
Type K 


applied loads 

UNBALANCE LOADS 

(Forward looking aft) 


Coordinate system 
fixed to rotor — 


Z (vertical direction) 



l/j(unbalance vector for 
/ birth event i) 

^•^(Phase angle for birth event i) 


Q (angular displacement of rotor) 
r7777-^ Y (horizontal direction 

^ Fixed global coordinate system 


NOTE - to have unbalance loads, must have rotor speed input 
(type J sheet) - otherwise the unbalance load input is ignored. 

If unbalance loads are desired, fill out the following (maximum of 
20 lines): 


Time of 

Point number on 

Magnitude 

Phase Angle ^ 

birth (seconds) 

unbalanced rotor 

(gm-in) 

(degrees) 


W 

UNBAL(1,1)= 



If restart run and the time of birth Is less than the restart time, then the 
unbalance load continues active for the restart run. Time of birth can be 
the same for different birth events. Also, the same point can be referenced 
more than once if desired. 


Send (include if this Is the last card of the run. otherwise omit.) 
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P cos cat and P sin cat Load Input (Type L Input Sheets) 

If P cos cat and P sin cat loads are not desired the Type L input sheet 
can be skipped altogether. If desired, the user must supply six values for 
each P cos cat or P sin cat load (CS array input) , as shown on the Type L 
input sheet. As for the unbalance load input, the same point can be ref- 
erenced more than once. In this case the loads are added together to get the 
total load. There can be from zero to a maximum of 20 P sin cat loads. 

Definitions are as follows: 

P = Force Amplitude (lb or in. -lb) 

^ = frequency (hertz) 
t = time (seconds) 
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Type L 


APPLIED LOADS 

P*cos (mt) and P*s1n (wt) LOADS 



Y (horizontal 
direction) 


GLOBAL 

DIRECTION 

NUMBER 

DIRECTION, 

1 

Z 

2 

0Y 

3 

1 Y 

4 

0Z 

5 

X 

6 

Ox 


X (axial direction, positive forward) 


If P*cos (u>t) or P*sin (wt) loads are desired, enter the following 
(maximum of 30 lines): 



TYPE 

AMPLITUDE 

FREQUENCY 

GLOBAL 

POINT 

(1=C0S 

P 

Ui 

DIRECTION 

NUMBER 

2=SIN) 

(lb or in-lb) 

(hertz) 

NUMBER 


CS(1,1)= 



Send (include if this is the last card of the run, otherwise omit). 




Time-Force History Loads (Type M-1 and M-2 Input Sheets ) 

If time-force history load input is not desired, sheets M-1 and M-2 can 
be skipped altogether. If desired the user must specify the point number, 
global direction number, and table number for each time-force history load 
using the NTF array (see sheet M-1) (values in the NTF array must be inte- 
gers). The table number must be a value between 1 and 10. There can be from 
zero to a maximum of 30 time-force history loads entered in the NTF array. 

Then, for each table number referenced in the NTF array, the array table 
(1, 1, NT) is required, where NT is the table number referenced in the NTF 
array. If the first time-force pair entry in the table is for a time other 
than zero, then the force is assumed to be zero up to the time of the first 
table entry. If the last time-force pair entry in the table is for a time 
less than the final time of the run, then the force is assumed constant and 
equal to the force for the last table entry for times greater than the last 
table entry time. A table may have only one time-force pair entry, in which 
case the force is zero till the time of the entry and afterward equal 
to the force given in the entry. There can be a maximum of ten time-force 
pairs in each table. 
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Page __ of 
NAMELIST 
Type M-1 


APPLIED LOADS 
TIME-FORCE HISTORY LOADS 


Z (vertical direction) 



X (axial direction, positive forward) 


GLOBAL 


DIRECTION 


NUMBER 

DIRECTION 

1 

Z 

2 

0y 

3 

Y 

4 

6Z 

5 

X 

6 

0X 


If time-force history loads are desired, enter the following (values must be 
integers) (maximum of 30 lines): 







Page of 

NAMELIST 
Type M-2 


APPLIED LOADS 

TIME-FORCE HISTORY LOADS (Continued) 


For each time-force table, enter the following (substitute actual table 
number in place of NT in subscript below) (maximum of 10 tables, and 
maximum of 10 time-force pairs in each table): 


TIME 

FORCE 

(seconds) 

(lb or in-lb) 


TABLE (1,1,NT)= 


$ END (Include if this is the last card of the run, otherwise omit.) 




Gyroscopic Loads - Type N Input 

This input identifies the gyroscopic loading locations. Up to 30 point 
numbers on the rotors can be entered along with the polar mass moment of in- 
tertia (lb-in. 2) values. Since this input models the cross-axis coupling 
forces associated with Coliolis evaluation, both the vertical and horizontal 
subsystems must be included for the rotor(s). 
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Page of 

NAMELIST 
Type N 


GYROSCOPIC LOADS 


Note - to have gyroscopic loads, must have rotor speed input (type 0 sheet) 
otherwise the gyroscopic load input is ignored. 

Note - for gyroscopic loading, both the vertical and horizontal subsystems 
must be included for the rotor(s). 

If gyroscopic loads are desired, enter the following (maximum of 30 lines): 


Point 

Polar Moment 

Number 

of Inertia I 

On Rotor 

(Ib-in'^) P 


V 

GYROd ,1 )= 


$END (Include if this is the last card of the run, otherwise omit.) 




Plot File Input (Type P-1 and P-2 Input Sheets) 


If no output plot file is wanted, the user should set variable IPLOT to 0 
and dispense with the rest of input sheets P-1 and P-2. If, however, IPLOT is 
set to 1 or omitted, then an output plot file will be produced, and the user 
should enter the desired input from pages P-1 and P-2, 

The plot file size would in many cases be excessive if all the data that 
conceivably might be desired were written onto the plot file. The NTF (see 
sheet P-1) and the NEPD (see sheet P-2) input arrays were added for this 
reason. 

The NPD array allows the user to select the point and direction pairs 
for which the coordinates, displacements, velocities, and modal forces are 
to be written onto the output plot file. If the user doesn’t want these 
values for any point and direction pairs, then the NPD array should be 
omitted. A maximum of 50 point and direction pairs may be specified. 

The NEPD array allows the user to select the physical connecting element 
number, point number on this element, and direction number for which the phys- 
ical connecting element forces are to be written onto the output plot file. 
Again, if the user doesn't want any such output, he should omit this array. 

A maximum of 50 element, point and direction trios may be specified. 

In addition, certain other data is always written to the output plot file 
(provided an output plot file is generated). This data includes the time, 
independent and dependent rotor speed (see discussion of Type J input sheet) , 
independent and dependent rotor angular displacement , and the relative dis- 
placement magnitude, clearance, and force magnitude for all rub elements 
(Type 3 physical connecting elements) present (if any) . 
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Page of 

NAMELIST 
Type P-1 


DATA FOR PLOT FILE 


Z (vertical direction) 



Y (horizontal 
direction) 


X (axial direction, positive forward) 


GLOBAL 


DIRECTION 


NUMBER 

DIRECTION 

1 

Z 

2 

0Y 

3 

Y 

4 

6Z 

5 

X 

6 

6X 



0=no plot file 
l=plot file produced 


IPL0T= 


(1 assumed) 


Enter the following points and directions (if any) for which the coordinates, 
displacements, velocities, and modal forces are to be written onto the plot 
file (values must be integer) (maximum of 50 lines): 



GLOBAL 

POINT 

DIRECTION 

NUMBER 

NUMBER 


NPD(1,1)= 
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Page of 

NAMELIST 
Type P-2 


DATA FOR PLOT FILE (Continued) 

Enter the following physical connecting elements, points, and directions 
(if any) for which the physical connecting element forces are to be written 
onto the plot file (values must be integer) (maximum of 50 lines): 


PHYSICAL 



CONNECTING 


GLOBAL 

ELEMENT 

POINT 

DIRECTION 

NUMBER 

NUMBER 

NUMBER 


NEPD(1,1)» 


Send (include if this is the last card of the run, otherwise omit.) 
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6.0 NASTRAN/TETRA INTERFACE: MODAL INPUT GENERATOR 


The TETRA program was written so that its input information is the type 
of information readily available to an engineer. As noted in Section 3.0, the 
required information consists of geometry, frequency, mode shapes and modal 
potential energy for each of the subsystems considered. Depending on the ap- 
plication, these data can be determined experimentally or obtained from an 
analysis program. At the General Electric Company, the VAST program is com- 
monly used to provide the information. This program is specifically tailored 
to analyze vibrations of axisymmetric structures and, as such, it represents 
a logical choice for the generation of this input data. If a finite element 
model is needed to provide the mode shapes, the MASS program, which has been 
developed by General Electric, is used. Both MASS and VAST have postproces- 
sors available to transform the program results into the format required by 
TETRA. These particular programs have been modified to conform to the TETRA 
requirements, but their basic structure has been used for many years at 
General Electric for several postprocessing activities. 

One of the terms of this contract was to code a similar postprocessor for 
NASTRAN level 17.5 as used by NASA on the UNIVAC 1100/42 computer system. 

This program has been written, debugged and successfully tested for the 
TETRA demonstrator case. Instructions on the use of this program are con- 
tained later in this section, but an overall description is given at this 
point . 

The information needed for the TETRA modal subsystem input can be found 
in the NASTRAN output. This fact guided us to write a program which reads 
this output, stores the variables which will be of interest and then, through 
a set of user responses in an interactive mode, the program selects the spe- 
cific data required, formats it correctly, and writes it on another file. 

This program can generate TETRA modal subsystem input for the vertical 
and horizontal plane subsystems (Numbers 1, 2, 4, 5, 7, and 8) (input sheets 
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C-1 , C“2 , and C-3) and the case torsional subsystem (Number 10) (input sheets 
C-7, C-8, and C-9)- Tb« program cannot, however, generate modal subsystem 
input for the rigid body r.ubt;yH*:ero^ USumbtr." 3. 6, and 9) (input sheets C-4, 
C-5, and C-6> or for the pylon Kubsystem (Number 11) (input sheets C-10, C-11, 
and C-12). 

Three fundamental rules were followed during the coding of this program. 

1. It w^as written in time-sharing ASCII for convenience in setting up 
TETRA files in an interactive fashion. Our experience at General 
Electric has shown that this method of operation is far more effi- 
cient than batch programming for this task. Thus, although NASTRAN 
runs may require large coro anc i.lme limits the postprocessor runs 
vil I be short . 

2. Simple, straight-forward FORTAN coding was used throughout. This 
decision was enforced as ai» aid for possible future enhancements. 

As the program is now written engineers with only moderate program- 
ming experience should have no trouble in modifying the routines to 
suit their own needs. There are only two exceptions to this. In 
the first part of the program, the NASTRAN output is read and copies 
directly to another file. This is done only to avoid BACKSPACE 
problems associated with the UNIVAC system. If the user has copied 
this file already, he/she need not copy it again. The second area 
of "awkward logic" is the coding required to read left justified 
integers in the NASTRAN output. This portion of the code reads the 
integers twice, once as an integer and once as a character array and 
actually counts the number of blanks in the field. Outside of these 
items, all coding is straight-forward. 

3. The user response questions were designed to permit a variety of 
applications. For example, although the NASTRAN run might calculate 
20 frequencies, the user might need only Modes 1, 2, and 6. This 
selection can be easily accomplished. 

In summary, the NASTRAN/TSTRA interface program represents a working, 
easily modifiable program based on the approaches we have found most satisfac- 
tory at the General Electric Company. 

The following section summarizes the NASTRAN 17.5 postprocessor operation 
in an interactive mode and also includes a sample computer-operator "conver- 
sation." 
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6.1 NASA 17 . 5 Generated Modal Data File 


Below is a sample interactive computer-operator "conversation" initiated 
by running the NASTRAN/TETRA interface program. From this conversation the 
interface program knows what values to obtain from the NASTRAN 17.5 output 
file, and after obtaining these values the interface program generates a 
file containing TETRA modal subsystem input. The below example demonstrates 
the use of the interface program to generate modal subsystem input for sub- 
system 7 (case vertical plane subsystem) and subsystem 8 (case horizontal 
plane subsystem) for the demonstrator model. Following the below computer- 
operator "conversation" is a listing of the resulting interface program 
output file which consists of TETRA modal subsystem input. 


RUN 


N A S T R A N -- T E T R A INTERFACE. PROGRAM 

SOME HELPFUL INFORMATION: 

THE ORIGINAL NASTRAN OUPUT MUST EE ASSIGNED 
TO FILE CODE 28 

TO ACCOMODATE THE UNI VAC SYSTEM REQUIREMENTS 
THIS FILE MUST EE COPIED TO ANOTHER .ASCII FILE 
THIS FUNCTION CAN EE ACCOMPLISHED BY THIS PROGRAM 
BY RESPONDING WITH A ONE (i) TO THE QUESTION CONCERNING 
WHETHER OR NOT YOU DESIRE AN OUTPUT READ/WRITE 

IF YOU EXERCISE THE READ/WRITE OPTION THE 
ORIGINAL NASTRAN OUTPUT FILE WILL EE COPIED TO 
FILE CODE 27 

THIS FILE CAN BE SAVED OFF ONTO A PERMANENT 
FILE FOR FUTURE USE 

ONCE THE FILE HAS BEEN SAVED THERE IS NO REQUIREMENT 
FOR FUTURE READ/WRITE(S) - WHEN FURTHER PROCESSING 
OF THE NASTRAN OUTPUT FILE IS DESIRED 

ALL OUTPUT FROM THIS INTERFACE PROGRAM IS ON 
TEMPORARY FILE CODE 29 
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u u 


LUCK 


TYPE 0 TO SKIP NASTRAN OUTPUT READ/WRITE 
> 1 . 

INPUT ? 

SUBSYSTEM NUMBER 
NUMBER OF FREQUENCIES DESIRED 
‘ J "■ h A L- T U R 
^ / ? o ? 1 o 

INPUT FREQUENCY NUMBERS DESIRED 
^ ^ ? A 

INPUT NUMBER OF GRID POINTS YOU WISH TO ELIMINATE 

>7 

j.N='iJT THE GRID POINTS TO EE ELIMINATED 


> 34 j 3 

5 1 1 W0 > i 0 3 t 

3. 0£ 1 1 03 ? 9000 

I NPUT 

IDENTIFIC 

ATI ON TITLE 

> case: 

VERTICAL 

PLANE SUBSY: 

i Hz.Kz. 

WERE 5 

JO I NTS 

THERE 

WERE 3 

EIGENVALUES 

THERE 

WERE 3 

EIGENVECTOR: 

WR I TT 

E 1 V T 0 THE 

TETRA INPUT i 


type i TO GENERATE ANOTHER SUBSYSTEM INPUT FILE 

>1 

INPUT: 

SUBSYSTEM NUMBER 
N*JME»t.R OF FREQUENCIES DESIRED 
Q-FACTOR 
^ U ? -U ? I 

input FREQUENCY NUMBERS DESIRED 
^ Z y -»* ? 4 

INPUT NU MEER OF GRID POINTS YOU WISH TO ELIMINATE 
>7 

INPUT THE GRID POINTS TO BE ELIMINATED 
^ SA ? ? i. 00 ? 1 0 1 1 i 0il ? 1 0S» ? V000 

INPUT IDENTIFICATION TITLE 

>CASE HORIZONTAL PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL 

there WERE 5 JOINTS 

THERE WERE 3 EIGENVALUES 

THERE WERE 3 EIGENVECTORS 

WRITTEN TO THE TETRA INPUT FILE . 

type i TO GENERATE ANOTHER SUBSYSTEM INPUT FILE 


foLLo\^it^o- /s fi of rue OUTPUT Fii-e fpo>^ -rms 

imTCRFACX Pfi/OCrP^m (Flt^ c6 t>€ : 

SEND 
$ L I S T Z 

TITLE= m:ASE vertical PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL 
ISUB= 7? 

X R z. P ^ 0 y 
Y R fc. F = 0 f 
Z! F»* t. F ^ 0 ? 

PTS(iyl)= 
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1 f 

0. » 

0. 

f 

0. 


z » 

-10.00001 

0. 

9 

0. 


3 9 

-90.00009 

0. 

9 

0. 


33 f 

-100,0000. 

0. 

9 

0. 


37 f 

-50. 0000 f 

0. 

t 

0. 



»;mode£(i . 1) = 


1 .2515911E 

0Z 9 

9.999440E 

01f 15f 

1. 

1 .588929E 

02. 

8.399970E 

01f 15f 

1. 

3.084352E 

04 9 

2.379934E 

06.15f 

0. 

VH( 1.1.1) 

= 





9 . •5’995S!;IE- 0 1 . 6 . 1 5 18 1 3E-07 . -4 . 999698E 00 » -6 . 1 82 1 94E-04 r 
9.999628E-01 > 6.235140E-07. 2.999931E 01t 9.999148E 01» 
9.99962iE-01.T8.431300E-07»-2.999929E 01» 9.999298E 01. 
9.999570E-0! .-6.347973E-07. 4.99969SE 00. 0. . 

1.000000E 00.-9.808070E-09. 3.039837E-06.-1 .499972E 03. 

V»H<1 . 1 .2) = 

-9.999996E-01 . 2.000014E-02. 8.087115E 00. -4 .571877E 01. 
-7.999982E-01; 2 .00001 2E-02 . - 1 . 093825E 01 . -2 . 9S4927E 02. 
7.999966E-01 . 2.000012E-02.-1.093826E 01. 2.070593E 02. 
1.000000E 00. 2 . 0000 1 4E-02 . 8.067120E 00. 0. . 

1 .9^;5323E-07. 2.000014E-02. 8. 7154£0E-01 . -4 . 571713E 01. 
vH(l?1.3)= 

1.000000E 00.-3. i21S30E-02.-3.000452E 05. 2.454122E 08. 

8. 1 i6338E-0i.-3..027587E-02.-9.671252E 05. 1.3S0302E 07. 
8.116836E-01 . 3.027587E-02. 9.871252E 05. 8.855028E 08. 
•.000000E 00. 3. 121830E-02. 3.000452E 05. 0. . 

-7. 185308E-01.-8.384795E-15. 0. . 8.S02237E 07. 

*END 

»LI3T2 

TITLE= 'CASE HORIZONTAL PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL 
ISUB= 8. 

XREF=0, 

YREF=0. 

ZREF=0, 

PTS (1.1)= 

1. 0. » 0. f 0. » 

2 . -10. 0000 » 0 . t 0 . . 

Si. -90.0000. 0. t 0. t 

33. -100.0000. 0. . 0. f 

37. -50.0000? 0. f 0. . 

X MODES «: 1 . 1 ) = 

i. 25091 IE 02. 9.999440E 01.15. 1. 

1.5S8929E 02. 6.399970E 01.15. 1. 

3.084352E 04. 2.379934E 08.15. 0. 

< 1 , 1 , 1 ) = 

9 . 999530E-01 . 8.15181 3E-07 • -4 . 999898E 00.-8.1821 94E-04 . 
9.999828E-01 . 8.235140E-07. 2.9999:I:iE 01. 9.999148E 01. 

7 . 9 Wrtic! i E-0 1 . -6 . 43 i SI00E-07 ? - 2. 9999Z9E 01. 9.99929SE 01. 

9 ^ 999570E-01 ? -8 . 347973E-07. 4.999899E 00 • 0. . 

1 , 000000E 00 » -9 . 303070E-09 . 3. 039S37E-08 . - 1 . 499972E 03 ? 

■3H< 1.1.2) = 

-9,999v98E-0! ? 2.000014E-02. 8.087115E 00 . -4 . 571877E 01. 
-7.999982E-01 ? 2.000012E-02. -1 .095825E 01 . -2 . 9S4927E 02. 
7.999988E-01 . 2 . 00001 2E-02 . - 1 . 095828E 01. 2.070593E 02. 
■.000000E 00. 2. 000014E-02. 8.087120E 00. 0. . 

1 . 945323E-07. 2.000014E-02. 8.715440E-01 .-4.571713E 01. 
vH(l . 1,3)= 


0 0 0 i'i 0 0 E v) 0 f 

-3.1Z1S30E-02. 

-3.00045ZE 

05. 

Z.65412ZE 

08 

116S36E-01, 

-3 . 027587E-02 ? 

-9.67iZ5ZE 

05. 

1 .33030ZE 

07 

1 1 68 - 0 1 4 

3 • 0Z75S7E-0Z 9 

9.67iZ5ZE 

05 f 

S.655026E 

08 

000000E 009 

3. 12l830E-0Z» 

3.00045ZE 

05 f 

0. 


1EJ5308E-01 9 

-S.364795E-15f 

0. 

t 

6.S0ZZ37E 

07 
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7 • 0 DEMONSTRATOR CASE 


Input Setup 

The results obtained with a demonstrator model that represents the engine 
system shown in Figure 54 will now be presented. The straight lines shown for 
the rotor and case in this figure represent segmented beam centerline elements 
that include both stiffness and mass properties. Free**free modal data for the 
rotor and case subsystems along with physical spring-damper connecting element 
data were input to the TETRA time transient computer program and solutions 
were obtained for various rotor speed and rotor imbalance conditions. 

Modal data based on planar finite element models of the total system 
shown in Figure 54 and on the free-free rotor and casing subsystems was com- 
puted with the VAST program. In the former case, a frequency domain modal 
analysis was used to establish the steady-state frequency response of the 
total system and in the latter case, the TETRA program was used to synthesize 
the free-free modal data and the physical connecting element data and to pre- 
dict the time domain response for the combined system. 

Spectral analysis and a comparison of the TETRA and frequency domain 
solutions indicated that TETRA had correctly synthesized the modal and physi- 
cal data to predict the time transient response for the combined system. 

The frequencies for the free-free subsystem modes used in the TETRA anal- 
ysis are shown in Table III. Eight generalized coordinates, corresponding to 
these frequencies, were used to define the time response for the total system, 
which has 40 degrees-of-freedom. Therefore, the dimensionality of the model 
has been reduced by 80%. 

Inspection of Table III shows that the smallest period of oscillation 
for the subsystem modes is equal to 1.0/510.62 = 0.001958 sec. Accounting 
for the stiffening effect of the physical connecting springs, a time step 
value that is considerably less than 0.001958 second must be used for the 
numerical integration that will be performed with the TETRA program. In prac- 
tice, the time step should be made equal to about 1/40 of the smallest period 
of oscillation for the subsystem modes. 

1/40 X .001958 = 4,895 x 10“5 
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(GYRO STATION) (GYRO STATION) 


PHYSICAL CONNECTING ELEMENTS 
Spring 1 0.5 X 10® Ib/in, 

Spring 2 0.5 x 10® Ib/in. 

Spring 3 1.0 x 10® Ib/in. 

Spring 4 1.0 x 10® Ib/in. 


POINT 1 IS THE GLOBAL ORIGIN 
(X=Y=2=0) FOR THE TETRA MODEL. 

PHYSICAL POINT NUMBERS 
ETC. IDENTIFIED. 


Figure 54, 40 DOF Finite Element Model for Demonstrator. 




On this basis a time step of = 50 x 10*"^ seconds was selected. 

Table III, Frequencies* for the Free-Free Modal 

Subsystems Used in the TETRA Analysis. 



Rotor 

Case 


cpm 

Hz 

cpm 

Hz 

99.3 

1.655 

125.1 

2.085 

120.8 

2.013 

159 

2.65 

10065 

167.75 

30367 

510.62 

16851 

280.85 



26425 

440.42 




* Free-Free Modal Data Obtained with the Following 
Boundary Conditions: 

Rotor was Supported on Each End by 300 Ib/in. Spring and Case 
was Supported on Each End by 100 Ib/in. Spring. 

The free-free modal data were used in the TETRA program to respresent 
both the vertical and horizontal planes and uncoupled springs (in the Z and 
Y directions) were used to connect the modal subsystems to each other and to 
the ground to model the configuration shown in Figure 54. The rigid body 
modes for the rotor and the case were approximated with the "soft spring" 
rigid body modes at 1.655, 2.013, 2.085, and 2.65 Hz. In addition, they were 
also defined in subsequent analyses with the true representations that were 
based on the mass properties. The results obtained with either approach were 
in excellent agreement. A modal Q-factor equal to 15 (3.33% critical damping) 
was used to represent the damping of each of the casing subsystem modes. The 
rotor subsystem modes were undamped. Proportional damping based on a physical 
Q-factor equal to 15 and a selected frequency of 3624 rpm (a critical speed 
with gyroscopic stiffening present) was used to represent the damping for the 
connecting spring elements. 

7.1 INPUT AND OUTPUT 

The TETRA input and representative output follows for constant 3000 rpm 
speed running with 5000 gm-in. sudden fan unbalance and a 10 mil rub element 
at the fan. The input listing is a part of the TETRA output ai\d directly fol- 
lows the cover sheet of the output. 
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TTTTTTTTT 

EEEEEEE 

TTTTTTTTT 

RRRRRRRR 

AAAAAA 

TTTTTTTTT 

EEEEEEE 

TTTTTTTTT 

RRRRRRRRR 

AAAAAAAA 

TTTTTTTTT 

EEE 

TTTTTTTTT 

RRR 

RRR 

AAA 

AAA 

TTT 

EEE 

TTT 

RRRRRRRRR 

AAA 

AAA 

TTT 

EEEEEEE 

TTT 

RRRRRRRR 

AAAAAAAAAAAA 

TTT 

EEEEEEE 

TTT 

RRRRRRR 

AAAAAAAAAAAA 

TTT 
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TTT 

RRRRRRRR 
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AAA 

TTT 
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TTT 

RRR RRRRR 
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TTT 
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TTT 

RRR 
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RRR 

RRRR 
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AAA 


G. BLACK DATE AUGUST 1980 CHARGE 19510 

BLDG. 500 MAIL DROP K-190 EXT. 3334 

ROTOR AND CASE VERTICAL AND HORIZONTAL MODEL 1 

WITH GYRO LOADS AT FAN AND TURBINE 

DATE 01/05/81 
TIME 16.94 
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( ( 


. 0959 , - . 00459 , 1 . 24 1 E6 , 2 . 250E7 . 
.6966, .00594.6.883E6, -2.488E5. 
1.0. -2.426E5. -9.743E6. 


. 6966 . - .00594 . - 1 . 026E6 . 1 . 768E4 . 
.0959. .00459,-1.241E6,2.250E7, 
-. 10361. .01856.0.0. 

$£MD 

$LIST2 

TITLE® ’rotor horizontal PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL l', 
ISUB=2. 

PTS( 1.1)= 

4.0. 0,0, 

5, -10.0,0. 

34.-30.0.0. 

38.-50.0.0. 

35, -70.0.0. 
6. -90.0.0. 

36, -100.0.0. 
XM0DES(1.1)= 

99.3. 192.0. 1. 

120.8.300.0. 1. 
10065.5.38E5.0.0. 
16851 ,9-90E6.0,0. 

26425, 1 .75E6,0,0, 
VH( 1. 1. D- 
-1. .02,0,0. 

-.8. .02,0,0, 

-.4, .02.0.0. 
0. .02,0.0, 
.4, .02.0.0. 
.8. .02.0,0, 

1. .02,0.0, 
VH( 1,1,2)= 
1.0. 0.0. 
1.0, 0.0. 

1.0, 0.0. 
1 .0.0.0. 
1 .0.0.0. 
1 .0.0.0. 

1,0. 0.0. 

VH( 1.1,3)= 

. 12215, - .034677,0,0. 

-.25282. -.029534, -2. 073E5, 1. 176E7, 

- .80295, -.016102, -1 .02 1E5, 1 .568E7, 
-1,0.3.519E4, 1.711E7. 

-.80295, .016102. 1 .586E5, 1 .558E7, 
-.25282, .029534, 2. 073E5, 1 . 176E7, 

. 12215. .034677,0,0. 

VH( 1,1.4)= 

.2748. . 1075.0.0, 

. 8258 . . 0707 . - 1 . 444E6 . -6 . 956E7 . 

.8584, .0174. -1 .892E6, -3.877E7, 
0, - .0017, - 1 .989E6, 1 .608E4, 
-.8584, .0174.-1 .723E6.3.846E7, 
- . 8258 . . 0707 , - 1 . 444E6 , 6 . 956E7 . 

- .2748, . 1075.0,0, 

VH( 1 , 1,5)= 

-. 10361 .- .01856,0.0. 

.0959. -.00459, 1 . 24 1 E6 . 2 . 250E7 , 


.6966, .00594.6.883E6.-2.488E5. 
1,0.-2.426E5.-9.743E6, 

. 6966 . - . 00594 , - 1 . 026E6 . 1 . 768E4 , 



8 







*CASE VERTICAL PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL tV 


. 0959 . . 00459 . - 1 . 24 1 E6 , 2 . 250E7 , 

10361, .01856,0,0, 

SEND 


SLIST2 

TITLE* CASE VERTICAL PLANE SUBS 
ISUB=7, 

PTSd. 1) = 


1.0. 0,0. 

2 . - 10 . 0 . 0 . 
37.-50.0.0. 
3.-90, 0.0, 


33. - 100.0.0. 
XMOOESd.D- 

125. 1.100.0. 1. 

159.64.0. 1. 


30637.2. 38E6. 15.0, 

VH(1. 1. D* 

1 . 0 , 0 . 0 , 

1 . 0 . 0 . 0 . • 


1 . 0 , 0 . 0 . 

1 . 0 . 0 . 0 . 

1 . 0 . 0 . 0 . 

VM(1.1,2)« 


- 1 . . 02 , 0 , 0 . 

-. 8 , . 02 . 0 , 0 , 

O. ,02.0,0. 

. 8 . . 02 , 0 , 0 . 


1 , . 02 . 0 . 0 , 

VHd.1,3)>* 

1. -.03122,0.0. 

.6117.-.03028.-9.6716E5, 1. 1229E7, 


- .7189.0,4.3121E5,6.8026E7. 
.6117, .03028,9.6716E5, 1. 1229E7, 

1. .03122.0.0. 

SEND 


$LISt2 , 

TITLE- CASE HORIZONTAL PLANE SU 
I5UB*8. 

PTSd.D* 


1 , 0 . 0 . 0 . 

2 . - 10 . 0 . 0 . 

37.-50.0,0. 

3. -90.0.0, 


33.-100,0.0, 

XMODES( 1,1)= 

125. 1. 100,0. 1. 

159,64.0d^ 


30637, 2. 38E6, 15,0, 

VHd,1.l)« 

1 . 0 . 0 . 0 . 

1 . 0 , 0 , 0 , 


1 . 0 . 0 . 0 . 

1 . 0 . 0 . 0 . 

1.0. 0.0, 

VH(1,1.2)= 


- 1 . . 02 , 0 , 0 . 
-. 8 . . 02 . 0 . 0 . 
0 . . 02 . 0 , 0 . 
. 8 . . 02 , 0 . 0 , 


1 , . 02 . 0 , 0 . 

VH( 1 . 1 .3)- 
1 . - .03122,0.0. 


CASE HORIZONTAL PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL 1^ 



( 



.61 17. - .03028. "9.67 16E5. 1 . 1229E7. 
-.7185.0.4.3121E5.6.8026E7. 

■6117, .03028,9.6716E5, 1 . 1229E7. 

1 , ,03122,0.0. 

$£ND 

$LIST3 

ITYPE=5. 

ILEM=1 , 

JT=2,39. 

YS=5E5. 

ZS==5E5. 

IDAMP=1. 

0ELEM=15. 

OFREO*60.4. 

$Cr^D 

$LIST3 

ITYPE=5. 

ILEM*2. 

JT=3,40, 

YS=5E5, 

ZS=5E5. 

IDAMP=1. 

QELEM=15. 

OFREO=60.4. 

$END 

SLIST3 

ITYPE°5. 

ILEM«3. 

JT«5,2, 

YS«1E6, 

2S°1E6, 

IDAMP«1 , 

0ELEM«15. 

0FRE0“60.4. 

$EH0 

$LIST3 

ITYPE=5. 

ILEM»4. 

JT°6,3, 

YS=1E6, 

ZS=*1E6. 

IDAMP«1 . 

QELEM^IS, 

0FRE0=60.4. 

$LIST3 

1TYPE = 3, 

ILEM«5, 

JT=4. 1 , 

SK=1E6. 

DBAND°10. 

CC»0.0, 

$CND 

$LIST4 

DELTA=^. 00005, 

TFINAL=.512, 

IPRMUL=1000, 

IPLMUL^IO, 

IROTI^I . 

BEGTIM=0. , 

BEGRPM=3000. . 

TRHIS(1.1)» 
















152 






154 








w CO colco CO CO ^ ^ ^ ^ ^ ^ Ito in in tnjin in in 


( ( 


3 . 

5 

-0.25282 

-0.02953 

-2.073E 05 

3 

34 

-0.80295 

-0.01610 

-1.021E 05 

3 

38 

-1.00000 

0. 

3.519E 04 

3 

35 

-0.80295 

0.01610 

1.586E 05 

3 

6 

-0. 25282 

0.02953 

2.073E 05 

3 

36 

0. 12215 

0.03468 

0. 

4 

4 

0.27480 

0. 10750 

0. 

4 

5 

0.82580 

0.07070 

-1.444E 06 

4 

34 

0.85840 

0.01740 

-1.892E 06 

4 

38 

0. 

-0.00170 

-1.989E 06 

4 

35 

-0.85840 

0.01740 

-1.723E 06 

4 

6 

-0.82580 

0.07070 

-1 .444E 06 

4 

36 

-0.27480 

0. 10750 

0. 

5 

4 

-0. 10361 

-0.01856 

0. 

5 

5 

0.09590 

-0.00459 

1.241E 06 

5 

34 

0.69660 

0.00594 

6.883E 06 

5 

38 

1.00000 

0. 

-2.426E 05 

5 

35 

0.69660 

-0.00594 

-1.026E 06 

5 

6 

0.09590 

0.00459 

-1.241E 06 

5 

36 

-0. 10361 

0.01856 

0. 


I 


CM O 



-6.956E 07 
“3.877E 07 
1.608E 04 


3.846E 07 
6.956E 07 
O. 
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,} 

DATA FOR MODAL SUBSYSTEM 2 


ROTOR HORIZONTAL PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL 1 


NUMBER 

GLOBAL 

OF SUBSYSTEM DIRECTIONS* 2 
DIRECTIONS FOR SUBSYSTEM* 3 4 








COORDINATES OF REFERENCE POINT 

RELATIVE 

TO GLOBAL SYSTEM (IN.) 





X= 0. 

Y= 0. 


Z= 0 








COORDINATES 

OF POINTS 

ON SUBSYSTEM (INCHES) 




■ *• - • - • 

POINT 

LOCAL COORDINATE SYSTEM 

GLOBAL COORDINATE SYSTEM 



NUMBER 

X 

Y 

Z 

X 

Y 


z 




4 

0. 

0. 

0. 

0. 

0. 


0. 




5 

-10.000 

0. 

0. 

-10.000 

0. 


0. 




34 

-30. OCX) 

0. 

0. 

-30.000 

0. 


0. 




38 

-50.000 

0. 

0. 

-50.000 

0. 


0. 




35 

-70.000 

0. 

0. 

-70.000 

0. 


o. 




6 

-90.000 

0. 

0. 

-90.000 

0. 


0. 




36 

-100.000 

0. 

0. 

-100.000 

0. 


0. 




NUMBER 

OF SUBSYSTEM 

POINTS® 7 









LOCAL 

GENERALIZED 




MODE TYPE 



GENERALIZED 

GENERALIZED 

MODE 

COORDINATE 

FREQUENCY 

POTENTIAL 0 

0«FLEXIBLE 

GENERALIZED 

STIFFNESS 

DAMPING VALUE 

NUMBER 

NUMBER 

RPM 

ENERGY FACTOR 

1»RIGID BODY 

WEIGHT-LB 

LB/IN 

(LB'SEO/IN 

1 

6 

99. 

1 .920E 

02 0. 


1 


1.372E 03 

0. 

0. 

2 

7 

121 . 

3.000E 

02 0. 


1 


1.449E 03 

0. 

0. 

3 

8 

10065. 

5 . 380E 

05 0. 


0 


3.743E 02 

1.076E 06 

0. 

4 

9 

16851 . 

9.900E 

06 0. 


0 


2.457E 03 

1.980E 07 

0. 

5 

10 

26425. 

1 . 750E 

06 0. 


0 


1.766E 02 

3.5CK)E 06 

0. 

NUMBER 

OF SUBSYSTEM 

MODES* 5 









THE MODE SHAPES AND 

CORRESPONDING FORCES FOR THIS SUBSYSTEM ARE 

- 




(SIGNS 

IN THETA-2 DIRECTION CHANGED TO 

OBTAIN RIGHT 

HAND COORDINATE SYSTEM) 






MODAL DISPLACEMENTS 


MODAL 

FORCES 



LOCAL 

GENERALIZED 


GLOBAL 

DIRECTION 


GLOBAL 

DIRECTION 




COORDINATE 

POINT 

3 

4 


3 


4 



NUMBER 

NUMBER 

NUMBER 

Y 

THETA-Z 


Y 


THETA-Z 




6 

4 

-1.00000 

-0.02000 


0. 


0. 



1 

6 

5 

-0.80000 

-0.02000 


0. 


0. 



1 

6 

34 

-0.40000 

-0.02000 


0. 


0. 



1 

6 

38 

0. 

-0.02000 


0. 


0. 



1 

6 

35 

0.40000 

-0.02000 


0. 


0. 



1 

6 

6 

0 . 80000 

-0.02000 


0. 


0. 



1 

6 

36 

1 . 00000 

-0.02000 


0. 


0. 



2 

7 

4 

1.00000 

0. 


0. 


0. 



2 

7 

5 

1.00000 

0. 


0. 


0. 



2 

7 

34 

1 .00000 

0. 


0. 


0. 



2 

7 

38 

1.00000 

0. 


0. 


0. 



2 

7 

35 

1 .00000 

0. 


0. 


0. 



2 

7 

6 

1.00000 

0. 


0. 


0. 



2 

7 

36 

1 .ooooo 

0. 


0. 


0. 
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(■ 


( 


0.03468 O. O. 

0.02953 “2.073E 05 -1.176E 07 

0.016 10 -1.02 IE 05 -1.568E 07 

0. 3.519E 04 -1 .71 IE 07 

“0.01610 1.586E 05 “1.558E 07 

-0.02953 2.073E 05 - 1 . 176E 07 

-0.03468 0. 

-0. 10750 0. O. 

-0.07070 -1.444E 06 6.956E 07 

-0.01740 “1.892E 06 3.877E 07 

0.00170 -1.989E 06 -1.608E 04 

-0.01740 -1.723E 06 -3.846E 07 

-0.07070 -1.444E 06 -6.956E 07 

-0.10750 O. q, 

0.01856 0. 0. 

0.00459 1.241E 06 -2.250E 07 

-0.00594 6.883E 06 2.488E 05 

0. -2.426E 05 9.743E 06 

0.00594 -1.026E 06 - 1 . 768E 04 

-0.00459 “1.241E 06 -2.250E 07 

-0.01856 O. q. 
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DATA FOR MODAL SUBSYSTEM 7 

CASE VERTICAL PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL 1 


NUMBER OF SUBSYSTEM DIRECTIONS^ 2 
GLOBAL DIRECTIONS FOR SUBSYSTEM* 1 2 


COORDINATES OF REFERENCE POINT RELATIVE TO GLOBAL SYSTEM (IN.) 
X* 0. Y= 0. Z= O. 


POINT 

NUMBER 



COORDINATES OF POINTS ON SUBSYSTEM (INCHES 
LOCAL COORDINATE SYSTEM GLOBAL COORDINATE SYSTEM 

X Y Z X Y Z 


00 

00 

00 


- 100.000 



- 100.000 


NUMBER OF SUBSYSTEM POINTS* 


LOCAL 

MODE 

NUMBER 


GENERALIZED 

COORDINATE 

NUMBER 


FREQUENCY POTENTIAL 
RPM ENERGY 


1 

1 1 

125. 

1.000E 02 


2 

12 

159. 

6.400E 01 


3 

13 

30637. 

2.380E 06 

1 


MODE TYPE 


GENERALIZED 

GENERALIZED 

0*FLEXIBLE 

GENERALIZED 

STIFFNESS 

DAMPING VALUE 

R 1»RIGID BODY 

WEIGHT -LB 

LB/ IN 

(LB-SEO/IN 

1 

4.503E 02 

0. 

0. 

1 

1.784E 02 

0. 

0. 

0 

1.787E 02 

4.760E 06 

9.891E 01 


NUMBER OF SUBSYSTEM MODES* 3 

THE MODE SHAPES AND CORRESPONDING FORCES FOR THIS SUBSYSTEM ARE- 


LOCAL GENERALIZED 
MODE COORDINATE POINT 


MODAL DISPLACEMENTS 
GLOBAL DIRECTION 


MODAL FORCES 
GLOBAL DIRECTION 


NUMBER 

NUMBER 

NUMBER 

Z 

THETA-Y 

z 

THETA-Y 

1 

11 

1 

1.00000 

0. 

0. 

0. 

1 

1 1 

2 

1 . 00000 

0. 

0. 

0. 


1 

1 1 

37 

1.00000 

0 

0. 

0. 

1 

1 1 

3 

1 . 00000 

0 

0. 

0. 

1 

1 1 

33 

1.00000 

I 

0. 

0. 


12 

1 

- 1 . 00000 

0.02000 

0. 

12 

2 

-0.80000 

0.02000 

0. 

12 

37 

0. 

0.02000 

0. 

12 

3 

0.80000 

0.02000 

0. 



2 

12 

33 

1 . 00000 

0.02000 

0. 

0. 

3 

13 

1 

1 . 00000 

-0.03122 

0. 

0. 

3 

13 

2 

0.61170 

-0.03028 

-9.672E 05 

1. 123E 07 

3 

13 

37 

-0.71850 

0. 

4.312E 05 

6.803E 07 

3 

13 

3 

0.61170 

0.03028 

9.672E 05 

1 . 123E 07 

3 

13 

33 

1 . 00000 

0.03122 

0. 

0. 
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DATA FOR MODAL SUBSYSTEM 8 


CASE HORIZONTAL PLANE SUBSYSTEM FOR DEMONSTRATOR MODEL 1 


NUMBER 

OF SUBSYSTEM 

DIRECTIONS 

= 2 








GLOBAL 

DIRECTIONS FOR SUBSYSTEM* 3 4 








COORDINATES OF REFERENCE POINT 

RELATIVE TO GLOBAL SYSTEM (IN.) 






x= 

3. 

Y= 0. 


z= 0 . 








COORDINATES 

OF POINTS ON SUBSYSTEM (INCHES) 






POINT 

LOCAL COORDINATE SYSTEM 

GLOBAL COORDINATE SYSTEM 





NUMBER 

X 

Y 

Z 

X 

Y 

Z 





1 

0. 

0. 

0. 

0. 

0. 

0. 





2 

-10.000 

0. 

0. 

-10.000 

0. 

0. 





37 

-50.000 

0. 

0. 

-50.000 

0. 

0. 





3 

-90.000 

0. 

0. 

-90.000 

0. 

0. 





33 

-100.000 

0. 

0. 

100.000 

0. 

0. 





NUMBER 

OF SUBSYSTEM 

POINTS* 5 









LOCAL 

GENERALIZED 




MODE TYPE 



GENERALIZED 

GENERALIZED 


MODE 

COORDINATE 

FREQUENCY 

POTENTIAL 

0 

0=FLEXIBLE 

GENERALIZED 

STIFFNESS 

DAMPING VALUE 


NUMBER 

NUMBER 

RPM 

ENERGY 

FACTOR 1 

*RIGID BODY WEIGHT-LB 


LB/IN 

(LB-SEO/IN 


1 

14 

125. 

1.000E 02 

0. 

1 

4.503E 02 


0. 

0. 


2 

15 

159. 

6.400E 01 

0. 

1 

1.784E 02 


0. 

0. 


3 

16 

30637 . 

2.380E 06 

15. 

0 

1.787E 02 


4.760E 06 

9.891E 01 


NUMBER 

OF SUBSYSTEM 

MODES* 3 









THE MODE SHAPES AND 

CORRESPONDING FORCES 

FOR THIS SUBSYSTEM ARE- 





(SIGNS 

IN THETA-Z DIRECTION CHANGED TO OBTAIN RIGHT 

HAND COORDINATE SYSTEM) 








MODAL DISPLACEMENTS 

MODAL FORCES 





LOCAL 

GENERALIZED 


GLOBAL DIRECTION 

GLOBAL DIRECTION 





MODE 

COORDINATE 

POINT 

3 

4 

3 

4 





NUMBER 

NUMBER 

NUMBER 

Y 

THETA-Z 

Y 

THETA 

-Z 




1 

14 

1 

1 .ooooo 

0. 

0. 

0. 





1 

14 

2 

1 . ooooo 

0. 

0. 

0. 





1 

14 

37 

1.00000 

0. 

0. 

0. 





1 

14 

3 

1.00000 

0. 

0. 

0. 





1 

14 

33 

1.00000 

0. 

0. 

0. 





2 

15 

1 

1 .ooooo 

-0.02000 

0. 

0. 





2 

15 

2 

-0.80000 

-0.02000 

0. 

0. 





2 

15 

37 

0. 

-0.02000 

0. 

0. 





2 

15 

3 

0.80000 

-0.02000 

0. 

0. 





2 

15 

33 

1 .ooooo 

-0.02000 


0. 





3 

16 

1 

1 .00000 

0.03122 

0. 

0. 






16 

2 

0.61170 

0.03028 

-9.672E 

05 -1.123E 

07 





16 

37 

-0.71850 

0. 

4.312E 

05 -6.803E 

07 





16 

3 

0.61170 

-0.03028 

9.672E 

05 -1.123E 

07 





16 

33 

1.00000 

-0.03122 

0. 

0. 







TOTAL NUMBER OF SUBSYSTEMS= 4 


TOTAL NUMB E R OF MODES OR GENERALIZED COORDINATES^ 16 

SUMMARY OF THE MODES OR GENERALIZED COORDINATES- 

GENERALIZED 

COORDINATE GENERALIZED GENERALIZED GENERALIZED 
NUMBER WEIGHT STIFFNESS DAMPING VALUE 


I 

r 

8 



1 .372E 

03 

0. 

0. 

2 

1 .449E 

03 

0. 

0. 

3 

3.743E 

02 

1.076E 06 

0. 

4 

2.457E 

03 

1.980E 07 

0. 

5 

1 . 766E 

02 

3.500E 06 

0. 

6 

1 . 372E 

03 

0. 

0. 

7 

1 .449E 

03 

0. 

0. 

8 

3.743E 

02 

1.076E 06 

0. 

9 

2.457E 

03 

1.980E 07 

0. 

10 

1 .766E 

02 

3.500E 06 

6 . 

1 1 

4.503E 

02 

0. 

0. 

12 

1 .784E 

02 

0. 

0. 

13 

1.787E 

02 

4.760E 06 

9.891E 01 

14 

4.503E 

02 

0. 

0. 

15 

1.784E 

02 

0. 

0. 

16 

1 .787E 

02 

4.760E 06 

9.891E 01 

SUMMARY OF 

THE COORDINATES FOR THE PHYSICAL 

POINTS- 

POINT 

COORDINATES ( INCHES ) -GLOBAL SYSTEM 

NUMBER 


X 

Y 

z 

1 


0. 

0. 

0. 

2 


- 10.000 

0. 

0. 

3 


-90.000 

0. 

0. 

4 


0. 

0. 

0. 

5 


-10.000 

0. 

0. 

6 


-90.000 

0. 

0. 

33 


- 100.000 

o. 

o. 

34 



0. 

0. 

35 


-70.000 

0. 

0. 

36 


- 100.000 

0. 

0. 

37 



0. 

0. 

38 


-50.000 

0. 

0. 

39 


-10.000 

0. 

0. 

40 


-90.000 

0. 

0. 








PHYSICAL CONNECTING ELEMENT NUMBER 1 


ELEMENT TYPE= 5 

NUMBER OF END POINTS* 2 

POINT NUMBER AT I END* 2 

POINT NUMBER AT J END* 39 

NUMBER OF DIRECTIONS FOR POINT AT I END* 5 
GLOBAL DIRECTIONS FOR POINT AT I END* 1 2 

NUMBER OF DIRECTIONS FOR POINT AT J END* 5 
GLOBAL DIRECTIONS FOR POINT AT J END* 1 2 

SPRING CONSTANT IN X DIRECTION* 0^ 

SPRING constant IN Y DIRECTION* 5.000E 05 
SPRING CONSTANT IN Z DIRECTION* 5.000E 05 
SPRING CONSTANT IN THETA-Y DIRECTION* 0. 
SPRING CONSTANT IN THETA-Z DIRECTION* 0. 


O-FACTOR* 


FREQUENCY* 


60.4 HERTZ 


DAMPING CONSTANTS (CALCULATED BASED ON ABOVE O-FACTOR AND FREQUENCY) - 


DAMPING CONSTANT IN X DIRECTION* O. 

DAMPING CONSTANT IN Y DIRECTION* 8.783E 01 
DAMPING CONSTANT IN Z DIRECTION* 8 . 783E 01 
DAMPING CONSTANT IN THETA-Y DIRECTION* 0. 
DAMPING CONSTANT IN THETA-Z DIRECTION* 0. 
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PHYSICAL CONNECTING ELEMENT NUMBER 


ELEMENT TYPE= 


NUMBER OF END POINTS=' 2 
POINT NUMBER AT I END= 3 


POINT NUMBER AT 0 END= 40 

NUMBER OF DIRECTIONS FOR POINT AT I END= 5 
GLOBAL DIRECTIONS FOR POINT AT I END= 1 2 


NUMBER OF DIRECTIONS FOR POINT AT J END*= 5 
GLOBAL DIRECTIONS FOR POINT AT J END= 1 2 


SPRING CONSTANT IN X DIRECTION= 

SPRING CONSTANT IN Y DIRECTION* 5.000E 05 
SPRING CONSTANT IN 2 DIRECTION* 5.000E 05 
SPRING CONSTANT IN THETA-Y DIRECTION* 


SPRING CONSTANT IN THETA-Z DIRECTION 


O-FACTOR* 15.0 


FREQUENCY*- 


60.4 HERTZ 


DAMPING CONSTANTS (CALCULATED BASED ON ABOVE O-FACTOR AND FREQUENCY )- 


DAMPING CONSTANT IN X DIRECTION* 


DAMPING CONSTANT IN Y DIRECTION 
DAMPING CONSTANT IN Z DIRECTION* 


DAMPING CONSTANT IN THETA-Y DIRECTION* 
DAMPING CONSTANT IN THETA-Z DIRECTION 


8.783E 01 
8.783E 01 
ON* O. 
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( 


PHYSICAL CONNECTING ELEMENT NUMBER 3 


( 


ELEMENT TYPE= 5 

NUMBER OF END POINTS= 2 

POINT NUMBER AT I END= 5 
POINT NUMBER AT J END* 2 

NUMBER OF DIRECTIONS FOR POINT AT I END* 5 

GLOBAL DIRECTIONS FOR POINT AT 1 END* 12345 

NUMBER OF DIRECTIONS FOR POINT AT J END* 5 

GLOBAL DIRECTIONS FOR POINT AT J END* 12345 

SPRING CONSTANT IN X DIRECTION* 

SPRING CONSTANT IN Y DIRECTION* 1.000E 06 

SPRING CONSTANT IN Z DIRECTION* 1.000E 06 

SPRING CONSTANT IN THETA-Y DIRECTION* 


SPRING CONSTANT IN THETA-Z DIRECTION* O. 


O-FACTOR* 


FREQUENCY* 


60.4 HERTZ 


DAMPING CONSTANTS (CALCULATED BASED ON ABOVE Q-FACTOR AND FREQUENCY )- 


DAMPING CONSTANT IN X DIRECTION* 


DAMPING CONSTANT IN Y DIRECTION 
DAMPING CONSTANT IN Z DIRECTION* 


DAMPING CONSTANT IN THETA-Y DIRECTION* 
DAMPING CONSTANT IN THETA-Z DIRECTION* 


0 . 


1.757E 02 
1,757E 02 
ON* 0. 


S 
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PHYSICAL CONNECTING ELEMENT NUMBER 4 


ELEMENT TYPE= 5 


NUMBER OF END POINTS' 2 

POINT NUMBER AT I END= 6 ^ 

POINT NUMBER AT 0 END= 3 

NUMBER OF DIRECTIONS FOR POINT AT I END= 5 

GLOBAL DIRECTIONS FOR POINT AT I END' 12345 


NUMBER OF DIRECTIONS FOR POINT AT J END' 5 

GLOBAL DIRECTIONS FOR POINT AT J END' 12345 


SPRING constant IN X DIRECTION' 0- 
SPRING CONSTANT IN Y DIRECTION' 1 .OOOE 06 
SPRING CONSTANT IN Z DIRECTION' 1.000E 06 

SPRING CONSTANT IN THETA-Y DIRECTION' 0. 

SPRING CONSTANT IN THETA-2 DIRECTION' 0. 

O-FACTOR' 15.0 FREQUENCY' 60.4 HERTZ 


DAMPING CONSTANTS (CALCULATED BASED ON ABOVE O-FACTOR AND FREOUENCY)- 

DAMPING CONSTANT IN X DIRECTION' 0. 

DAMPING CONSTANT IN Y DIRECTION' 1.757E 02 
DAMPING CONSTANT IN Z DIRECTION' 1.757E 02 
DAMPING CONSTANT IN THETA-Y DIRECTION' 0. 

DAMPING CONSTANT IN THETA-Z DIRECTION' 0. 



( 


PHYSICAL CONNECTING ELEMENT NUMBER 5 


ELEMENT TYPE= 3 

NUMBER OF END POINTS= 2 

POINT NUMBER AT I END= 4 

POINT NUMBER AT 0 END= 1 

NUMBER OF DIRECTIONS FOR POINT AT I END= 2 
GLOBAL DIRECTIONS FOR POINT AT I ENO= 1 3 


NUMBER OF DIRECTIONS FOR POINT AT J END= 2 
GLOBAL DIRECTIONS FOR POINT AT d END= 1 3 

LOCAL RADIAL SPRING RATE= 1 .OOE 06 LB/ IN 
RADIAL DEAD BAND= 10.0 MILS 
LOCAL DAMPING COEFFICIENT^ O. 


(L8-SEC)/IN 




o 


I 

c 

I 

H 

I 


\ 

/ 


NUMBER OF TYPE 3 PHYSICAL CONNECTING ELEMENTS* 1 
NUMBER OF TYPE 4 PHYSICAL CONNECTING ELEMENTS* 0 


NUMBER OF TYPE 5 PHYSICAL CONNECTING ELEMENTS* 4 
TOTAL NUMBER OF PHYSICAL CONNECTING ELEMENTS* 5 


THIS RUN IS NOT A RESTART RUN. 

TIME STEP* 0.000050 SECONDS . 

FINAL TIME* 0.512000 SECONDS 
PRINT MULTIPLE* 1000 

PLOT MULTIPLE* 10 

INDEPENDENT ROTOR NUMBER (ONE FOR WHICH SPEED-TIME HISTORY IS INPUT)* 1 
BEGINNING TIME FOR FIRST SEGMENT* 0. SECONDS 

BEGINNING SPEED FOR FIRST SEGMENT- 3000. RPM 


ENDING ACCEL 

TIME RATE 

SEGMENT SECONDS RPM/SEC 


1 0.600000 O. 


TOTAL NUMBER OF SPEED SEGMENTS FOR INDEPENDENT ROTOR SPEED-TIME 
HISTORY* 1 

SUMMARY OF UNBALANCE LOAD INPUT - 


BIRTH POINT MAGNITUDE PHASE ANGLE 
TIME (SEC.) NUMBER GM-IN DEGREES 

0. 4 5000. O. 


TOTAL NUMBER OF UNBALANCE BIRTH EVENTS* 1 
TOTAL NUMBER OF P*COS(WT) AND P*SIN(WT) LOADS* 0 
TOTAL NUMBER OF TIME -FORCE HISTORY LOADS* O 
SUMMARY OF THE GYROSCOPIC LOAD INPUT- 


POLAR MOMENT 

POINT OF INERTIA 

NUMBER LB-IN**2 


4 184205. 

36 184205. 


TOTAL NUMBER OF GYRO LOAD LOCATIONS* 2 


THIS RUN PRODUCES A PLOT FILE (FILE CODE 23). 


TIMES. ROTOR SPEEDS, AND ROTOR ANGULAR DISPLACEMENTS 
ARE WRITTEN ONTO THE PLOT FILE. 


} 

NUMBER OF TYPE 1 PHYSICAL CONNECTING ELEMENTS* 0 

NUMBER OF TYPE 2 PHYSICAL CONNECTING ELEMENTS* 0 







( 


( 


( 


0 > 


DISPLACEMENTS, VELOCITIES. MODAL FORCES AND COORDINATES ARE WRITTEN 
ONTO THE PLOT FILE FOR THE FOLLOWING POINTS AND DIRECTIONS- 

GLOBAL 

POINT DIRECTION 

NUMBER NUMBER DIRECTION 


1 1 Z 

1 3 Y 

4 1 Z 

4 3 Y 

37 1 Z 

37 3 Y 

38 1 Z 

38 3 Y 

TOTAL NUMBER OF POINTS AND DIRECTIONS FOR DISPLACEMENT. VELOCITY 
MODAL FORCE. AND COORDINATE PLOT FILE OUTPUT= 8 

THE RELATIVE DISPLACEMENT MAGNITUDE. CLEARANCE, AND FORCE MAGNITUDE 
IS WRITTEN TO THE PLOT FILE FOR ALL TYPE 3 PHYSICAL CONNECTING 
ELEMENTS (RUB ELEMENTS)(IF ANY). 


PHYSICAL CONNECTING ELEMENT FORCES ARE WRITTEN ONTO THE PLOT FILE FOR THE 
FOLLOWING PHYSICAL CONNECTING ELEMENTS. POINTS, AND DIRECTIONS- 


ELEMENT 

POINT 

DIRECTION 


NUMBER 

NUMBER 

NUMBER 

DIRECTION 

3 

2 

1 

2 

3 

2 

3 

Y 

3 

5 

1 

2 

3 

5 

3 

Y 

4 

3 

1 

Z 

4 

3 

3 

Y 

4 

6 

1 

Z 

4 

6 

3 

Y 


TOTAL NUMBER OF ELEMENTS, POINTS, AND DIRECTIONS FOR ELEMENT FORCE PLOT 
FILE OUTPUT* 8 

SUMMARY OF THE CONNECTIONS BETWEEN THE PHYSICAL CONNECTING ELEMENTS 
AND THE MODAL SUBSYSTEMS- 

ELEMENT 1 IS CONNECTED TO SUBSYSTEM 7 AT POINT 2 

ELEMENT 1 IS CONNECTED TO SUBSYSTEM 8 AT POINT 2 

ELEMENT 2 IS CONNECTED TO SUBSYSTEM 7 AT POINT 3 

ELEMENT 2 IS CONNECTED TO SUBSYSTEM 8 AT POINT 3 

ELEMENT 3 IS CONNECTED TO SUBSYSTEM 1 AT POINT 5 

ELEMENT 3 1$ CONNECTED TO SUBSYSTEM 2 AT POINT 5 

ELEMENT 3 IS CONNECTED TO SUBSYSTEM 7 AT POINT 2 

ELEMENT 3 IS CONNECTED TO SUBSYSTEM 8 AT POINT 2 

ELEMENT 4 IS CONNECTED TO SUBSYSTEM 1 AT POINT 6 

ELEMENT 4 IS CONNECTED TO SUBSYSTEM 2 AT POINT 6 

ELEMENT 4 IS CONNECTED TO SUBSYSTEM 7 AT POINT 3 

I ELEMENT 4 IS CONNECTED TO SUBSYSTEM 8 AT POINT 3 

« ELEMENT 5 IS CONNECTED TO SUBSYSTEM 1 AT POINT 4 

I ELEMENT 5 IS CONNECTED TO SUBSYSTEM 2 AT POINT 4 

g ELEMENT 5 IS CONNECTED TO SUBSYSTEM 7 AT POINT 1 

I ELEMENT 5 IS CONNECTED TO SUBSYSTEM 8 AT POINT 1 
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) 


POINT 

4 

OF 

SUBSYSTEM 

2 

IS 

A 

GYRO 

LOAD 

LOCATION 

POINT 

36 

OF 

SUBSYSTEM 

1 

IS 

A 

GYRO 

LOAD 

LOCATION 

POINT 

36 

OF 

SUBSYSTEM 

2 

IS 

A 

GYRO 

LOAD 

LOCATION 


SUMMARY OF THE GYRO LOAD LOCATIONS ON THE MODAL SU8SYSTEMS- 
POINT 4 OF SUBSYSTEM 1 IS A GYRO LOAD LOCATION 


( 

SPEED 

r 

SEGMENT NUMBER 

TIME* 

= 1 

0. SECONDS 

( 


( 

ROTOR 

PROPERTIES 

FOR 

INDEPENDENT ROTDR (ROTOR 1)- 




SPEED* 

3000. 

RPM 







ACCELERATION* 


0. RPM/SEC 






ANGULAR DISPLACEMENT 

* 0. 

REVOLUTIONS 







DISPLACEMENTS 

IN GIVEN DIRECTION 



POINT 

X 


V 

Z 

THETA-X 

theta-y 

THETA-Z 


NUMBER 

INCHES 


INCHES 

INCHES 

RADIANS 

RADIANS 

RADIANS 


1 

0. 


O. 

0. 

0. 

0. 

0. 


2 

0. 


0. 

0. 

0. 

0. 

O. 


3 

0. 


0. 

O. 

O. 

0. 

0. 


4 

0. 


0. 

0. 

0. 

0. 

0. 


5 

0. 


0. 

0. 

o. 

0. 

0. 


6 

0. 


0. 

0. 

0. 

0. 

0. 


33 

0. 


0. 

0. 

0. 

0, 

0. 


34 

0. 


0- 

0. 

0. 

0. 

0. 


35 

0. 


0. 

0, 

0. 

O. 

0. 


36 

0. 


o. 

0. 

0. 

0. 

0. 


37 

0. 


0. 

0. 

0. 

0. 

0. 


38 

0. 


0. 

0. 

0. 

0. 

0. 


39 

0. 


0. 

0. 

0. 

0. 

0. 


40 

0. 


0. 

0. 

0. 

0. 

0. 






VELOCITIES 

IN GIVEN DIRECTION 



POINT 

X 


Y 

z 

THETA-X 

theta-y 

THETA-Z 


NUMBER 

IN/SEC 


IN/SEC 

IN/SEC 

RAD/SEC 

RAO/SEC 

RAO/SEC 


1 

0. 


o. 

0. 

0. 

0. 

0. 


2 

0. 


o. 

0. 

o. 

O- 

0. 


3 

0. 


0. 

0. 

o. 

0. 

O. 


4 

0. 


0. 

0- 

0. 

0. 

0. 


5 

0. 


0. 

0. 

o. 

0. 

0. 


6 

0. 


0. 

0. 

0. 

0. 

0. 


33 

o. 


0. 

o. 

o. 

o. 

O. 


34 

0. 


0. 

0. 

0. 

0. 

0. 


35 

0. 


0. 

0. 

0. 

0. 

0. 


36 

0. 


0. 

0. 

0. 

0. 

0. 


37 

0, 


0. 

0. 

0. 

0. 

0. 


38 

0. 


0. 

0. 

0. 

0. 

0. 


39 

0. 


O- 

0. 

0. 

0. 

0. 


40 

0. 


0. 

0. 

0. 

0. 

0. 





FORCES 

CONTRIBUTED 

BY THE SUBSYSTEM 

MODE SHAPES 



POINT 

X 


Y 

z 

THETA-X 

theta-y 

THETA-Z 


NUMBER 

POUNDS 


POUNDS 

POUNDS 

IN-LB 

IN-LB 

IN-LB 


1 

0. 


0. 

0. 

0. 

0. 

0. 


2 

0. 


0. 

0. 

0. 

0. 

0. 


3 

o. 


O. 

0. 

O. 

0. 

O. 


4 

0. 


0. 

0. 

0- 

0. 

0. 


5 

0. 


0. 

0. 

0. 

0. 

0. 


6 

0. 


0. 

0. 

0. 

0. 

0. 


33 

0. 


0. 

0. 

0. 

0. 

0. 


34 

0. 


0. 

0. 

0. 

0. 

0. 


35 

0. 


0. 

imimsii 

0. 

0. 

0. 



i: 
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; 



36 

0. 

O. 

37 

0. 

0. 

36 

0. 

0. 


I 

o. o. o. 

O. O. 0, 

g. 


o. 


o. 

o. 


) 


THE FOLLOWING VALUES ARE FOR THE TYPE 3 PHYSICAL CONNECTING ELEMENTS. 
FORCES ARE THOSE THAT THE ELEMENT EXERTS ON THE ROTOR OR CASE- 


RELATIVE DEAD I -END JEND FORCE IN Y DIRECTION FORCE IN Z DIRECTION "ibRCE' 

ELEMENT DISPLACEMENT BAND CLEARANCE POINT POINT I END J END I END J ENO MAGNITUDE 

NUMBER INCHES INCHES INCHES NUMBER NUMBER POUNDS POUNDS POUNDS POUNDS POUNDS 


5 “ O. ' 0.0100 0.0100 4 1 0. 0. O. 0. "" 67 


THE FORCES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 
S PRING-DAMPER ELEMENTS) EXERT ON THE ENGINE COMPONENTS OR GROUND ARE- 






FORCE IN 

GIVEN DIRECTION 



ELEMENT 


POINT 

X 

Y 

Z 

THETA-Y 

THETA-Z 

NUMBER 

END 

NUMBER 

POUNDS 

POUNDS 

POUNDS 

IN-LB 

IN-LB 

1 

I 

2 

0. 

0. 

0. 

0. 

0. 

1 

J 

39 

0. 

0. 

0. 

0. 

0. 

2 

I 

3 

0. 

0. 

0. 

0. 

0. 

2 

J 

40 

0. 

0. 

0. 

0. 

0. 

3 

I 

5 

0. 

0. 

0. 

0. 

O. 

3 

J 

2 

0. 

0. 

0. 

0. 

0. 

4 

I 

6 

0. 

0. 

0. 

0. 

0. 

4 

J 

3 

0. 

0. 

0. 

0. 

0. 

THE GYROSCOPIC 

FORCES ACTING ON 

THE ROTOR(S) ARE- 







POLAR MOMENT 

Y-AXIS 

Z-AXIS 




POINT 

ROTOR 

OF INERTIA 

MOMENT 

MOMENT 




NUMBER 

NUMBER 

LB-1N**2 

IN-LB 

IN-LB 




4 

1 

184205. 

0. 

0. 




36 

1 

184205. 

0. 

0. 





SUMMARY OF UNBALANCE FQRCES- 


BIRTH PHASE 

TIME POINT ROTOR MAGNITUDE ANGLE FORCE (LB.) FORCE (LB.) 

SECONDS NUMBER NUMBER GM-IN DEGREES Y-PIRECTION Z-PIRECTION 


0. 41 5000. O. 2815.522 O. 

GENERALIZED 

GENERALIZED FORCE DUE 

COORDINATE TO APPLIED 

NUMBER FORCES ONLY 


1 O. 

2 0 . 

3 O. 

4 0. 


5 

6 

7 

8 

6. 

-2815.522 

2815,522 

343.916 




9 

773.705 




10 

-291.716 




1 1 

0. 








( 


( 


( 


12 O. 

13 O. 

14 0 . 


15 

O. 









16 

0. 









GENERALIZED 




GENERALIZED 

GENERALIZED 

GENERALIZED 



COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

WEIGHT 

STIFFNESS 

DAMPING VALUE 

GENERALIZED 


NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

POUNDS 

LB/IN 

(LB-SEO/IN 

ACCELERATION 


1 

0. 

0. 

0. 

1 .372E 

03 

0. 

0. 

0. 


2 

0. 

0. 

0. 

1 .449E 

03 

0. 

0- 

0. 


3 

0. 

0. 

0. 

3.743E 

02 

1.076E 06 

0. 

0. 


4 

0. 

0. 

0. 

2.457E 

03 

1.980E 07 

0. 

0. 


5 

0. 

0. 

0. 

1 .766E 

02 

3.500E 06 

0. 

0. 


6 

0. 

0. 

“2815.522 

1 .372E 

03 

0. 

0. 

-792.8347 


7 

0. 

o. 

2815.522 

1 .449E 

03 

0. 

0. 

750.9274 


8 

0. 

O. 

343.916 

3.743E 

02 

1.076E 06 

0. 

355.0783 


9 

0. 

0. 

773.705 

2.457E 

03 

1.9B0E 07 

0. 

121.6797 


10 

0. 

0. 

-291.716 

1.766E 

02 

3.500E 06 

0. 

-638.2326 


11 

0. 

0. 

0. 

4 . 503E 

02 

0. 

0. 

0. 


12 

0- 

0- 

0. 

1 .784E 

02 

0. 

0- 

0. 


13 

O. 

0. 

0. 

1.787E 

02 

4.760E 06 

9.891E 01 

0. 


14 

0. 

0. 

0. 

4.503E 

02 

0. 

0. 

0. 


15 

0. 

0. 

0. 

1.784E 

02 

0. 

0. 

O. 


16 

0. 

o. 

0. 

1.787E 

02 

4.760E 06 

9.891E Oi 

O. 




ROTOR PROPERTIES FOR INDEPENDENT ROTOR (ROTOR 1)- 
SPEED= 3000. RPM 


ACCELERATION= 0. RPM/SEC 

ANGULAR DISPLACEMENT® 2.50000218 REVOLUTIONS 


POINT 

NUMBER 


X 

INCHES 


DISPLACEMENTS IN GIVEN DIRECTION 
Y Z THETA-X THETA-Y 

INCHES INCHES RADIANS RADIANS 



0 

-0.03277097 

0.00831543 

0 

-0-00023358 

0 

-0.02860382 

0.00596802 

0 

-0.00023343 

0 

0.00473958 

-0.01233554 

0 

-0.00022416 

0 

-0.04511452 

0.01 107266 

0 

-0.00031388 

0 

-0.03852580 

0.00805975 

0 

-0.00031162 

0 

0.00700484 

-0.01765176 

0 

-0.00035026 

0 

0.00890828 

-0.01456402 

0 

-0.00022401 

0 

-0.02587625 

0.CX)203800 

0 

-0.00031369 

0 

-0.00327497 

-0.01072466 

0 

-0.00033705 

0 

0.01207768 

-0.02122501 

O 

-0.00035569 

0 

-0.01 193478 

-0.00338748 

0 

-0.00022879 

0 

-0.01413032 

-0.00416926 

0 

-0.00032337 


POINT 

NUMBER 


X 

IN/SEC 


IN/SEC 


VELOCITIES IN GIVEN DIRECTION 


THETA-X THETA-Y 

IN/SEC RAO/SEC RAD/SEC 


0. 145031 



K 7.907604 

4.907780 0 

0.198268 

>. 1.729063 

-3.755285 0 

0.145098 

>. 2.053834 

-4.363334 0 

0.200898 

0. 

0. 0 

0. 


POINT 


NUMBER 


X 


POUNDS 


FORCES CONTRIBUTED BY THE SUBSYSTEM MODE SHAPES 
Y Z THETA-X THETA-Y 


POUNDS POUNDS IN-LB IN-LB 


0 . 

- 1.935 


O. 

’148. 1 19 


3 

0 

1.935 

148. 119 

0 

4 

O 

0. 

0. 

0 

5 

0 

420.714 

458.600 

0 

6 

0 

-952.426 

- 156.844 

0 

33 

0 

0. 

0- 

0 

34 

O 

1473.846 

1 196.969 

0 

35 

0 

-860.611 

-65.747 

3 


0 . 

1719.704 


1719.704 

0 . 

2368. 106 
12167.978 



THETA-Z 

RADIANS 

^0,00041673 
-0.00041673 
-0.0004 1685 
-0.00065406 
-0.00064084 
-0.00053106 
-0.00041685 
-0.00060806 
-0.00054421 
-0.00053139 
-0.00041679 
-0.00057262 
0 . 

O. 


THETA-Z 

RAD/SEC 

-0.080491 
-0.080416 
-0.075585 
-O. 1 13788 
-0. 1 10744 
-O.’l 15872 
-0.075510 
-0. 108675 
-O. 1 14362 
-0. 1 15298 
-0.078000 
-0. 1 10877 


IHETA-Z 

IN-LB 

0 . 

-22.461 

-22.461 

0 . 

29773.079 

4159.59 6 

0 . 

36917.246 

22445.777 









( 

36 

37 

38 

0. 

0. 

O. 

0. 

0.863 

-490.301 

0. 

66,039 

151.269 


0. 

0. 

0, 

( 

0. 

10418.075 

-12872.401 

0. 

- 136,073 
34733.307 

« 

4 

( 

THE FOLLOWING VALUES 

ARE FOR THE TYPE 

3 PHYSICAL CONNECTING ELEMENTS. 





l| FORCES ARE THOSE THAT THE ELEMENT EXERTS ON THE ROTOR OR 

CASE- 






RELATIVE 

DEAD 

I -END 

JEND 


FORCE IN Y DIRECTION 

FORCE IN Z 

DIRECTION 

FORCE 

ELEMENT 

DISPLACEMENT 

BAND CLEARANCE 

POINT 

POINT 


I END 

J END 

I END 

J END 

MAGNITUDE 

NUMBER 

INCHES 

INCHES INCHES 

NUMBER 

NUMBER 


POUNDS 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

5 

0.0126 

0.0100 -0.0026 

4 

1 


2584.069 

-2584.069 

-577.213 

577.213 

2647.751 


THE FORGES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 
SPRING-DAMPER ELEMENTS) EXERT ON THE ENGINE COMPONENTS OR GROUND ARE- 


FORCE IN GIVEN DIRECTION 


ELEMENT 


POINT 

X 

Y 


Z 

THETA-Y 

THETA-Z 



NUMBER 

END 

NUMBER 

POUNDS 

POUNDS 


POUNDS 

IN-LB 

IN-LB 



1 

I 

2 

0. 

14433.403 


-2144.638 

0. 

0. 



1 

J 

39 

0. 

- 14433,403 


2144.638 

0. 

0. 



2 

I 

3 

0. 

-2786.382 


5987.578 

0. 

0. 



2 

J 

40 

0. 

2786.382 


-5987.578 

0. 

0. 



3 

I 

5 

0. 

10030.745 


-1511.401 

0. 

0. 



3 

J 

2 

0. 

-10030.745 


1511.401 

0. 

0. 



4 

1 

6 

0. 

-2609.586 


5138.481 

0. 

0. 



4 

J 

3 

0. 

2609.586 


-5138.481 

0. 

0. 



THE GYROSCOPIC 

FORCES ACTING ON THE 

ROTOR(S) ARE- 









POLAR MOMENT 


Y-AXIS 

Z-AXIS 






POINT 

ROTOR 

OF INERTIA 


MOMENT 

MOMENT 






NUMBER 

NUMBER 

LB-IN»»2 


IN- LB 

IN-LB 







4 1 184205 . 1704 1.671 33690.387 

36 1 184205 . 17267 . 804 29693.79 1 


SUMMARY OF UNBALANCE FORCES- 


BIRTH 

TIME POINT 

SECONDS NUMBER 

ROTOR MAGNITUDE 

NUMBER GM-IN 

PHASE 

ANGLE 

DEGREES 

FORCE (LB.) 
Y-DIRECTION 

FORCE (LB. ) 
Z-DIRECTION 


0. 4 

1 5000. 

0. 

-2815.522 

-0.001 



GENERALIZED 





GENERALIZED 

FORCE DUE 





COORDINATE 

TO APPLIED 





NUMBER 

FORCES ONLY 





1 

0.001 





2 

-0.001 





3 

-0.000 





4 

-0.000 





5 

0.000 





6 

2815.522 





7 

-2815.522 





8 

-343.916 





9 

-773.705 






291 .715 





i1 

0 . 










) 




13 

1 £ 

15 

16 


O. 

0 . 

0 ^ 

0 . 

0 , 


GENERALIZED 




GENERALIZED 

GENERALIZED 

GENERALIZED 


COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

WEIGHT 

STIFFNESS 

damping' value” 

GENERALIZED " 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

POUNDS 

LB/IN 

(LB-SEO/IN 

ACCELERATION 

1 

-0.01617755 

10.053262 

6583.309 

1.372E 03 

0. 

0. 

1853.8220 

2 

-0.00497969 

-5.039456 

3049.866 

1.449E 03 

0. 

0. 

813.4293 

3 

-0.00067527 

-0.486346 

-979.663 

3.743E 02 

1.076E 06 

0. 

-261.2904 

4 

-0.00010449 

0.098095 

-1961.822 

2.457E 03 

1.980E 07 

0. 

16.8276 

5 

0.00013516 

0. 189776 

411.839 

1.766E 02 

3.500E 06 

0. 

- 133.9766 

6 

0.02B64669 

5.546710 

-11148.495 

1.372E 03 

0. 

0. 

-3139.3522 

7 

-0.01626238 

2.338521 

7189.706 

1.449E 03 

0. 

0. 

1917.5658 

8 

-0.00189543 

0. 113431 

-1765.899 

3.743E 02 

1.076E 06 

0. 

282.4671 

9 

0.00018411 

0.033574 

3560.983 

2.457E 03 

1.980E 07 

0. 

- 13.2750 

10 

0.00023662 

-0. 171256 

809.847 

1.766E 02 

3.500E 06 

0. 

-40.0935 

11 

-0.00327744 

-3.753732 

793.073 

4.503E 02 

0. 

0. 

680.5414 

12 

-0.01143973 

7.254919 

608.654 

1.784E 02 

0. 

0. 

1318.2907 

13 

0.00015315 

0.002161 

709.255 

1.787E 02 

4.760E 06 

9.891E 01 

-43. 1302 


14 

-0.01 193335 

1.671748 

1641.794 

4.503E 02 

0. 

0. 

1408.8335 

15 

0.02083962 

3 . 900020 

-1079.494 

1.784E 02 

0. 

0. 

-2338.0892 

16 

0.00000200 

-0.079770 

0.891 

1.787E 02 

4.760E 06 

9.891E 01 

- 1 . 6009 






66 66 

0.02746222 

0.02343173 

-0.00872697 

0.03889252 

-0.00475721 

-0.00470091 

-0.00399713 

-0.00653721 

5 

6 

0. 

o. 

0.03248879 

-0.01268604 

-0.00626488 

-0.00566416 

33 

0. 

-0.01273616 

-0.00387748 

34 

o. 

0.02015537 

-0.00575915 

35 

0. 

-0.00224095 

-0.00542387 

36 

0. 

-0.01790092 

-0.00584569 

37 

0. 

0.00731590 

-0.00445754 

38 

0. 

0.00859837 

-0.00544643 



ROTOR PROPERTIES FOR INDEPENDENT ROTOR (ROTOR 1)- 
SPEED= 3000. RPM 


ACCELERATION^ 0. RPM/5EC 

ANGULAR DISPLACEMENT = 5,00000525 REVOLUTIONS 


POINT 

NUMBER 


POINT 

NUMBER 


POINT 




X 

INCHES 


IN/SEC 



DISPLACEMENTS IN GIVEN DIRECTION 
V Z THETA'X THETA-y 

INCHES INCHES RADIANS RADIANS 


0.00000625 
0.00000633 
0.00001 127 
0.00001866 


0.00001916 
-0.00001 130 
0.00001134 
0.CKX)01542 


-0.00000339 

-0.00001369 

0.00CX)0880 

0.00000677 



IN/SEC 
2.001 101 


1.971567 

1.627004 

2.672810 

2.560899 


2.317475 

1.570397. 

2.362936 

2.238117 


VELOCITIES IN GIVEN DIRECTION 


THETA-X 
IN/SEC RAD/SEC 


THETA-V 

RAD/SEC 

-0. 133852 


O. 133831 
-O. 132466 
-0.214916 
-0.210981 


-0. 176755 
-O. 132445 
-0.200439 
-O. 179594 


O. 178139 
-0. 133149 
-O. 188636 


FORCES CONTRIBUTED BY THE SUBSYSTEM MODE SHAPES 


8.928284 


7,588047 

-3.063852 

12-745767 

10.567732 


-4.413607 
-4.386590 
6.382647 
- 1 .043631 


; 0. 

2.373811 

-6.090387 ( 

0. 

1.845658 

2.232126 ( 

1 0. 

f 0. 

2.246493 

O. 

2.512109 ( 

O. ( 


THETA-X 


THETA-Y 



THETA-Z 

RADIANS 

0.00040284 
0.00040281 
0.000401 15 
0.00063848 
0.00062634 
0.(XX)54221 
0.000401 13 
0. 00059776 
0"00054793 
0.00054586 
0.00040198 
0.(XX)S6675 
0 . 

0 . 


THETA-Z 

RAD/SEC 

0.003219 

O.CK>3251 

0.005363 

0.008544 

0.008345 

-0.002897 

0.005395 

0.006393 

-0.000572 

-0.003352 

0.(X)4307 

0.002977 


THETA-Z 

IN-LB 

0 . 

-307.91 7 

-307.917 

0 . 

-26832.642 

3005.351 

O. 

-31278. 178 
-14505.715 
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) ) 


36 

0. 

0. 

0. 0. 

0. 

o. 



37 

0. 

1 1.824 

35.177 0. 

5549.427 

- 1865.379 



38 

0. 

534.765 

27.539 0. 

-10438.742 

-27269.947 



THE FOLLOWING VALUES 

ARE FOR THE TYPE 

3 PHYSICAL CONNECTING ELEMENTS. 




FORCES 

ARE THOSE THAT 

THE ELEMENT EXERTS ON THE ROTOR OR 

CASE- 





RELATIVE 

DEAD 

I -END JEND 

FORCE IN Y 

DIRECTION 

FORCE IN Z DIRECTION 

ri)RCE 

ELEMENT 

DISPLACEMENT 

BAND CLEARANCE 

POINT POINT 

I END 

J END 

I END J END 

MAGNITUDE 

NUMBER 

INCHES 

INCHES INCHES 

NUMBER NUMBER 

POUNDS 

POUNDS 

POUNDS POUNDS 

POUNDS 

5 

0.0116 

0.0100 -0.0016 

4 1 

-1549.398 

1549.398 

241.282 -241.282 

1ogF'^7'2 


THE FORCES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 
SPRIN G-DAMPER ELEMENTS) EXERT ON THE ENGINE COMPONENTS OR GROUND ARE- 


FORCE IN GIVEN DIRECTION 

ELEMENT POINT X Y Z THETA-Y THETA-2 

NUMBER END NUMBER PO UNDS POUNDS POUNDS IN-LB LN-LB_ 


112 O. -11889.037 1683.969 O 

1 J 39 0. 11889.037 -1683.969 O 

2 I 3 4220.579 2 267.676 O 

0 . 
O. 
O. 
O. 

4 J 3 O. -3037.776 -1904.136 O. O, 

THE GYROSCOPIC FORCES ACTING ON THE ROTOR(S) ARE- 


POLAR MOMENT Y-AXIS Z-AXI5 
POINT ROTOR OF INERTIA MOMENT MOMENT 
NUMBER NUMBER LB-IN««2 IN-LB IN-LB 


4 1 184205. -1279.643 -32187.216 

36 1 184205. 502.057 -26679.177 

SUMMARY OF UNBALANCE FORCES- 


BIRTH PHASE 

TIME POINT ROTOR MAGNITUDE ANGLE FORCE (LB.) FORCE (LB.) 

SECONDS NUMBER ^MBER GM-IN DEGREES Y-DIRECTION Z-DIRECTION 


O. 4 1 5000. O. 2815.522 0.019 

GENERALIZED 

GENERALIZED FORCE DUE 

COORDINATE TO APPLIED 

NUMBER FORCES ONLY 


2 

J 

40 

0. 

-4220.579 

-2267.676 

0. 

3 

I 

5 

0. 

-9160.585 

1040.527 

0. 

3 

J 

2 

0. 

9160.585 

-1040.527 

0. 

4 

I 

6 

0. 

3837.776 

1904. 136 

0. 


0 . 

0 . 

0 . 


1 

t 


3 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 


-0.019 

0.019 

0.002 

0.005 

- 0.002 

2815.522 

2815.522 

343.916 

773.705 

-291.716 

O. 




( ■ ’ ( 


12 

0. 




13 

0. 




14 

0. 




15 

0. 




‘ 16 

0. 




GENERALIZED 




GENERALIZED 

COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

WEIGHT 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

POUNDS 

1 

0.00033499 

“9.437936 

434.035 

1.372E 03 

2 

, “0.00611250 

3.241049 

3185.963 

1.449E 03 

3 

“0.00053598 

0.599544 

-653.211 

3.743E 02 

4 

-0.00003920 

“0.072269 

“730.449 

2.457E 03 

5 

0.00013009 

-0. 129396 

290.460 

1.766E 02 

6 

-0.02845566 

-0. 148541 

10309 . 892 

1.372E 03 

7 

0.01029312 

2.493996 

“4056.685 

1.449E 03 

8 

0.00146076 

0. 198019 

1309.367 

3.743E 02 

9 

-0.00021448 

0.003487 

-4057.978 

2.457E 03 

10 

-0.00023399 

“0.049484 

“743.870 

1.766E 02 

1 1 

-0.00439892 

2.248315 

765.700 

4.503E 02 

12 

0.00043986 

-6.657437 

17.361 

1.784E 02 

13 

0.00008158 

0.022532 

374.689 

1.787E 02 

14 

0.00733561 

1 .820610 

-796.251 

4.503E 02 

15 

-0.02009919 

-0.215352 

939.607 

1.784E 02 

16 

0.00002742 

-0.034861 

114.564 

1.787E 02 


GENERALIZED 

GENERALIZED 


STIFFNESS 

DAMPING VALUE 

GENERALIZED 

. LB/IN 

(LB-SEO/IN 

ACCELERATION 

0. 


0. 

122.2217 

0. 


0. 

849.7279 

1 .076E 

06 

0- 

-78.9803 

1 .980E 

07 

0. 

7. 1772 

3.500E 

06 

0. 

“360.7131 

O. 


0. 

2903.2063 

0. 


0. 

-1081.9581 

1 .076E 

06 

0- 

-270.9261 

1 .980E 

07 

0. 

29.6704 

3.500E 

06 

0. 

164.3333 

0. 


0. 

657.0524 

0. 


0. 

37.6020 

4.760E 

06 

9.891E 01 

-34.2763 

0. 


0. 

-683.2680 

0. 


0. 

2035. 1062 

4.760E 

06 

9.891E 01 

-27.0606 




0. 

-1.522249 

-6.562371 

0. 

0.087676 

0. 

-1.236602 

-5.687009 

0. 

Q. 087693 

0. 

1.052557 

1.372801 

0. 

0.088802 

0. 

-2 . 140803 

-9.686391 

0. 

0. 152856 

0. 

-1.764106 

-8. 186424 

0. 

0. 147781 

0. 

1.555216 

2.049304 

0. 

0. 1 19804 

0. 

1.339200 

2.262391 

0. 

0.088820 

1 0. 

- 1.008783 

-5.353252 

b. 

0. 137825 

> 0. 

0.637592 

-0.292410 

0. 

0. 122937 

> 0. 

2.030940 

3.205218 

0. 

0. 1 19439 

' 0. 

-0.093730 

-2. 181473 

0. 

0.088248 

\ 0. 

\ 0. 

-0.217100 

0. 

-2.736721 

0. 

0. 

0. 

0. 129159 
0. 


POINT 


NUMBER 


X 


POUNDS 


FORCES CONTRIBUTED BY THE SUBSYSTEM MODE SHAPES 
Y Z THETA-X THETA-Y 


POUNDS POUNDS IN-LB IN- LB 


0 . 

“62.729 


0 . 

3.664 



0 . 

-42.538 


-42.538 

0 . 

1930.647 

-1761.441 


0 . 

1274.031 

-778.288 


THETA-Z 

RADIANS 

-6.00028367 

-0.00028369 

-0.00028761 

-0.00048880 

-0.00047714 

-0.00038470 

-0.00028767 

- 0.00044832 

-0.00039385 

-0.00038662 

-0.00028665 

-0.00041756 

0 . 

0 . 


THETA-Z 

RAD/SEC 

-0. 028574 

-0.028576 

-0.028653 

-0.039846 

■0.0 39515 

-0.046193 

-0.028655 

-0.040005 

-0.044113 

-0.046959 

-0.028614 

-0.041774 

0 . 

O. 


THETA-Z 

IN-LB 

0 . 

-728.307 

-728.307 

0 . 

26251 . 171 
575.666 
O. 

32342.404 

17864.532 



( 

36 

37 

38 

0. 

0. 

0. 

• t 

0. 

27.968 

•482.042 

O. 0. 
-1.634 0. 
54.395 0. 

( 

O. 

-257.698 

308.477 

0. 

-4412. 128 
29767.813 

( 


/ 

THE FOLLOWING VALUES 
FORCES ARE THOSE THAT 

ARE FOR THE TYPE 3 PHYSICAL CONNECTING ELEMENTS. 
THE ELEMENT EXERTS ON THE ROTOR OR CASE- 






RELATIVE 

DEAD 

I -END UEND 

FORCE IN Y 

DIRECTION 

FORCE IN Z 

DIRECTION 

FORCE 

ELEMENT 

DISPLACEMENT 

BAND CLEARANCE 

POINT POINT 

I END 

J END 

I END 

U END 

magnitude 

NUMBER 

INCHES 

INCHES INCHES 

NUMBER NUMBER 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

5 

0.0104 

0.0100 -0.0004 

4 1 

351.603 

-351.603 


68.836 

358.278 


THE FORCES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 
SPRING-DAMPER ELEMENTS) EXERT ON THE ENGINE COMPONENTS OR GROUND ARE- 


FORCE IN GIVEN DIRECTION 


ELEMENT 

NUMBER 

END 

POINT 

NUMBER 

X Y 

POUNDS POUNDS 

Z 

POUNDS 

THETA-Y 

IN-LB 

THETA-Z 

IN-LB 


1 

I 

2 

0. 9284.354 

- 1390.968 

0. 

0. 


1 

J 

39 

0. -9284.354 

1390.968 

0. 

0. 


2 

I 

3 

0. -2342.675 

1209.079 

0. 

0. 


2 

J 

40 

0. 2342.675 

-1209.079 

0. 

0. 


3 

I 

5 

0, 8159.802 

-1167.863 

0. 

O. 


3 

J 

2 

0. -8159.802 

1167.863 

0. 

0. 


4 

I 

6 

0. -2294.221 

1 1 13.588 

0. 

0. 


4 

J 

3 

0. 2294.221 

-1113.588 

0. 

0. 


1 THE GYROSCOPIC FORCES ACTING ON THE 

ROTOR(S) ARE- 







POLAR MOMENT 

Y-AXIS 

Z-AXIS 




POINT 

ROTOR 

OF INERTIA 

MOMENT 

MOMENT 




NUMBER 

NUMBER 

LB-IN^*2 

IN-LB 

IN-LB 




4 

1 

184205. 

5967.537 22892.702 




36 

1 

184205. 

7032.891 17887.896 




SUMMARY 

OF UNBALANCE FORCES- 






BIRTH 



PHASE 





TIME 

POINT 

ROTOR MAGNITUDE 

ANGLE FORCE ( LB . ) 

FORCE (LB. ) 




SECONDS 

NUMBER 

NUMBER GM-IN 

DEGREES Y -DIRECT ION 

Z-OIRECTION 




0. 

4 

1 5000. 

0. -2815.522 

-0,037 






GENERALIZED 






II GENERALIZED 

FORCE DUE 






COORDINATE 

TO APPLIED 






NUMBER 


FORCES ONLY 






1 


0.037 






2 


-0.037 






3 


-0.004 






4 


-0.010 






5 


0.004 






6 


2815.522 






7 


-2815.522 






8 


-343.916 






9 


-773.705 






10 


291.716 






» I 


0. 
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12 

13 

14 

odd 



15 

0. 



16 

0 . 



GENERALIZED 




COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

1 

-0.00582716 

6.463391 

2154.041 

2 

0.00075243 

-3. 185382 

-123. 148 

3 

0.000015S6 

-0.470272 

42.252 

4 

-0.00002654 

0.063997 

-505.402 

5 

-0.00000437 

-0.021611 

21.704 

6 

0.02089362 

2.089969 

-6714.912 

7 

-0.01020442 

-0.072127 

3401.662 

8 

-0.00160249 

0. 117351 

-1610.352 

9 

0.00018456 

0.014912 

3571.934 

10 

0.00024142 

-0.027621 

910.685 

1 1 

0.00056314 

-2. 168310 

-58.778 

12 

-0.00402517 

4.412381 

186.041 

13 

-0.00000379 

0.018320 

-9.226 

14 

-0.00696519 

-0.092807 

724.495 

15 

0.01428245 

1.430725 

-586.601 

16 

0.00006486 

0.001283 

306.646 


|i 

r 

I 

If 

?! 


I 


) 


GENERALIZED 

GENERALIZED 

GENERALIZED 



WEIGHT 

STIFFNESS 

DAMPING VALUE 

GENERALIZED 


POUNDS 

LB/IN 

(LB-SEO/IN 

ACCELERATION 


1.372E 03 

0. 

0. 

606.5656 


1.449E 03 

6 . 

6 . 

-32.8447 


3.743E 02 

1.076E 06 

0. 

26.3343 


2.457E 03 

1.9B0E 07 

0. 

3. 1560 


1.766E 02 

3.500E 06 

0. 

80.9803 


1.372E 03 

6 . 

0. 

- 1890.8806 


1.449E 03 

0. 

0. 

907.2569 


3.743E 02 

1.076E 06 

0. 

117.6279 


2.457E 03 

1.980E 07 

0. 

- 12.9405 


1.766E 02 

3.500E 06 

0. 

143.7962 


4.503E 02 

0. 

0. 

-50.4378 


1.784E 02 

0. 

0. 

402.9487 


1.787E 02 

4.760E 06 

9.891E 01 

15. 1248 


4.503E 02 

6 . 

0. 

621.6937 


1.784E 02 

O. 

O. 

- 1270.9598 


1.787E 02 

4.760E 06 

9.891E 01 

-4.7825 





( 


T I ME = 0 . 2000000 SECONDS 


( 


SPEED SEGMENT NUMBER- 1 


ROTOR PROPERTIES FOR INDEPENDENT ROTOR (ROTOR t)- 
SPEED= 3000 . RPM 


ACCELERATION^’ 0. RPM/SEC 

ANGULAR DISPLACEMENT* 10.00001156 REVOLUTIONS 


POINT 

NUMBER 


X 

INCHES 


DISPLACEMENTS IN GIVEN DIRECTION 
Y 2 THETA-X THETA-Y 

INCHES INCHES RADIANS RADIANS 



0. 

0.02728557 

-0.00532500 

0. 

0.02344330 

-0.00433106 

0. 

-0.00730348 

0.00354902 

0. 

0.03863449 

-0.00727169 

0. 

0.03247054 

-0.00591165 

0. 

-0.01070009 

0.00512456 

1 0. 

-0.01114790 

0.00452510 

0. 

0.02062662 

-0.00320969 

» 0. 

“0.00076177 

0.00228454 

; 0. 

-0.01566360 

0.00656417 

0. 

0.00807360 

-0.0<X)36040 

t 0. 

0.00957046 

-0.00049036 


0.00009922 

0.00009920 

0.00009780 

0.00013899 


0.00013797 

0.00014283 

0.00009778 

0.00013735 


0.00014059 

0.00014378 

0.00009850 

0.00013846 


O. 

o. 


POINT 

NUMBER 


IN/SEC 


VELOCITIES IN GIVEN DIRECTION 


Y 2 THETA-X THETA-Y 

IN/SEC IN/SEC RAD/SEC RAD/SEC 


1.739153 


8.663948 




1 . 393604 
- 1 . 326227 
2.364448 
1.690633 

7.437793 0 
-2.372012 O 
12.248238 0 
10.304116 .0 


-1.897122 

-3.483391 0 


- 1.660636 

-3.598308 0 


0.974188 

6.557749 O 


-0.907108 

-0.274314 0 


-2.391637 

-5.091035 0 


0.014608 

2.533133 0 


0.052904 

0. 

3.039241 0 

0. 0 


-0. 122617 


O. 122617 
’O. 122628 
“O. 194153 
-0. 190481 


-0. 166855 
-O. 122628 
-O. 182016 
O. 167940 


0. 168279 
0. 122623 
0. 173656 


FORCES CONTRIBUTED BY THE SUBSYSTEM MODE SHAPES 


THETA-X 


THETA -Y 




THETA-Z 

RADIANS 

0.00038425 

0.00038425 

0.00038442 

0.00061479 

0.00060234 

0.00051731 

0.00038442 

0.00057320 

0.00052288 

0.00052112 

0.00038433 

0,00054382 

O. 

o. 


THETA-Z 

RAD/SEC 


0.034446 


0.034432 

0.033564 

0.050795 

0.049246 


0.048028 

0.033550 

0.047624 

0.048144 


0.047518 

0.033998 

0.047614 





TH ETA -Z 
IN-LB 


0 . 

31 . 190 

31 . 190 

0 . 

-27456.883 

3290.225 

O. 

-31770.476 

-14491.698 
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) 


36 

0. 

0. 


0. 

0. 

0, 

0. 



37 

0. 

- 1 . 

198 

-9.924 

0. 

-1565.515 

188.953 



38 

0. 

548. 

401 

-41.802 

0. 

1837.753 

-27534.555 



THE FOLLOWING VALUES 

ARE FOR 

THE 

TYPE 3 PHYSICAL CONNECTING ELEMENTS. 




FORCES 

ARE THOSE THAT 

THE ELEMENI 

EXERTS ON THE ROTOR OR 

CASE- 





RELATIVE 

DEAD 


I -END 

JEND 

FORCE IN Y DIRECTION 

FORCE IN Z DIRECTION 

FORCE 

ELEMENT 

DISPLACEMENT 

BAND CLEARANCE POINT 

POINT 

I END 

J END 

I END J END 

MAGNI TUOE 

NUMBER 

INCHES 

INCHES 

INCHES NUMBER 

NUMBER 

POUNDS 

POUNDS 

POUNDS POUNDS 

POUNDS 

5 

0.01 15 

0.0100 

-0. 

0015 4 

1 

-1492.869 

1492.869 

256.073 -256.073 

1514.672 

THE FORCES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 




SPRING- 

DAMPER ELEMENTS) EXERT 

ON 

THE ENGINE COMPONENTS OR 

GROUND ARE- 









FORCE IN GIVEN DIRECTION 

» 


ELEMENT 

POINT 


X 

Y 

Z 

THETA-Y 

THETA-Z 


NUMBER 

END NUMBER 


POUNDS 

POUNDS 

POUNDS 

IN-LB 

IN-LB 


1 

1 

2 


0. 

-11844.054 

1512.238 

0. 

0. 


1 

J 

39 


0. 

M844.054 

-1512.238 

0. 

0. 


2 

I 

3 


0. 

3768.227 

-1566.169 

0. 

0. 


2 

J 

40 


0. 

-3768.227 

1566.169 

0. 

0. 


3 

I 

5 


0. 

-9114.552 

1077.071 

0. 

0. 


3 

J 

2 


0. 

9114.552 

-1077.071 

0. 

0. 


4 

I 

6 


0. 

3496.899 

-1380.307 

0. 

0. 


4 

J 

3 


0, 

-3496.899 

1380.307 

0. 

0. 


THE GYROSCOPIC FORCES 

ACTING 

ON THE ROTOR(S) 

ARE- 






POLAR MOMENT Y-AXIS 2-AXIS 


POINT ROTOR OF INERTIA MOMENT MOMENT 

NUMBER NUMBER LB-IN*«2 IN-LB IN-LB 


4 I 184205 . -7607 . 362 -29077 . 600 

36 1 184205. -7116.568 -25202.552 

S UMMARY OF U NBALANCE FORCES- 

I 

BIRTH PHASE 

TIME POINT ROTOR MAGNITUDE ANGLE FORCE (LB.) FORCE (LB.) 

SECONDS NUMBER NUMBER GM-IN DEGREES Y-OIRECTION 2-DlRECTION 


0. 4 1 SOOO, O. 2815.522 0.055 

GENERALIZ ED 

GENERALIZED FORCE DUE 

COORDINATE TO APPLIED 

NUMBER FORCES ONLY 


1 -0.055 

2 0.055 

3 0.007 

_J 0.015 

5 -0.006 

6 -2815.522 

7 2815.522 

_8 343.916 

9 773.705 

10 -291.716 

11 0 . 


I 

r 





( 

12 

0. 



( 



( 

13 

0, 







14 

0. 







15 

0. 







16 

0. 







GENERALIZED 




GENERALIZED 

GENERALIZED 

GENERALIZED 


COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALI ZED 

WEIGHT 

STIFFNESS 

DAMPING Value 

GENERALIZED 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

POUNDS 

LB/IN 

(LB-SEO/IN 

ACCELERATION 

1 

0.00692529 

“8.688663 

-2516.508 

1.372E 03 

0. 

0. 

“708.6342 

2 

“0.00036806 

3.521061 

-47. 109 

1.449E 03 

b. 

0. 

-12.5644 

3 

0.00008765 

0.410984 

124.969 

3.743E 02 

1.076E 06 

0. 

31.6487 

4 

0.00002679 

-0.069238 

516.864 

2.457E 03 

1.980E 07 

0. 

-2. 1454 

5 

“0.00003465 

-0.070835 

-46.509 

1.766E 02 

3.500E 06 

0. 

163.5546 

6 

“0.02720978 

-2.381923 

9852. Ill 

1.372E 03 

6. 

0. 

2774.2979 

7 

0.01128086 

-0.007515 

-4295.000 

1.449E 03 

0. 

0. 

- 1 145.5191 

a 

0.00147599 

0.009714 

1447.442 

3.743E 02 

1.076E 06 

0. 

-145.2921 

9 

-0.00022101 

-0.014120 

-4215.955 

2.457E 03 

1.980E 07 

0. 

25. 1742 

10 

-0.00023441 

0.070133 

-747.694 

1.766E 02 

3.500E 06 

0. 

159. 1817 

1 1 

-0.00037694 

2.533002 

-6.768 

4.503E 02 

O. 

o. 

-5.8078 

12 

0.00492505 

“6. 131128 

-240.751 

1.784E 02 

0. 

0. 

-521 .4445 

13 

“0.00002301 

-0.000182 

-103.573 

1.787E 02 

4.760E 06 

9.891E 01 

12.9511 

14 

0.00807161 

0.024914 

-965.306 

4.503E 02 

0. 

0. 

-828.3349 

15 

-0.01921673 

-1.699894 

907.795 

1.784E 02 

0. 

0. 

1966.2032 

16 

-0.00000278 

0.014344 

-10.796 

1.787E 02 

4.760E 06 

9.891E 01 

2. 1765 



184 


SPEED SE 

SMENT NUMBER® 

TIME= 0 

1 

.2500000 SECONDS 

y 





ROTOR PROPERTIES 

FOR 

INDEPENDENT 

ROTOR (ROTOR 

D- 






SPEED® 

3000. 

RPM 









ACCELERA 

TION» 


0. RPM/SEC 








ANGULAR 

DISPLACEMENT = 

12.50001454 REVOLUTIONS 









DISPLACEMENTS IN 

GIVEN DIRECTION 





POINT 

X 


Y 

2 

THETA-X 

THETA-Y 

THETA-Z 




NUMBER 

INCHES 


INCHES 

INCHES 

RADIANS 

RADIANS 

RADIANS 





QjHHIH 


-0.02557192 

0.00464358 

0. 

-0.00006024 

-0.00035631 





0. 


-0.02201033 

0.00404220 

0. 

-0.00006025 

-0.00035633 





0. 


0.00654277 

-0.000B1002 

0. 

-0.00006105 

-0.00035750 







-0.03661259 

0.00644784 

0. 

-0.00009014 

-0.00058006 







-0.03078181 

0.00553558 

0. 

-0.00008945 

-0.00056801 





0. 


0.00961066 

-0.00121169 

0. 

-0.00008085 

-0.00048042 




33 

0. 


0.0101 1946 

-0.00142170 

0. 

-0.00006107 

-0.00035752 




34 

0. 


-0.01959914 

0.00374459 

0. 

-0.00008727 

-0.00053897 




35 

0. 


0.00040890 

0.00037883 

0. 

-0.00008208 

-0.00048692 




36 

0. 


0.0141984 1 

-0.00199962 

0. 

-0.00008081 

-0.00048392 




37 

0. 


-0.00775964 

0.00163375 

0. 

-0.00006065 

-0.00035691 




36 

0. 


-0.00921793 

0.00202322 

0. 

-0.00008451 

-0.00050891 




39 

0. 


0. 

0. 

0. 

0. 

0. 




40 

0. 


O. 

0. 

0. 

0. 

0. 







VELOCITIES IN 

GIVEN DIRECTION 





POINT 

X 


Y 

2 

THETA-X 

THETA-Y 

THETA-Z 




NUMBER 

IN/SEC 


IN/SEC 

IN/SEC 

RAD/SEC 

RAD/SEC 

RAD/SEC 




1 

0. 


-1.569837 

-8.00964B 

0. 

0. 111385 

-0.019882 





0 . 


-1.371367 

-6.896129 

0. 

0.111389 

-0.019886 





o. 


0.230613 

2.025686 

o. 

0. 1 1 1656 

-0.020163 





0. 


-2. 160661 

-11.500158 

o. 

O. 182766 

-0-028637 





0. 


-1.853177 

-9.654720 

0. 

0. 178982 

-0.029028 







0.344183 

2.970336 

0. 

0. 148627 

-0.026855 





■H 


0.432637 

3. 142621 

0. 

0. 1 1 1660 

-0.020168 





bh 


-1.241126 

-6. 120173 

0. 

O. 169400 

“0.028751 




35 

0. 


-0. 147072 

0. 133914 

0. 

0. 151217 

-0.026578 




36 

0. 


0.593904 

4.380524 

0. 

0. 149671 

-0.027627 




37 

0. 


-0.576464 

-2.441073 

0. 

0. 1 1 1523 

-0.020025 




38 

0. 


-0.667301 

-2.859568 

0. 

0. 159056 

-0.027565 




39 

0. 


0. 

O. 

0. 

0. 

0. 




40 

0. 


0. 

0. 

0. 

0. 

0. 







FORCES CONTRIBUTED BY 

THE SUBSYSTEM 

MODE SHAPES 





POINT 

X 


Y 

Z 

THETA-X 

THETA-Y 

THETA-Z 




NUMBER 

POUNDS 


POUNDS 

POUNDS 

IN-LB 

IN-LB 

IN-LB 




1 

O, 


o. 

0. 

O. 

O. 

0. 




2 

0. 


-18.802 

12.840 

0. 

-149.073 

-218.292 




3 

0. 


18.602 

-12.840 

0. 

- 149.073 

-218.292 




4 

0. 


0. 

0. 

0. 

0. 

0. 




5 

0. 


326.207 

-54.787 

0. 

1720.214 

26878.249 




6 

O. 


-936.554 

80.294 

0. 

491 . 465 

-2523.240 




33 

0. 


0. 

0. 

0. 

O. 

0. 




34 

0. 


1507.079 

-200.567 

O. 

2704-374 

32127.619 




35 

0. 


-866.668 

69. 137 

O. 

1999.420 

15585.857 













186 


12 

0. 



13 

0. 



14 

0. 



15 

0. 



16 

0. 



GENERALIZED 




COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

1 

-0.00422616 

7.958366 

1355.665 

2 

0.00220272 

-3.485190 

-964.582 

3 

0.00015017 

-0.528935 

166.804 

4 

-0.00000883 

0.065590 

- 174.493 

5 

-0.00002934 

0.096687 

-57.098 

6 

0.02546358 

1.378709 

-8907. 140 

7 

-0.01099476 

-0.771710 

4124. 138 

6 

-0.00152237 

-0.045883 

-1487.218 

9 

0.00021134 

0.005192 

4059.547 

10 

0.00025447 

0.058526 

815.480 

1 1 

0.00162421 

-2.437912 

-223.793 

12 

-0.00303264 

5.576134 

126.915 

13 

-0.00001328 

0.004399 

-62.551 

14 

-0.00774567 

-0.573176 

894.318 

15 

0.01784569 

1 .001237 

-806.902 

16 

0.00001944 

0.004576 

92. 130 




I 

c 


I 




k 


} 


] 


GENERALIZED GENERALIZED GENERA L I ZED 

WEIGHT^ STIFFNESS DAMPING^ VALUE' GENERALIZED 
POUNDS LB/IN (LB-SEO/IN ACCELERATION 


1.372E 

03 

0. 


0. 

381 .7474 


1 449E 

03 

0. 


0. 

-257.2635 


3.743E 

02 

1 .076E 

06 

0. 

5.3916 


2.457E 

03 

1.980E 

07 

0. 

0.0608 


1.766E 

02 

3 . 500E 

06 

0, 

99.7477 


1.372E 

03 

0. 


0. 

-2508. 1994 


1.449E 

03 

0. 


0. 

1099.9483 


3.743E 

02 

1 .076E 

06 

0. 

155.7435 


2.457E 

03 

1.980E 

07 

0. 

- 19.6533 


1.766E 

02 

3 . 500E 

06 

0. 

- 164.4296 


4.503E 

02 

0. 


0. 

- 192.0381 


1 .784E 

02 

0. 


0. 

274.8864 


1 .787E 

02 

4.760E 

06 

9.891E 01 

0.4475 


4.503E 

02 

0. 


0. 

767.4202 



1 .784E 02 
1.787E 02 


O. 

9.891E 01 


1747, 
- 1 


6792 

8524 


O. 

4.760E 06 



( ' ■ c 

TIME" 0.3000000 SECONDS 
SPEED SEGMENT NUMBER" 1 


ROTOR PROPERTIES FOR INDEPENDENT ROTOR (ROTOR 1)- 

SPEED» 3000. RPH 

ACCELERATION" 0. RPM/SEC 

ANGULAR DISPLACEMENT" 15.00001776 REVOLUTIONS 




DISPLACEMENTS IN 

GIVEN DIRECTION 



POINT 

X 

Y 

Z 

THETA-X 

THETA-Y 

THETA-Z 


NUMBER 

INCHES 

INCHES 

INCHES 

RADIANS 

RADIANS 

RADIANS 


1 

0. 

0.02486154 

-0.00534319 

0. 

0.00008801 

0.00034511 


2 

0. 

0.02141232 

-0.00446294 

0. 

0.00008801 

0.00034513 


3 

0. 

-0.00625826 

0.00257237 

0. 

0.00008788 

0.00034663 


4 

0. 

0.03573390 

-0.00744054 

0. 

0.00012886 

0.00056692 


5 

0. 

0.03004666 

-0.00616757 

0. 

0.00012758 

0.00055433 


6 

0. 

-0.00920327 

0.00374897 

0. 

0-00012573 

0.00046534 


33 

0. 

-0.00972668 

0.00345095 

0. 

0.00008787 

0.00034665 


34 

0. 

0.01916034 

-0.00365991 

0. 

0.00012556 

0.00052474 


35 

0. 

-0.00027704 

0.00127204 

0. 

0.00012460 

0.00047251 


36 

0. 

-0.01364711 

0.00499987 

0. 

0.00012662 

0.00046823 


37 

0. 

0.00760992 

-0.00094243 

0. 

0.00008794 

0.00034588 


38 

0. 

0.00907131 

-0.001 18867 

0. 

0.00012452 

0.00049461 


39 

0. 

0. 

0, 

0. 

0. 

0. 


40 

0. 

0. 

0. 

0. 

0. 

0. 




VELOCITIES IN 

GIVEN DIRECTION 



POINT 

X 

Y 

Z 

THETA-X 

THETA-Y 

THETA-Z 


NUMBER 

IN/SEC 

IN/SEC 

IN/SEC 

RAO/SEC 

RAD/SEC 

RAO/SEC 


1 

0. 

1.714370 

7.764889 

0. 

-0. 107436 

0.028594 


2 

0. 

1.428434 

6.691115 

0. 

-0. 107443 

0.028594 


3 

0. 

-0.859353 

- 1.922986 

0. 

-0. 107909 

0.028600 


4 

0. 

2.375689 

11. 143891 

0. 

-O. 176688 

0.040666 


5 

o. 

1.968059 

9.387116 

.o. 

-0. 172278 

0.040644 


6 

0. 

-1.263389 

-2.823882 

0. 

-0. 144715 

0.042179 


33 

0. 

-1. 145363 

-3.002737 

0. 

-0. 107917 

0.028600 


34 

0. 

1 . 156224 

6.030862 

0. 

-0. 162815 

0.040563 


35 

0. 

-0.449564 

-0.017487 

0. 

-0. 147340 

0,041072 


36 

0. 

-1.680513 

-4.217023 

0. 

-0. 145026 

0.042978 


37 

0. 

0.284667 

2.394302 

0. 

-0. 107676 

0-020597 


38 

0. 

0.351976 

2.905867 

0. 

-0. 153853 

0.040604 


39 

0. 

0. 

0. 

0. 

0. 

0. 


40 

0. 

0. 

0. 

0. 

0, 

0. 




FORCES CONTRIBUTED BY 

THE SUBSYSTEM 

MODE SHAPES 



POINT 

X 

Y 

Z 

THETA-X 

THETA-Y 

THETA-Z 


NUMBER 

POUNDS 

POUNDS 

POUNDS 

IN-LB 

IN-LB 

IN-LB 


1 

0. 

0. 

0. 

0. 

0. 

0, 


2 

0. 

23.911 

2.073 

O. 

'24.071 

277.610 


3 

0. 

-23.911 

-2.073 

0. 

-24.071 

277.610 


4 

0. 

0. 

0. 

0. 

0. 

0. 


5 

0. 

-289.025 

-39.589 

0. 

-2477.952 

-27794.224 


6 

0. 

896.396 

-45.704 

0. 

1630.753 

1463.907 


33 

0. 

0. 

0. 

0. 

0. 

0. 


34 

0. 

- 1275.207 

-71.258 

0. 

-1628.371 

-32368.677 


35 

0. 

832.686 

-53.020 

0. 

654.873 

- 15913.777 




188 


36 

0. 

0. 

0. 

0. 

/ 

0. 

0. 

37 

0. 

- 10.661 

-0.924 

0. 

- 145.825 

1681.780 

38 

0. 

525.952 

-59, 175 

0. 

-501 .461 

-28505.856 


THE FOLLOWING VALUES ARE FOR THE TYPE 3 PHYSICAL CONNECTING ELEMENTS. 
FORCES ARE THOSE THAT THE ELEMENT EXERTS ON THE ROTOR OR CASE- 


RELATIVE 

DEAD 


1-END 

JEND 

FORCE IN Y DIRECTION 

FORCE IN Z DIRECTION 

FORCE 

ELEMENT DISPLACEMENT 

BAND 

CLEARANCE 

POINT 

POINT 

I END 

J END 

I END 

J END 

MAGNITUDE 

NUMBER INCHES 

INCHES 

INCHES 

NUMBER 

NUMBER 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

5 0.0111' 

0.0100 

-0.0011 

4 

1 

-1053.3B0 

1053.380 

203.204 

-203.204 

1072.801 


THE FORCES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 
S PRING-DAMPER ELEMENTS) EXERT ON THE ENGINE COMPONENTS OR GROUND ARE- 


ELEMENT 

NUMBER 

END 

POINT 

NUMBER 

FORCE IN GIVEN DIRECTION 
X Y Z 

POUNDS POUNDS POUNDS 

THETA-Y 

IN-LB 

THETA-Z 

IN-LB 


1 

I 

2 

0. -10831.625 

1643.765 

0. 

0. 


1 

J 

39 

0. 10831.625 

-1643.765 

0. 

0. 


2 

I 

3 

0. 3204.609 

-1117.283 

0. 

0. 


2 

J 

40 

0. -3204.609 

1117.283 

0. 

0. 


3 

I 

5 

0. -8729. 129 

1231.021 

0. 

0. 


3 

J 

2 

0. 8729.129 

-1231.021 

0. 

0. 


4 

I 

6 

0. 3015.988 

- 1018.337 

0. 

0. 


4 

J 

3 

0. -3015.988 

1018.337 

0. 

0. 


THE GYROSCOPIC FORCES ACTING ON THE 

ROTOR(S) ARE- 







POLAR MOMENT 

Y-AXIS 

Z-AXIS 




POINT 

ROTOR 

OF INERTIA 

MOMENT 

MOMENT 




NUMBER 

NUMBER 

LB-IN»»2 

IN-LB 

IN-LB 




4 

1 

184205. 

-6090.348 -26461.887 




36 

1 

184205. 

-6436.726 -21719.919 




SUMMARY 

OF UNBALANCE FORCES- 






BIRTH 



PHASE 





TIME 

POINT 

ROTOR MAGNITUDE 

ANGLE FORCE (LB.) 

FORCE (LB. ) 




SECONDS 

NUMBER 

NUMBER GM-IN 

DEGREES Y-DIRECTION 

Z-DIRECTION 




0. 

4 

1 5000. 

0. 2815.522 

0.089 






GENERALIZED 







GENERALIZED FORCE DUE 

COORDINATE TO APPLIED 


NUMBER FORCES ONLY 


I 

t 


1 

2 

3 

4 

5 

6 
7 

9 
10 
1 1 


-0.089 
0.089 
0.01 1 
0.025 
-0.009 
-2815.522 
2815.522 
343.9 16 
773.705 
-291 .716 


0 . 









( 

12 

0. 

A 


13 

0. 



14 

0. 



15 

0. 



16 

0. 



GENERALIZED 




COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

1 

0.00622832 

-7.698082 

-2253.322 

2 

-0.00121685 

3.403180 

415.977 

3 

-0.00003087 

0.470751 

-40.950 

4 

0.00002953 

-0.064140 

566.725 

5 

-0.00000269 

-0.026562 

-7.096 

6 

-0.02474829 

-2.031166 

8597.588 

r 

0.01083237 

0.348386 

-3951.000 

8 

0.00154083 

-0.022973 

1495.205 

9 

-0.00021031 

-0.011151 

-4035.337 

10 

-0.00022022 

-0.019383 

-818.477 

1 1 

“0.00094397 

2.388772 

110.595 

12 

0.00439707 

-5.383813 

-206. 147 

13 

-0.00000214 

-0.007696 

-11.253 

14 

0-00759215 

0.284599 

-860.495 

15 

-0.01729411 

-1.429867 

779.514 

16 

-0.00002472 

-0-000095 

-117.337 



00 3 

<0 i 

I 




( 


GENERALIZED 

GENERALIZED 

GENERALIZED 



WEIGHT 

STIFFNESS 

DAMPING VAtUE 

GENERALIZED 


POUNDS 

LB/IN 

(LB-SEO/IN 

ACCELERATION 


1 .372E 

03 

0. 


0. 

-634.5224 


1 .449E 

03 

0. 


0. 

1 10.9452 


3.743E 

02 

1.076E 

06 

0. 

-7.9858 


2.457E 

03 

1.980E 

07 

0. 

-2.8370 


1.766E 

02 

3.500E 

06 

0. 

5.0933 


1.372E 

03 

0. 


0. 

2421.0314 


f .449E 

03 

O. 


o. 

- 1053. 7708 


3.743E 

02 

1 .076E 

06 

o. 

- 168.0105 


2.457E 

03 

1.980E 

07 

0. 

20.2520 


1.766E 

02 

3.500E 

06 

0. 

-104.3555 


4.503E 

02 

O. 


o. 

94.9023 


1 .784E 

02 

0. 


0. 

-446.4968 


1.787E 

02 

4 . 760E 

06 

9.891E 01 

-0.6230 


4.503E 

02 

0. 


0. 

-738.3959 


1.784E 

02 

0. 


0. 

1688.3584 


1.787E 

02 

4 . 760E 

06 

9.891E 01 

0.7615 






TIME= 0.3500000 SECONDS 

...SPIIP.. S.E -.MENT NUMBER^ 1 

POTOR PP';t5ERTIES FOR INDEPENDENT ROTOR (ROTOR 1)- 

SPEED«_ 3000. RP M 

ACCELEPt *rON= O. RPM/SEC 

ANGULAR irSPLACEMENTs 17.50002074 REVOLUTIONS 


DISPLACEMENTS IN GIVEN DIRECTION 


POINT 

NUMBER 


X 

INCHES 


1 0 . 

2 0 . 

3 0. 

4 0. 



Y 

Z 

THETA-X 

THETA-Y 

INCHES 

INCHES 

RADIANS 

RADIANS 

-0.02588359 

0.00504401 

0. 

-0.00008165 

-0.02225120 

0.00422724 

0. 

-0.00008164 

0.00685484 

-0.00229398 

0. 

-0.00008139 

-0.03696257 

0.00696776 

0. 

-0.0001 1831 

-0.03105714 

0.00580273 

0. 

-0.00011711 

0.01005826 

-0.00333814 

0. 

-0.0001 1598 

0.01049896 

-0.00310752 

0. 

-0.00008139 

-0.01972735 

0.00350565 

0. 

-0.0001 1540 

0.00063819 

-0.00104252 

0. 

-0.0001 1502 

0.01475840 

-0.00449619 

O. 

-0.0001 1661 

-0.00771828 

0.00096108 

0. 

-0.00008152 

-0.00918230 

0.00123256 

O. 

-0.00011473 


POINT 

NUMBER 


IN/SEC 


IN/SEC 


VELOCITIES IN GIVEN DIRECTION 


THETA-X THETA-Y 

IN/SEC RAD/SEC RAD/SEC 





-1.324966 
0.741222 
-2. 189696 
-1.817932 

-6.984489 
2. 145690 
-11.601 606 
-9.748956 


1 .075170 

3. 151010 


0.998440 

3.289017 


-1.094759 

-6. 189541 


0.344630 

0.200908 


1 .438521 

4.624574 


-0.288257 

-2.426440 


-0.379043 

0. 

-2.877580 

0. 


0. 1 13962 


0. 113967 
0. 1 14288 
O. 184148 
0. 180531 


0. 154176 
O. 1 14292 
0. 171704 
0. 155879 


0. 155476 
0. 1 14127 
0. 162515 
O. 


POINT 


NUMBER 


X 


POUNDS 


FORCES CONTRIBUTED BY THE SUBSYSTEM MODE SHAPES 


THETA-X 


THETA-Y 


) 


THETA-Z 

RADIANS 

-0.00036335 

-0.00036337 

-0.00036428 

-0.00058882 

-0.00057645 

-0.00049084 

-0.00036430 

-0.00054742 

-0.00049687 

-0.00049434 

-0.00036383 

-0.00051803 

0 . 

0 . 


THETA-Z 

RAD/SEC 

-0.025912 

-0.025910 

-0.025745 

-0.038727 

-0. 03 7852 

-0.035903 

-0.025742 

-0.036725 

-0.036197 

-0.035604 

-0.025827 

-0.036312 


T HETA-Z 

IN-LB 

O. 

- 169.669 

- 169.669 

0 . 

27300.832 

-2703.071 

0 . 

31703.978 
14838. 167 





( • ( 


36 

37 

38 

odd 

0. 

6.516 

-536.916 

0. 

1.797 

46.849 

666 

O. 

283.476 

149.768 

0. 

- 1027.868 
27642. 124 




1 THE FOLLOWING VALUES 

ARE FOR THE TYPE 

3 PHYSICAL CONNECTING ELEMENTS. 





j FORCES ARE THOSE THAT THE ELEMENT EXERTS ON THE ROTOR OR 

CASE- 






RELATIVE 

DEAD 

I -END 

JENO 

FORCE IN V DIRECTION 

FORCE IN Z 

DIRECTION 

rORCE 

ELEMENT 

DISPLACEMENT 

BAND CLEARANCE 

POINT 

POINT 

I END 

J END 

I END 

J END 

MAGNITUDE 

NUMBER 

INCHES 

INCHES INCHES 

NUMBER 

NUMBER 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

5 

0.0112 

0.0100 -0.0012 

4 

1 

1226.405 

-1226.405 

-212.953 

212.953 

1 2*44 ."756 


THE FORCES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 
SPRING-DAMPER ELEMENTS) EXERT ON THE ENGINE COMPONENTS OR GROUND ARE- 


ELEMENT 

NUMBER 

END 

POINT 

NUMBER 

X 

POUNDS 

FORCE IN GIVEN DIRECTION 
V Z THETA-y 

POUNDS POUNDS IN- LB 

THETA-Z 

IN-LB 


1 

I 

2 

O. 

1 124 1.978 

- 1500. 143 

0. 

0, 


1 

J 

39 

0. 

-11241.978 

1500. 143 

0. 

0. 


2 

I 

3 

0. 

-3492.524 

958.528 

0. 

0. 


2 

J 

40 

0. 

3492.524 

-958.528 

0. 

0. 


3 

I 

5 

O. 

8892. 535 

-1089.867 

0. 

0. 


3 

U 

2 

0. 

-8892.535 

1089.867 

0. 

0. 


4 

I 

6 

0. 

-3262.083 

867.553 

0. 

0. 


4 

J 

3 

0. 

3262.003 

-867. 5S3 

0. 

0. 


II THE GYROSCOPIC FORCES ACTING ON 

THE ROTOR(S) ARE- 







POLAR MOMENT 


Y-AXIS Z-AXIS 




POINT 

ROTOR 

OF INERTIA 


MOMENT MOMENT 




NUMBER 

NUMBER 

LB-IN^^2 


IN-LB IN- 

LB 





4 1 184205. 5799.931 27579.121 

36 1 184205. 5332.254 23285.045 


SUMMARY OF UNBALANCE FORCES- 


BIRTH PHASE 

TIME POINT ROTOR MAGNITUDE ANGLE FORCE (LB.) FORCE (LB.) 

SECONDS NUMBER NUMBER GMIN DEGREES Y-PIRECTION Z-DIRECTIQN 


0. 4 1 5000. 0. -2815.522 -0.107 


GENERALIZED 

GENERALIZED FORCE DUE 

COORDINATE TO APPLIED 

NUMBER FORCES ONLY 


CD 


I 

r 

s 


1 

2 

3 

5 

6 

7 

8 
9 

10 


O. 107 
-O. 107 
-0.013 
-0.029 
0.01 1 
2815.522 
-2815.522 
-343.916 
-773.705 
291.716 
O. 





I 


12 

0. 



13 

0. 



14 

0. 



15 

0. 



16 

0. 



GENERALIZED 




COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

1 

-0.00573884 

8. 131451 

2001.640 

2 

0.00123520 

-3.423439 

-435.375 

3 

0.00001689 

-0.459588 

13.963 

4 

-0.00002499 

0.066817 

-478.276 

5 

0.00001425 

0.086271 

-7.925 

6 

0.02591975 

1.816256 

-9151.860 

7 

-0.01089690 

-0.375441 

4041.335 

8 

-0.00148504 

-0.028073 

-1468.696 

9 

0.00021567 

0.007815 

4132.696 

10 

0.00022956 

-0.031675 

784.307 

11 

0.00096408 

-2.422637 

-106.348 

12 

-0.00407576 

5.706362 

188.048 

13 

0.00000417 

0.005293 

17.637 

14 

-0.00770742 

-0.290209 

892.597 

15 

0.01819128 

1 .291367 

-837.503 

16 

0.00001511 

-0.002718 

69.789 


c 


I 


S 




( 

SPEED 

•' * 

TIME= 0 

SEGMENT NUMBER^ 1 

.4000000 SECONDS 

( 


» * 


ROTOR 

PROPERTIES 

FOR INDEPENDENT 

ROTOR (ROTOR 

n- 





SPEED*^ 

3000. 

RPM 







ACCELERATION^ 

0. RPM/SEC 







' ANGULAR DISPLACEMENT:* 20.00002408 REVOLUTIONS 







DISPLACEMENTS IN 

GIVEN DIRECTION 




POINT 

X 

Y 

2 

THETA-X 

THETA-Y 

THETA-2 



NUMBER 

INCHES 

INCHES 

INCHES 

RADIANS 

RADIANS 

RADIANS 



1 

0. 

0.02541246 

-0.66503269 

0. 

0.00007793 

0.00035327 



2 

0. 

0.02188145 

-0.00425343 

0. 

0.00007793 

0.00035329 



3 

0. 

-0.00643449 

0.00198239 

0. 

0.00007796 

0.00035461 



4 

0. 

0.03640190 

-0.00699374 

0. 

0.0001 1467 

0.00057624 



5 

0. 

0.03061435 

-0.00585301 

0. 

0.0001 1359 

0.00056402 



6 

0. 

-0.00945753 

0.00289451 

0. 

0.00010953 

0.00047636 



33 

0. 

-0-00998246 

0.00276208 

0. 

0.00007796 

0.00035463 



34 

0. 

0.01952040 

-0.00360907 

0. 

0.00011148 

0.00053486 



35 

0. 

-0.00032740 

0.00074312 

0. 

0.00010905 

0.00048300 



36 

O. 

-0.01400781 

0.00397735 

0. 

0.0001 1023 

0.00047967 



37 

0. 

0.00775255 

-0.00113625 

0. 

0-00007795 

0.00035395 



38 

0. 

0.00922351 

-0.00141531 

0. 

0.00010980 

0.00050490 



39 

0. 

0. 

0. 

0. 

6. 

0. 



40 

0. 

0. 

0. 

o. 

O. 

0. 





VELOCITIES IN 

GIVEN DIRECTION 




POINT 

X 

Y 

2 

THETA-X 

THETA-Y 

THETA-Z 



NUMBER 

IN/SEC 

IN/SEC 

IN/SEC 

RAD/SEC 

RAO/SEC 

RAD/SEC 



1 

0. 

1.586089 

7.962738 

0. 

-0. 1 10448 

0.024349 



2 

0. 

1.342701 

6.858897 

0. 

-0, 1 10456 

0.024350 



3 

0. 

-0.608431 

-1.998010 

o. 

-0. 1 10966 

0.024428 



4 

0. 

2.202479 

1 1 .416348 

0. 

-0. 180144 

0.035118 



5 

0. 

1.844659 

9.595497 

o. 

-0. 176662 

0,035140 



6 

0. 

-0.888825 

-2.939413 

0. 

-0. 148619 

0.034614 



33 

0. 

-0.852826 

-3. 108395 

0. 

-O. 1 10974 

0.024430 



34 

0. 

1. 134515 

6. 100765 

0. 

-0. 167682 

0.034834 



35 

0- 

-0.225318 

-O. 108612 

o. 

-0. 150651 

0.034093 



36 

0. 

- 1 . 226124 

-4.351037 

0. 

-0. 149935 

0.035207 



37 

0. 

0.368861 

2.441652 

o. 

-0. 110711 

0.024389 



38 

0. 

0.443878 

2.868068 

0. 

-0. 157921 

0.034315 



39 

0. 

0. 

0. 

0. 

0. 

0. 



40 

0. 

0. 

0. 

0. 

0. 

0. 





FORCES CONTRIBUTED BY 

THE SUBSYSTEM 

MODE SHAPES 




POINT 

X 

Y 

Z 

THETA-X 

THETA-Y 

THETA-Z 



NUMBER 

POUNDS 

POUNDS 

POUNDS 

IN-LB 

IN-LB 

IN-LB 



1 

0. 

0. 

0. 

0. 

0. 

0. 



2 

0. 

21.137 

-0.530 

0. 

6. 152 

245.411 



3 

0. 

-21. 137 

0.530 

0. 

6. 152 

245.4 1 1 



4 

0. 

0. 

0. 

0. 

0. 

0. 



5 

0. 

-309.821 

-4.998 

0. 

-2237.058 

-27153.837 



6 

0. 

919.629 

-64.930 

0. 

1131. 500 

2221.680 



33 

0. 

0. 

0. 

0. 

0. 

0. 



34 

0, 

-1415.958 

46.394 

0. 

-2049.314 

-32103.534 



35 

o. 

852.587 

-65.591 

0. 

- 168.965 

- 15579.615 







194 


36 

37 
36 

o o o 

0. 

>9,424 

532.036 

0. 

0.236 

-53.642 

odd 

O. 

37.267 
- 1328.414 

0, 

1486.715 

-28379.419 




THE FOLLOWING VALUES 

ARE FOR THE TYPE 

3 PHYSICAL CONNECTING ELEMENTS. 





j FORCES ARE THOSE THAT THE ELEMENT EXERTS ON THE ROTOR OR 

CASE- 






RELATIVE 

DEAD 

I -END 

JEND 

FORCE IN Y DIRECTION 

FORCE IN 2 

DIRECTION 

FORCE 

ELEMENT 

DISPLACEMENT 

BAND CLEARANCE 

POINT 

POINT 

I END 

J END 

I END 

J END 

MAGNITUDE 

NUMBER 

INCHES 

INCHES INCHES 

NUMBER 

NUMBER 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

POUNDS 

5 

0.01 12 

0.0100 -0.0012 

4 

1 

-1144.962 

1 144.962 

204 .317 

-204,317 

i 163704 9 


THE FORCES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 
SPRING-DAMPER ELEMENTS) EXERT ON THE ENGINE COMPONENTS OR GROUND ARE- 


ELEMENT POINT 

NUMBER END NUMBER 

FORCE IN GIVEN DIRECTION 

X Y Z THETA-Y 

POUNDS POUNDS POUNDS IN-L8 

THETA-Z 

IN-LB 


1 I 

2 

0. -11058.660 

1524.269 

0. 

0. 


1 V 

39 

0. 11058.660 

-1524.269 

0. 

0. 


2 1 

3 

0. 3270.683 

-815.704 

0. 

0. 


2 V 

40 

0. -3270.683 

815.704 

0. 

0. 


3 1 

5 

0. -8821.072 

1118.851 

0. 

O. 


3 w 

2 

0. 8821.072 

-1118.851 

0. 

D. 


4 1 

6 

0. 3072.299 

-746.747 

0. 

0. 


4 c 

3 

0. -3072.299 

746.747 

0. 

0. 


THE GYROSCOPIC FORCES ACTING ON THE 

ROTOR(S) ARE- 






POLAR MOMENT 

Y-AXIS 

Z-AXIS 




POINT ROTOR OF INERTIA 

MOMENT 

MOMENT 




NUMBER NUMBER LB- IN* *2 

IN-LB 

IN-LB 




4 1 

184205. 

-5259.468 -26979.465 




36 1 

184205. 

-5272.765 -22455. 169 




SUMMARY OF UNBALANCE FORCES' 






BIRTH 


PHASE 





TIME POINT ROTOR MAGNITUDE 

ANGLE FORCE (LB.) 

FORCE (LB. ) 




SECONDS NUMBER NUMBER GM-IN 

DEGREES Y -DIRECT ION 

Z-DIRECTtON 




0. 

4 1 5000. 

0. 2815.522 

0. 125 





GENERALIZED 






GENERALIZED 

FORCE OUE 






COORDINATE 

TO APPLIED 






NUMBER 

FORCES ONLY 






1 

-0. 125 






2 

0. 125 






3 

0.015 






4 

0.034 






5 

-0.013 






6 

-2815.522 






7 

2815.522 






8 

343.916 






9 

773.705 






10 

-291 .716 






1 1 

0. 










12 

0. 




13 

0. 




14 

0. 




15 

0. 




16 

0. 




GENERALIZED 




GENERALIZED 

COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

WEIGHT 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

POUNDS 

1 

0.00549220 

-7.901605 

-1907.565 

1.372E 03 

2 

-0.00149829 

3.462017 

576.546 

1.449E 03 

3 

“0.00007063 

0.490791 

-69.564 

3.743E 02 

1 4 

0.00002421 

-0.065184 

464.576 

2.457E 03 

! 5 

0.00001235 

-0. 103159 

14.256 

1.766E 02 

6 

-0.02526288 

-1.716432 

8832.830 

1.372E 03 

7 

0.01098620 

0.483799 

-4078.214 

1.449E 03 

8 

0.00152150 

0.013084 

1500.575 

3.743E 02 

9 

-0.00021115 

-0.007754 

“4048.253 

2.457E 03 

1 10 

-0.00024119 

-0.026837 

-808.365 

1.766E 02 

1 1 

-0.00113585 

2.435598 

132. 144 

4.503E 02 

1 12 

0.00389739 

“5.535567 

-175. 183 

1.784E 02 

i 13 

0.00000055 

-0.008426 

1.497 

1.787E 02 

1 14 

0.00773685 

0.367929 

-894.241 

4.503E 02 

1 15 

-0.01769746 

“1.219457 

803.816 

1.784E 02 

1 

-0.00002186 

“0.001297 

- 102.418 

1.787E 02 


GE NERALIZED GENERALIZED 

STIFFNESS DAMPING VALUE GENERALIZED 
LB/IN (LB-SEO/IN ACCELERATION 


0. 


0. 

-537. 1593 


0. 


0. 

153.7706 


1 .076E 

06 

0. 

6.6436 


t .980E 

07 

0. 

-2.3351 


3.500E 

06 

0. 

-63.3680 


0. 


0. 

2487.2741 


0. 


0. 

- 1087.7001 


1.076E 

06 

0. 

-140.9936 


1 . 980E 

07 

o. 

20.8482 


3.500E 

06 

0. 

78.3281 


0. 


0. 

1 13.3940 


0. 


0. 

-379.4313 


4 . 760E 

06 

9.891E 01 

-0.6004 


0. 


0. 

“767.3538 


0. 


0. 

1740.9946 


4.760E 

06 

9.891E 01 

. 3.7633 





19(5 



! 


TIME® 0 

.4500000 SECONDS 

) 



) 

SPEED 

SEGMENT NUMBER 

= 1 







ROTOR 

PROPERTIES 

FOR 

INDEPENDENT 

ROTOR (ROTOR 

D* 





SPEED® 

3000. 

RPM 








II ACCELERATION® 


0. rpWsec 







ANGULAR DISPLACEMENT 

22.49989796 REVOLUTIONS 








DISPLACEMENTS IN 

GIVEN DIRECTION 




POINT 

X 


Y 

Z 

THETA-X 

THETA-Y 

THETA-Z 



NUMBER 

INCHES 


INCHES 

INCHES 

RADIANS 

RADIANS 

RADIANS 



1 

0. 


-0.02538712 

0.00523581 

0. 

-0.00008409 

-0.00035406 




0. 


-0.02184803 

0.00439483 

0. 

-0.(XX)08409 

-0.00035408 




0. 


0.00652519 

-0.00232857 

0. 

-0.00008400 

-0.00035526 




0. 


-0.03636473 

0.00726953 

0. 

-0.00012296 

-0.00057776 



5 

0. 


-0.03056947 

0.00605131 


-0.00012180 

-0.00056524 



6 

0. 


0.00958412 

-0.00339448 

0. 

-0.0001 1949 

-0.00047720 



33 

o. 


0.01007941 

-0.00316844 

0. 

-0.00008400 

-0.00035527 



34 

0. 


-0.01946490 

0.00365181 


-0.00011982 

-0.00053583 



35 

0. 


0.00042362 

-0.00104597 

0. 

-0.00011847 

-0.00048409 



36 

O. 


0.01414736 

-0.00458040 

o. 

-0.0001204 1 

-0.00048023 



37 

0. 


-0.00768733 

0.00103122 

0. 

-0.00008404 

-0.00035467 



38 

0. 


-0.00915290 

0.00129327 

0. 

-0.0001 1861 

-0.00050592 



39 

0. 


0. 

0. 

0. 

0. 

0. 



40 

0. 


0. 

0. 

0. 

0. 

0. 

o q 





VELOCITIES IN 

GIVEN DIRECTION 


- 


POINT 

X 


Y 

z 

THETA-X 

THETA-Y 

THETA-Z 



NUMBER 

IN/SEC 


IN/SEC 

IN/SEC 

RAD/SEC 

RAD/SEC 

RAO/SEC 

s| 


1 

0. 


- 1.641356 

-7.968073 

0. 

0.111035 

-0.026378 




0. 


- 1.377564 

-6.858186 

0. 

0. 111041 

-0.026378 

it; Q 



0. 


0.732146 

2.039868 

0. 

0. 111410 

-0.026365 



o. 


'2.280028 

-11. 405145 

0. 

O. 181249 

-0.038402 


5 

0. 


-1.897346 

-9.594643 

0. 

0. 177093 

-0.038138 


6 

0. 


1.069723 

2.994161 

0. 

0. 149686 

-0.037594 

< c/3 


33 

0. 


0.995781 

3. 154495 

.0. 

0. 1114 16 

-0.026365 


34 

0. 


-1. 141672 

-6. 128393 

0. 

0. 167766 

-0.037623 



35 

0. 


0.334277 

0. 107195 

0. 

0. 151977 

-0.037202 



36 

0. 


1.441705 

4.430858 

0. 

0. 150388 

-0.037957 



37 

0. 


-0.322439 

-2.417278 

0. 

0. 1 1 1226 

-0.026371 



38 

0. 


-0.399426 

-2.902432 

0. 

0. 158611 

-0.037250 



39 

0. 


0. 

0. 

0. 

O. 

0. 



40 

0. 


0. 

O. 

0. 

0. 

0. 






FORCES CONTRIBUTED BY 

THE SUBSYSTEM 

MODE SHAPES 




POINT 

X 


Y 

Z 

THETA-X 

THETA-Y 

THETA-Z 



NUMBER 

POUNDS 


POUNDS 

POUNDS 

IN-LB 

IN-LB 

IN-LB 



1 

0. 


0. 

0. 

0. 

0. 

0. 



2 

0. 


- 18.837 

-1.386 

0. 

16.092 

-218.707 



3 

0. 


18.837 

1 . 386 

O. 

16.092 

-218.707 



4 

0. 


0. 

0. 

0. 

0. 

0. 



5 

0. 


286.090 

25. 167 

0. 

2299. 156 

27642.973 



6 

0. 


-896.338 

56.233 

0. 

-1621 .987 

-1753.708 



33 

0. 


0. 

0. 

0. 

0. 

0. 



34 

0. 


1281.743 

9.046 

0. 

1720.860 

32157.419 



35 



-833.369 

60.894 


-461 .446 

15626.812 
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( 

36 

37 

38 

0. 

O. 

0. 

0. 

8.399 
-527 . 725 

O. 

0.618 

58.884 


0. 

0. 

0. 

( 

0. 

97.486 

740.263 

0. 

-1324.940 

28247.012 


( 

THE FOLLOWING VALUES 

ARE FOR THE TYPE 

3 PHYSICAL CONNECTING ELEMENTS. 




1 FORCES ARE THOSE THAT THE ELEMENT EXERTS ON THE ROTOR OR 

CASE- 





RELATIVE 

DEAD 

I -END 

JEND 


FORCE IN Y DIRECTION 

FORCE IN Z DIRECTION 

FORCE 

ELEMENT 

DISPLACEMENT 

BAND CLEARANCE 

POINT 

POINT 


I END 

J END 

I END vJ END 

MAGNITUDE 

NUMBER 

INCHES 

INCHES INCHES 

NUMBER 

NUMBER 


POUNDS 

POUNDS 

POUNDS POUNDS 

POUNDS 

5 

0.0112 

0.0100 -0.0012 

4 

1 


1144.928 

- 1 144.928 

-212.110 212.110 

i t 6 4T4”l0' 


THE FORCES THAT THE TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 
SPRING-DAMPER ELEMENTS) EXERT ON THE ENGINE COMPONENTS OR GROUND ARE- 


FORCE IN GIVEN DIRECTION 


ELEMENT 


POINT 


X Y 


Z 

THETA-Y 

THETA-Z 


NUMBER 

END 

NUMBER 


POUNDS POUNDS 


POUNDS 

IN-LB 

IN-LB 


1 

1 

2 


0. 11045.01 1 


1595.033 

0. 

0. 


1 

0 

39 


0. -11045.01 1 


1595.033 

0. 

0. 


2 

1 

3 


0. -3326.904 


985. 1 14 

0. 

0. 


2 

J 

40 


0. 3326.904 


-985. 114 

0. 

0. 


3 

I 

5 


0. 8812.753 

- 

1 175.773 

0. 

0. 


3 

J 

2 


0. -8812.753 


1175.773 

0. 

0. 


4 

I 

6 


0. -31 18.223 


898.274 

0. 

O. 


4 

u 

3 


0. 3118.223 


-898.274 

0. 

0. 


THE GYROSCOPIC FORCES ACTING ON THE 

ROTOR(S) ARE- 








POLAR 

MOMENT 

Y-AXIS 

Z-AXIS 





POINT 

ROTOR 

OF INERTIA 

MOMENT 

MOMENT 





NUMBER 

NUMBER 

LB- 

IN**2 

IN-LB 

IN-LB 





4 

1 


184205. 

5751.313 27144.989 





36 

1 


184205. 

5684.619 22523. 1 12 

• 




SUMMARY 

OF UNBALANCE FORCES- 







BIRTH 

) 



PHASE 






TIME 

POINT 

ROTOR 

MAGNITUDE 

ANGLE FORCE (LB.) 

FORCE 

(LB. ) 




SECONDS 

NUMBER 

NUMBER 

GM-IN 

DEGREES Y-DIRECTION 

Z-DIRECTION 





0. 4 1 5000. O. -2815.521 2.142 


GENERALIZED 

GENERALIZED FORCE DUE 

COORDINATE TO APPLIED 

NUMBER FORCES ONLY 


I 

t 

8 


1 

2 

3 

.A 

5 

6 

7 

8 
9 

10 
1 1 


-2. 142 


2. 142 
0.262 
0.589 
- 0.222 


2815.521 
-2815.521 
-34 3.916 
‘ -773.705 
291.716 
0 . 
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I 


12 

0 . 

13 

0 . 

14 

o. 

15 

0. 

16 

0 . 


GENE I ZED 


COORD lUATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

1 

-0.00593271 

7.936141 

2097.925 

2 

0.00133908 

-3.424323 

-487.467 

3 

0.00003997 

-0.471789 

42.197 

4 

-0.00002819 

0.066009 

-541.084 

5 

-0.00000584 

0.050102 

-6.095 

6 

0.02531411 

1.863206 

-8867.550 

7 

-0.01089946 

-0.417135 

4023.938 

8 

-0.00152487 

-0.010356 

-1483.481 

9 

0.00021130 

0.008513 

4054.201 

10 

0.00022168 

0.007352 

804.978 

1 1 

0.00103225 

-2.412892 

-120.310 

12 

-0.00420212 

5.561284 

192.770 

13 

0.00000143 

0.006103 

8.769 

14 

-0.00767333 

-0.322585 

878.649 

15 

0.01773326 

1.318569 

-807.823 

16 

0.00001948 

-0.000203 

92.894 


9 

r 

8 , 


f 


GENERALIZED 

GENERALIZED 

GENERALIZED 



WEIGHT 

STIFFNESS 

OAMPiNG VALUE 

GENERALlZeD 


POUNDS 

LB/IN 

(LB-5EO/IN 

ACCELERATION 


1.372E 

03 

0. 

0. 

590.7635 


1 .449E 

03 

0. 

0. 

- 130.0123 


3.743E 

02 

1.076E 06 

0. 

-0.8325 


2.457E 

03 

1.980E 07 

0. 

2.6712 


1.766E 

02 

3.500E 06 

0. 

31.3885 


1.372E 

03 

0. 

6 . 

-2497.0512 


1.449E 

03 

0. 

0. 

1073.2240 


3.743E 

02 

1.076E 06 

0. 

162.3856 


2.457E 

03 

1 . 980P 07 

0. 

-20.3865 


1 . 766E 

02 

3.500E 06 

0. 

63.6412 


4.503E 

02 

0. 

0. 

-103.2387 


1.784E 

02 

0. 

0. 

417.5234 


1.787E 

02 

4.760E 06 

9.891E 01 

2 . 9060 


4.503E 

02 

0. 

6. 

753.9741 


1 .784E 

02 

o. 

O. 

-1749.6741 


1.787E 

02 

4.760E 06 

9.891E 01 

0.4400 




( ■ ( 

TIME= 0.^000000 SECONDS 


SPEE D SEGMENT NUMB ER=» t 


ROTOR PROPERTIES 

FOR 

INDEPENDENT 

ROTOR (ROTOR 

D- 




II SPEED = 

3000. 

RPM 







ACCELERATION= 


0. RPM/SEC 






ANGULAR 

DISPLACEMENT= 

24.99974918 REVOLUTIONS 







DISPLACEMENTS IN 

GIVEN DIRECTION 



POINT 

X 


Y 

z 

THETA-X 

THETA-Y 

THETA-Z 


NUMBER 

INCHES 


INCHES 

INCHES 

RADIANS 

RADIANS 

RADIANS 


1 

0. 


0.02551177 

-0.00512239 

0. 

0.00008136 

0.00035611 


2 

0. 


0.02195211 

-0.00430869 

0. 

0.00008136 

0.00035613 


3 

0. 


-0.00658546 

0.00219568 

0. 

0.00008125 

0.00035731 


4 

0. 


0.03651526 

-0.00710485 

0. 

0.0001 1906 

0.00057969 


5 

0. 


0.03069712 

-0.00592832 

0. 

0.0001 1781 

0.00056740 


6 

0. 


-0.00967227 

0.00320016 

0. 

0.0001 1522 

0.00048053 


33 

0. 


-0.01016019 

0 . 00300808 

0. 

0.00008125 

0.00035733 


34 

0. 


0.01954141 

-0.00361285 

0. 

0-00011578 

0.00053832 


35 

0. 


-0.00045457 

0.00092791 

0. 

0.0001 1443 

0.00048703 


36 

0. 


-0.01426733 

0.00434566 

0. 

0.0001 1591 

0.00048383 


37 

0. 


0.00770913 

-0.00105427 

0. 

0.00008130 

0.00035672 


38 

0. 


0.00917577 

-0.00133505 

o. 

0.0001 1461 

0.00050863 


39 

0. 


0. 

0. 

0. 

6. 

0. 


40 

0. 


0. 

0. 

0. 

0. 

0. 





VELOCITIES IN 

GIVEN DIRECTION 



POINT 

X 


Y 

2 

THETA-X 

THETA-Y 

THETA-Z 


NUMBER 

IN/SEC 


IN/SEC 

IN/SEC 

RAD/SEC 

RAD/SEC 

RAD/SEC 


1 

0. 


1.599841 

8.012369 

0. 

-O. 111856 

0.025413 


2 

0. 


1 . 345630 

6.894241 

0. 

-0. 111861 

0.025412 


3 

0. 


-0.684586 

-2.068456 

0. 

-0. 1 12206 

0.025344 


4 

0. 


2.220163 

1 1 .470458 

0. 

-0. 182110 

0.037864 


5 

0. 


1.854763 

9.641096 

0. 

-0. 178289 

0.037109 


6 

0. 


-0.995855 

-3.038490 

“0. 

-0. 150824 

0.035598 


33 

0. 


-0.937930 

-3. 191003 

o. 

-O. 1 12211 

0.025343 


34 

0. 


1. 141820 

6. 133135 

0. 

-0. 169155 

0.036132 


35 

0. 


-0.276457 

-O. 147535 

o. 

-0. 152884 

0.035734 


36 

0. 


- 1 . 355066 

-4.479548 

0. 

-0. 151886 

0.035424 


37 

0. 


0.329037 

2.420463 

0. 

-0. 112034 

0.025378 


38 

0. 


0.436063 

2.875520 

0. 

-0. 159752 

0.035782 


39 

0. 


0. 

O. 

o. 

O. 

0. 


40 

0. 


0. 

O. 

0. 

0. 

O. 





FORCES CONTRIBUTED BY 

THE SUBSYSTEM 

MODE SHAPES 



POINT 

X 


Y 

Z 

THETA-X 

THETA-Y 

THETA-Z 


NUMBER 

POUNDS 


POUNDS 

POUNDS 

IN-LB 

IN-LB 

IN-LB 


1 

0. 


0. 

0. 

0. 

0. 

0. 


2 

0. 


18.762 

1.626 

0. 

- 18.877 

217.830 

r 


o. 


- 18.762 

- 1 .626 

O. 

-18.877 

217.830 



0. 


0. 

0. 

0. 

0. 

0. 



0. 


-295.807 

-33.684 

0. 

-2426.568 

-27209. 151 



0. 


907.485 

-42.419 

0. 

1239.421 

2256.436 





0. 

0. 

0. 

0. 

0. 





-1357.509 

-71.050 

0. 

-1700.477 

-31897.414 





843.077 

-48.460 

0. 

337.975 

-15327.313 
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36 

37 

38 

0. 

0. 

o. 

0. 

-8.365 

530.845 

0. 

-0.725 

-53.035 

0. 

0. 

0. 

0. 

-114.357 

-705.069 

0. 

1319.631 

-28050.409 


/ 

THE FOLLOWING VALUES 

ARE FOR THE 

TYPE 3 PHYSICAL CONNECTING ELEMENTS. 




FORCES 

ARE THOSE THAT THE ELEMENT 

EXERTS ON THE ROTOR OR CASE- 





RELATIVE 

DEAD 

I -END 

JEND 

FORCE IN Y DIRECTION 

FORCE IN Z DIRECTION 

FORCE 

ELEMENT 

DISPLACEMENT 

BAND CLEARANCE POINT 

POINT 

I END 

J END 

I END U END 

MAGNITUDE 

NUMBER 

INCHES 

INCHES INCHES NUMBER NUMBER 

POUNDS 

POUNDS 

POUNDS POUNDS 

POUNDS 

5 

0.0112 

0,0100 -0. 

0012 4 

1 

-1161.945 

1161.945 

209.344 -209.344 

1 180.653 

THE FORCES THAT THE 

TYPE 5 PHYSICAL CONNECTING ELEMENTS (UNCOUPLED POINT 




SPRING- 

DAMPER ELEMENTS) EXERT ON 

THE ENGINE COMPONENTS OR 

GROUND ARE- 








FORCE IN GIVEN DIRECTION 



ELEMENT 

POINT 

X 

Y 

Z 

THETA-Y 

THETA-Z 


NUMBER 

END NUMBER 

POUNDS 

POUNDS 

POUNDS 

IN-LB 

IN-LB 


1 

I 

2 

0. 

-11094.245 

1548.797 

0. 

0. 


1 

J 

39 

0. 

11094.245 

-1548.797 

0. 

0. 


2 

1 

3 

0. 

3352.860 

-916. 161 

0. 

0. 


2 

J 

40 

0. 

-3352.860 

916. 161 

0. 

0. 


3 

I 

5 

0. 

-8834.451 

1137.095 

0, 

0. 


3 

J 

2 

0. 

8834.451 

-1137.095 

0. 

0. 


4 

I 

6 

0. 

3141.491 

-834.069 

' 0. 

0. 


4 

J 

3 

0. 

-3141.491 

834.069 

0. 

0. 



' THE GYROSCOPIC FORCES ACTING ON THE ROTOR(S) ARE- 


POLAR MOMENT Y-AXIS Z-AXIS 
POINT ROTOR OF INERTIA MOMENT MOMENT 
NUMBER NUMBER LB-IN««2 IN-LB IN-LB 


4 1 184205. -5670.702 -27273.964 

36 1 184205. -5305.378 -22747.371 


SUMMARY OF UNBALANCE FORCES- 


BIRTH PHASE 

TIME POINT ROTOR MAGNITUDE ANGLE FORCE (LB.) FORCE (LB.) 

SECONDS NUMBER NUMBER GM-IN DEGREES Y-DIRECTION 2-DIRECTION 


O. 4 1 5000. 0. 2815.517 -4.809 


3 

r 

I 

K 

8 

K 


GENERALIZED 

COORDINATE 

NUMBER 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


GENERALIZED 
FORCE DUE 
TO APPLIED 
FORCES ONLY 


4.809 

-4.809 

-0.587 

-1 .321 

0.498 

-2815.517 

2815.517 

343.9 15 

773.704 

-291.716 

O. 







1 12 

O. 



13 

0. 



14 

0. 



15 

0. 



16 

0. 



GENERALIZED 




COORDINATE 

GENERALIZED 

GENERALIZED 

GENERALIZED 

NUMBER 

DISPLACEMENT 

VELOCITY 

FORCE 

1 

0.00573250 

-7.993238 

-2000.988 

2 

-0.00137468 

3.429242 

507.561 

3 

-0.00004335 

0.476910 

-38.959 

4 

0.00002635 

-0.066357 

504.065 

5 

-0.00000372 

-0.076811 

14.649 

6 

-0.02544950 

. -1.789783 

8927.608 

7 

0.01091575 

0.432514 

-4039.388 

8 

0.00150688 

0.021673 

1484.309 

9 

-0.00021180 

-0.007892 

-4058.038 

10 

-0.00023309 

0.025221 

-801.295 

11 

-0.00105548 

2.416374 

120.267 

12 

0.00406523 

-5.601686 

-185.691 

13 

-0.00000168 

-0.005691 

-7.721 

14 

0.00769519 

6.3^9839 

-886.480 

15 

-0.01783598 

-1.268885 

814.985 

16 

-0.00001940 

0.001117 

-91.077 



GENERALIZED 

GENERALIZED 

GENERALIZED 




WEIGHT 

STIFFNESS 

DAMPING VALUE 

GENERALIZED 



POUNDS 

LB/IN 

(LB-SEC)/IN 

ACCELERATION 



1.372E 

03 

0. 


0. 

-563.4668 



1.449E 

03 

6. 


0. 

135.37 14 



3.743E 

02 

1.076E 

06 

0. 

7.9380 



2.457E 

03 

1 .980E 

07 

0. 

-2.7821 



1 .766E 

02 

3 . 500E 

06 

o. 

60.5520 



1 .372E 

03 

0. 


o. 

2513.9631 



1.449E 

03 

0. 


0. 

-1077.3447 



3.743E 

02 

1.076E 

06 

0. 

-141.5476 



2.457E 

03 

1 .980E 

07 

0. 

21.3253 



1.766E 

02 

3 . 500E 

06 

0. 

31 .7929 



4 . 503E 

02 

0. 


0. 

103.2016 



1 .784E 

02 

0. 


0. 

-402. 1903 



1.787E 

02 

4 . 760E 

06 

9.891E 01 

1.8247 



4.503E 

02 

0. 


0. 

-760.6942 



1.784E 

02 

0. 


0. 

1765. 1854 



1.787E 

02 

4 . 760E 

06 

9.891E 01 

2.4907 
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OUTPUT PLOT FILE IS A BINARY SEQUENTIAL FILE 
LISTING OF AT LEAST PART OF OUTPUT PLOT FILE- 



) 

8 

1 

1 

1 8 1025 

1 0. 

3 0. 

0. 

0. 

0. 

0. 





4 

1 0. 

0. 

0. 





4 

3 0. 

0. 

0. 





37 

1 -50.000 

0- 

0. 





37 

3 -50.000 

0. 

0. 





38 

1 '-50.000 

b. 

0. 





38 

3 -50.000 

0. 

0. 





5 

5 5 







3 

2 1 







3 

2 3 







3 

5 1 







3 

5 3 







4 

3 1 







4 

3 3 







4 

6 1 







4 

6 3 







0. 

3000. 

0 0. 

0. 


0. 




0. 

0. 


0. 





0. 

o. 


0. 





0. 

0. 


0. 





0. 

0. 


0. 





0. 

0. 


0. 





0. 

0. 


0. 





0. 

0. 


0. 





0. 

0. 


0. 





3. 0.01000 0. 







0. 

0. 

0. 






0. 

0. 

0. 






0. 

0. 

0. 






0. 

0. 

0, 






0. 

0. 

0. 






0. 

0. 

0. 






0. 

0. 

0. 






0. 

0. 

0. 





0.0005000 3000. 

0 0. 

0. 

02500000 

0. 




0 . 00000094 

0.007831 


0. 





0.00002663 

0. 167883 


0. 





0.00001077 

0.064133 


0. 





0.00020520 

0.804002 


0. 





-0 . OOOOOOOO 

0.000040 


0. 118 





0.00000018 

0.003689 


3. 194 





-0.00000132 

-0.006846 


0.222 





-0.00002222 

-0.053297 


-7.959 





0.00018 0.00982 0. 







15.004 

15.004 

15. 

004 





230.258 

230.258 

230. 

258 





- 15.004 

- 15.004 

- 15. 

004 





-230. 258 

-230.258 

-230. 

258 





-1.264 

-1.264 

- 1 . 

264 





-21.225 

-21.225 

-21 . 

225 





1.264 

1 .264 

1 . 

264 





2 1 . 225 

21 . 225 

21 . 

225 




II 0.0010000 3000 

0 0. 

0. 

05000000 

0. 




0.00001672 

0.066912 


0. 





0,00022502 

0.657962 


0. 





0.00008429 

0.247373 


0. 









( 

0.00077849 
0 . 00000079 
0-00001721 

1.461820 

0.005175 

0.088271 

0. 

1.768 

20.892 

( 

1 I 

( 


-0-00000470 

-0.001465 

3.634 




0.00000761 

0.229270 

-27.652 




0.00056 

0.00944 0. 





72.870 

72.870 

72.870 




557,798 

557,798 

557.798 




-72.870 

-72.870 

-72.870 




-557.798 

-557.798 

-557.798 




-4.094 

-4.094 

-4.094 




- 12.203 

-12.203 

-12.203 




4.094 

4.094 

4,094 




12.203 

12.203 

12.203 




0.0015000 

3000.0 O. 

0.07500000 

0. 



0.00008096 

0.201962 

0. 




0.00069202 

1. 197454 

0. 




0.00027534 

0.531994 

0. 




0.00164013 

1.963227 

0. 




0.00000949 

0.035764 

6.474 




0.0001 1922 

0.346955 

41.671 




0.00000251 

0.032435 

19.325 




0-00023341 

0.657456 

-45.672 




0.00097 

0.00903 O. 





175.670 

175.670 

175,670 


1 


869.411 

869.411 

869.411 




- 175.670 

-175.670 

-175.670 




-869.411 

-869.411 

-869.411 




-1.333 

-1.333 

-1.333 




74.811 

74.811 

74.811 




1,333 

1.333 

1.333 




-74.811 

-74.811 

-74.811 




0.0020000 

3000.0 0. 

0. fOOOOOOO 

0. 



0.00023031 

0.405774 

0- 




0.00140063 

1.612629 

0. 




0.00063139 

0.905093 

0. 




0.00272247 

2.349578 

0. 




0.00004442 

0. 111813 

12.646 




0.00038136 

0.703936 

44.688 • 




0-00002594 

0-058152 

60.431 




0.00061597 

0.816088 

-31.739 




0.00138 

0.00862 O. 





322.392 

322.392 

322.392 




1135.017 

1135.017 

1135.017 




-322.392 

-322.392 

-322.392 




-1135.017 

-1135.017 

-1135.017 




10.818 

10.818 

10.818 




172.825 

172.825 

172.825 




- 10.818 

- 10.818 

-10.818 




-172.825 

-172.825 

-172.825 




0.0025000 

3000.0 0. 

0. 12499999 

0. 



0-00049786 

0.675256 

0. 




0.00227675 
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7.2 COMPUTED RESULTS 


For the plotted results that are presented, the following computing 
strategy was used for constant speed running. 

at = 50 microseconds 

Number of time integration points = 10,240 

Time frame ® 10,240 x 50 x 10"^ * 0.512 second. 

Figure 55 shows the TETRA computed displacement-time history at the mid- 
dle of the rotor (point 38) in the vertical direction corresponding to a sud- 
den 100 gm-in. fan unbalance at a constant 12,000 rpm speed. The number of 
computed values were decimated so that 1,024 points are shown and the time 
increment between points is equal to 50 x 10“5 seconds. The 12,000 rpm speed 
corresponds to super critical speed operation relative to a rotor dominated 
mode computed at 9,908 rpm for the total system. Overshoot is clearly in evi- 
dence for the transient response shown in Figure 55. The Fast Fourier Trans- 
form of this time history response that is shown in Figure 56 indicates that 
four modes are contributing to the response. The peak at 200 Hz (12,000 rpm) 
is associated with the driving force. The VAST predicted mode shapes for the 
total system sketched in Figure 57 show that the modes at 3617 rpm, 9908 rpm, 
13,983 rpm, and 26,473 rpm should be responsive at the mlddTe of the rotor for 
fan unbalance. This indicates that TETRA has correctly synthesized the modal 
and physical data to predict the time transient response for the total system. 
Figure 58 shows the displacement-time history at the middle of the rotor in 
the vertical direction corresponding to a sudden 100 gm-in. fan unbalance at 
9908 rpm critical speed operation. In this case, there is no overshoot and 
the transient response builds up to the steady-state level. This response be- 
havior is characteristic of operation at a critical speed. 

The results presented in Figures 55 and 58 do not include gyroscopic 
stiffening. Figure 59 shows the off-resonant response in two planes at the 
fan that reflects the effects of gyroscopic stiffening at both the fan and 
turbine. Figure 60 shows orbit plots of the data presented in Figure 59. 
Inspection of Figure 60 shows that the initial response is a noncircular 
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Figure 55. Demonstrator Model - 100 gm-in. Sudden Unbalance at the Fan at 12.000 rpm 
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Figure 56. Demonstrator Model - 100 g^-in. Sudden Unbalance at the Fan at 12,000 rpm. 
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Figure 58. Demonstrator Model - 100 gm-in. Sudden I 
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Figure 60. Orbit Plots at the Fan for NASA Demonstrator Model with Gyro 
100 gm-in. Sudden Fan Unbalance at 3000 rpm. 


whirl. This noncircular whirl transitions into a circular whirl after steady- 
state conditions are reached. Figures 51 through 55 show response plots at 
the fan and case that include both the effects of gyroscopic stiffening and 
the rub load path. Figures 67 through 70 show the steady-state frequency 
response for the total system with and without gyroscopic stiffening at the 
fan and the turbine locations. Comparison of the time response, once steady- 
state conditions have been reached, with the frequency response shows good 
agreement . 

Figures 71 through 76 show the results of analyses for a 2,000 rpm/sec 
accel rate that reflects use of the TETRA restart option. The restart option 
allows the utilization of the results obtained from a previous analysis to 
continue the analysis out to future time without the introduction of pseudo 
transients. For several different times (those for which output was printed) 
the TETRA program writes to a restart file data including values for each of 
the generalized coordinates for the current time step, for one time step 
earlier, and for two time steps earlier. The user chooses one of the times 
for which output was printed on the original run (usually the last such time) 
as the restart time for the new run. The data for that time from the restart 
file is then utilised to provide the initial data needed to continue the anal- 
ysis. The process can be repeated any number of times. 

The computing strategy used to generate the response data shown in Fig- 
ures 71 through 76 was as follows. The maximum speed for the accel was 
selected as 5,000 rpm so as to encompass the two gyro stiffened critical 
speeds located at 3292 rpm and 3624 rpm for steady-state operation. Modal 
Q-factor values of 15 were used for each of the casing subsystem modes and 
the rotor subsystem modes were undamped. Physical damping for the connecting 
springs based on a Q-factor equal to 15 and a 60.4 Hz frequency, corresponding 
to tbe 3624 rpm mode, was used. 

The analysis was accompanished with three restart segments for the fol- 
lowing cases. 

a. Large fan-case radial clearance (200 mils) ' 

b. Small fan-case radial clearance (10 mils) 
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Figure 63. Rotor and Case Orbit Plot for the NASA Demonstrator Model for 5000 Gm~In. 
Sudden Fan Unbalance at 3000 RFH (No Rub) • 
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Shown to time 0.05 second. 10 mil radial displacement dead band and 
1 X 10^ Ib/in. rub spring at the fan rotor-case. 
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Figure 64. Rotor and Case Orbit Plot for the NASA Demonstrator Model for 5000 Gm-In. 
Sudden Fan Unbalance at 3000 RPM (With Rub). 
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Shown to time 0,512 second. Radial displacement dead band 
exceeds the fan rotor-case relative dispiacemenc . 
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Figure 65. Response in the Vertical Direction at the Case and Fan 
Rotor for th^ NASA Demonstrator Model for 5000 Gm-In. 
Sudden Fan Unbalance at 3624 RPM (No Rub) . 
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Shown to time 0.512 second. 10 mil radial displacement dead band and 1 x 10 
lb/ in. rub spring at the fan rotor-case. 
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Figure 66. Response In the Vertical Direction at the Case and Fan 
Rotor for the NASA Demonstrator Model for 5000 Gm-ln. 
Sudden Fan Unbalance at 3624 RPM (With Rub) . 
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Figure 67, NASA Demonstrator - Steady-State Frequency Response, 100 gm-in. Pan (No Gyro) 
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Figure 69. NASA Demonstrator - Steady-State Frequency Response, 100 gm-in. Fan (Gyro a(: 

Fan and Turbine). 
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Figure 70, NASA Demonstrator - Steady-State Frequency Response, 100 gm-in. Fan (Gyro at 
Fan and Turbine)* 
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Figure 71. Response in the Vertical Direction of the Case and Fan 
Rotor for the NASA Demonstrator Model for 1000 to 3000 
RPM Accel Segment. Radial Displacement Dead Band Exceeds 
the Fan-Case Relative Displacement (No Rub) . 
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Figure 72. Response in the Vertical Direction at the Case and Fan 
Rotor for the NASA. Demonstrator Model for 3000 to 5000 
RPM Accel Segment. Radial Displacement Dead Band Exceeds 
the Fan-Case Relative Displacement (No Rub) . 
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Figure 73. Response in the Vertical Direction at the Case and Fan 
Rotor for NASA Demonstrator Model for 1000 to 3000 RPM 
Accel Segment. 10 Mil Radial Displacement Dead Band and 
1 X 10 Lb/In. Rub Spring at the Fan Rotor-Case (With Rub). 
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Figure 74. Response in the Vertical Direction at the Case and Fan 
Rotor for the NASA Demonstrator Model for 3000 to 5000 
RPM Accel Segment. 10 Mil Radial Displacement Dead Band 
and 1 X 106 Lb/ In. Rub Spring at the Fan Rotor-Case 
(With Rub). 
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Force Exerted by the Forward Frame/Bearing (Spring 3) 
on the Rotor at Point 5 in the Vertical Direction for 
the NASA Demonstrator Model for 3000 to 5000 RPM Accel 
Segment. Radial Displacement Dead Band Exceeds the Fan 
Rotor-Case Relative Displacement (No Rub) . 
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Figure 76. Force Exerted by the Forward Frame/Bearing (Spring 3) 
on the Rotor at Point 5 in the Vertical Direction for 
the NASA Demonstrator Model for 3000 to 5000 RPM Accel 
Segment. 10 Mil Radial Displacement Dead Band and 
1 X 10^ Lb/In. Rub Spring at the Fan Rotor Case (With Rub) . 
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Initial constant speed runs were made for 0.520 second of simulated 
time with the rotor speed set at 1,000 rpm and with 100 gm-in. fan unbalance 
to establish steady^state conditions. For the 2,000 rpm/sec accel rate, the 
time required to accel from 1,000 to 5,000 rpm is 


(5000-1000)rpm 
2000 rpm/sec 


2 seconds. 


With an integration time step of At = 50 microseconds, the number of time 
steps required is equal to 


2 seconds 


50 X 10 


-6 seconds 
time step 


40,000 time 


steps . 


For both cases a and b, two restart segments of 20,000 integration points each 
were used to compute the response for the 2 second accel from 1,000 to 
5,000 rpm. For plotting purposes, the computed data were decimated so that 
each of the restart segment plots represents 2,000 points with a time incre- 
ment between points equal to 50 x 10*5 seconds. This means that there are at 
least (60/5000) x 1/(50 x 10*^) =» 24 plotted points per cycle at the highest rotor 
speed (5,000 rpm). Although the plots are labeled 0 to 1.0 second for the 
accel from 1,000 to 3,000 rpm, the 0 time actually corresponds to 0.52 sec- 
ond. At a time of 0.6 second (labeled as 0.6 - 0.52 « 0.08 second), a large 
unbalance increment of 5000 gm-in, is added to represent sudden unbalance. 

It will be noted that the apparent critical speeds for the large clearance 
case (i.e., no rub) occur at higher rotor speeds than the steady-state values. 

The 3580 rpm speed corresponds to the 3292 rpm speed and the 3960 rpm corre- 
sponds to the 3624 rpm speed. The shifting of the peak response speeds is 
caused by the accel rate and this phenomenon is discussed in the literature^^^ 

The fan rub which is present for the small clearance case causes a stiffening 
effect that shifts the peak response to a higher speed. 


(1) Vibration During Acceleration Through a Critical Speed, F.M. 
Lewis: Trans. ASME, Vol. 54, pp253-261, 1932. 
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